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FOREWORD 


The  Applied  Physics  Laboratory  (APL)  of  The  Johns  Hopkins  University, 
with  field  activities  in  Florida,  Utah,  New  Mexico,  Virginia,  New  Jersey, 
California,  and  Heidelberg,  Germany,  is  located  in  Howard  County,  Mary- 
land, midway  between  Baltimore  and  Washington.  Its  work  is  carried  out 
under  a contractual  agreement  between  the  University  and  the  Government, 
with  about  80Cfi  of  the  effort  devoted  to  the  Department  of  the  Navy.  APL 
employs  a staff  of  about  2400,  more  than  1200  of  whom  arc  professional 
scientists  and  engineers.  Their  ideas  are  implemented  and  extended  through 
a network  of  associate  contractors  from  coast  to  coast. 

The  primary  mission  of  APL  is  to  serve  the  nation  through  the  application 
of  advanced  science  and  technology  to  the  solution  of  problems  of  importance 
in  both  the  defense  and  the  civil  sectors.  APL  conducts  board  programs  in 
basic  and  applied  research,  exploratory  and  advanced  development,  com- 
ponent engineering,  systems  engineering  and  integration,  and  test  and  evalua- 
tion of  operating  systems. 

The  Laboratory  was  established  in  1942  at  the  urgent  request  of  the  Office 
of  Scientific  Research  and  Development  to  demonstrate  that  science  and  tech- 
nology could  advance  national  security  more  effectively  through  an  active 
partnership  of  military  and  civilian  technologists  than  through  the  traditional 
buyer -supplier  relationship.  The  concept  was  validated  in  less  than  a year 
when  APL  developed  the  first  practical  VT  (variable  time)  fuze.  This  initial 
product  proved  to  be  a major  contributor  to  the  Allied  victory  in  World 
War  II. 

After  the  war,  at  the  request  of  the  Secretary  of  the  Navy,  the  University 
trustees  agreed  to  continue  the  Laboratory  as  an  important  national  security 
resource.  APL  quickly  became  a leader  in  guided  missile  and  space  technology 
as  well  as  a contributor  to  fundamental  knowledge  in  such  areas  as  high- 
altitude  physics,  spectroscopy,  flame  propagation,  the  origin  and  evolution  of 
the  universe,  supersonic  and  hypersonic  aerodynamics,  chemical  kinetics,  and 
magnetic  resonance. 

In  addition  to  major  contributions  to  national  defense,  effort  is  increasingly 
devoted  to  the  solution  of  pressing  problems  in  the  civilian  sector.  A few 
examples  of  the  areas  to  which  attention  has  been  devoted  in  recent  years  are 
biomedical  research  and  engineering;  transportation;  fire  research;  pollution 
of  the  biosphere;  energy,  including  ocean  thermal,  geothermal,  and  com- 
munity storage;  navigation  satellite  technology;  leak  detection  in  natural  gas 
distribution  lines;  computer  technology;  and  advanced  education. 


In  its  role  as  a University  laboratory,  APL  performs  research  to  advance 
the  basic  sciences  underlying  its  technological  development  programs.  In  this 
area,  APL  has  significantly  contributed  to  Johns  Hopkins’  100-year-old 
reputation  for  leadership  and  originality.  Through  unique  applications  of 
systems  engineering,  science,  and  technology  to  the  needs  of  society,  APL  has 
enhanced  the  University's  tradition  of  excellence  while  gaining  worldwide 
recognition  of  its  own. 

& This  report  contains  a representative  sample  of  APL  technical  accomplish- 
mentw primarily  in  the  civil  area.  All  results  of  these  programs  arc  fully 


available  to  the  public. 
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INTRODUCTION 


The  Laboratory's  involvement  in  space  programs  began  in  the  postwar  years 
when  V-2  and  Aerobee  rockets  carried  Geiger  tubes,  magnetometers,  and 
optical  spectrometers  high  above  the  earth's  surface.  The  flights  provided  the 
first  high-altitude  measurements  of  cosmic  rays,  the  geomagnetic  field,  and 
atmospheric  constituents  such  as  ozone,  and  were  conducted  by  pioneers  James 
A.  Van  Allen,  John  J.  Hopfield,  and  S.  Fred  Singer  (who  were  then  APL 
staff  members).  From  these  distinguished  beginnings,  the  APL  record  of 
accomplishments  includes  the  conception,  design,  and  development  of  the 
Satellite  Doppler  Navigational  System  for  the  Navy. 

The  satellite  activities  spawned  a multitude  of  "firsts”  by  APL,  including 
development  of  the  first  gravity-gradient  satellite,  the  first  photographs  (in 
black  and  white  and  in  color)  of  the  entire  earth  from  synchronous  satellites, 
the  first  solid-state  particle  detectors  flown  on  a satellite,  and  the  first  ex- 
tremely accurate  measurements  of  the  geomagnetic  field. 

The  Space  Development  Department  of  APL  has  continued  to  maintain 
a strong  balance  between  the  development  of  spacecraft  systems  and  hard- 
ware and  basic  research  activities  directed  toward  an  understanding  of  the 
earth  and  the  solar  system.  This  applied  and  research  expertise  has  been 
manifested  in  the  past  year  with  the  successful  launch  of  the  APL-built 
TRANSAT  satellite  (designed  to  test  elements  of  the  U.S.  Navy's  Satellite 
Missile  Tracking  (SATRACK)  System)  and  the  launches  of  the  Voyagers  1 
and  2 spacecraft,  bound  for  Jupiter  and  Saturn,  carrying  APL  scientific 
experiments  to  measure  charged  particles  in  the  solar  system. 

Research  programs  included  scholarly  investigations  of  ancient  astronomical 
records,  which  led  to  the  suggestion  that  Ptolemy  was  a fraud;  the  discovery 
of  plasma  acceleration  regions  behind  the  earth  that  can  generate  charged 
particles  with  energies  up  to  hundreds  of  thousands  of  electron  volts-  Investi- 
gation of  the  interaction  between  the  interplanetary  magnetic  field  and  the 
geomagnetic  field  (a  possible  link  between  solar  activity  and  terrestrial 
weather) ; and  the  discovery  of  the  fact  that  Jupiter  is  a prominent  source  of 


energetic  particles  in  the  vicinity  of  the  earth  (they  previously  were  thought 
to  originate  in  the  sun). 

While  contributing  to  a fundamental  understanding  of  the  physical  and 
chemical  processes  that  occur  in  the  vast  region  from  the  earth's  upper  atmo- 
sphere to  the  distant  planets,  API.  staff  members  have  advanced  the  state  of 
the  art  of  a wide  variety  of  spacecraft  instruments,  from  complete  satellites 
(including  attitude  control  and  microprocessor  systems)  to  complex  scien- 
tific detectors.  For  example,  the  energetic  particle  detectors  on  board  the 
Voyager  spacecraft  (which  were  designed,  tested,  and  constructed  at  APL) 
are  more  complicated  than  the  entire  Mariner  IV  spacecraft  that  traveled  to 
Mars. 

NASA  has  selected  scientists  and  engineers  from  APL  to  participate  in  a 
record  number  of  interplanetary  missions.  (These  missions  include  Voyagers 
1 and  2,  which  will  rendezvous  with  Jupiter  in  1979;  Galileo,  which  is 
scheduled  to  orbit  Jupiter;  and  Solar  Polar,  which  will  orbit  over  the  poles 
of  the  sun.)  Support  for  the  joint  APL 'Max- Planch  Institute  Active  Mag- 
netosphere Particle  Tracer  Experiment  (AMPTF)  was  confirmed  this  year 
by  NASA  and  the  Federal  Republic  of  Germany.  Hie  purpose  of  the  experi- 
ment is  to  create  an  artificial  ion  cloud  outside  the  earth's  magnetosphere  in 
order  to  investigate  the  mechanisms  responsible  for  the  formation  of  the 
Van  Allen  radiation  belts  (an  appropriate  study  linked  with  the  pioneering 
activities  of  Van  Allen  at  APL). 

In  the  incredibly  short  period  of  a few  weeks,  an  energetic  electron  de- 
tector was  designed,  assembled,  and  tested  at  APL  and  successfully  launched 
on  the  NASA  International  Ultraviolet  Explorer  (IUE)  spacecraft  this  year. 
APL  is  presently  building  for  NASA  an  earth-orbiting  satellite  called 
MAGSAT,  to  measure  the  earth's  magnetic  field  with  an  accuracy  not  pre- 
viously available  on  a worldwide  basis.  APL  has  also  built  principal  instru- 
ments for  the  SEASAT  satellite  that  will  measure  the  general  characteristics 
of  the  oceans  on  a global  scale  in  all  weather  conditions. 
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TRANSAT  TRANSLATOR 


A translator  has  been  designed,  using  microelec- 
tronic packaging  techniques,  to  translate  two  UHF 
input  spectra  and  a tracking  aid  ("pilot  carrier")  to  a 
wideband  spectrum  at  S band  tor  retransmission  to 
ground  receivers.  The  translator  has  been  incorporated 
into  a Navy  Navigation  Satellite.  The  combination, 
known  as  TRANSAT,  will  provide  a means  tor  early 
checkout  ot  surface  stations  and  data-processing 
facilities  of  the  Satellite  Missile  Tracking  (SATRACK) 
System  now  being  developed  and  evaluated  at  APL. 


BACKGROUND 

Designed  to  work  with  a constellation  of  Global 
Positioning  Satellites  (GPS)  that  transmit  signals  at 
1575.42  MHz,  the  missile-borne  translator  tracks  the 
signals  by  range  difference  and  doppler  methods.  The 
translator  also  receives  signals  at  393.855  MHz  from 
a ground  transmitter.  This  lower-frequency  signal  is 
used  to  measure  the  ionospheric  electron  density  as 
the  missile  transits  the  ionosphere;  it  is  also  used 
as  a range  safety  signal.  The  spectra  of  these  two  sig- 
nals, along  with  a pilot  carrier  signal  derived  from 
the  translator's  local  oscillator,  are  translated  to 
2273.5  MHz  (S  band)  and  transmitted  to  tracking 
ground  stations.  Figure  1 illustrates  the  system. 

By  orbiting  the  translator  as  an  add-on  to  a satel- 
lite, a permanent  test  vehicle  will  be  available  to 
simulate  a Trident  missile  flight  each  orbit.  Since  the 
satellite  chosen  in  this  case  is  a Navy  Navigation 


'T  GPS 
ll575.42  MHz 


Satellite,  it  can  be  commanded  to  its  normal  naviga- 
tion mode  when  the  translator  is  not  in  use. 

For  redundancy,  two  translators  were  designed, 
space  qualified,  and  incorporated  into  a penthouse 
structure  on  the  satellite.  Either  may  be  activated  by 
ground  command. 


DISCUSSION 

Under  ground  station  control,  a translator  is 
selected  and  commanded  on.  The  2273.5-MHz  sig- 
nal is  transmitted  through  four  antennas  mounted  on 
four  booms  equally  spaced  around  the  diameter  of 
the  satellite.  After  15  minutes,  the  translator  will 
automatically  turn  off.  However,  an  off  command 
from  the  ground  station  will  turn  it  off  regardless 
of  where  it  is  in  the  on  cycle. 

The  translator  is  shown  in  Fig.  2,  and  Fig.  3 is  a 
block  diagram  indicating  the  basic  operation. 

Design  ranges  for  mpu*  signals  for  both  the 
393.855-  and  1575.42-MHz  channels  (channels  1 
and  2,  respectively)  are  between  —130  and  — <50 
dBm.  The  channel  1 signal  is  amplified  and  mixed 
with  a local  oscillator  (LO)  derived  from  a 5-MHz 
source  to  provide  25.855  MHz  at  the  summing  am- 
plifier. Likewise,  channel  2 is  heterodyned  to  28.42 
MHz.  In  order  to  aid  the  tracking  of  these  signals,  a 
pilot  carrier  is  also  added  at  the  summing  amplifier. 
This  25-MHz  signal  is  coherent  with  the  other  LO 
signals  since  it  is  derived  from  the  same  5-MHz 
temperature-compensated  crystal  oscillator. 


Fig.  1 The  SATRACK  system. 


Fig.  2 The  TRANSAT  translator. 
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Fig.  3 Block  diagram  of  the  translator. 
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Fig.  4 Output  spectrum  of  the  translator. 


being  translated  to  S band  with  the  spectral  pattern 
showt.  in  Fig.  4 and  a power  level  of  —3  dBm. 

Frequency  synthesis  by  phase-locked  loop  tech- 
nique was  used  to  derive  the  mixing,  or  LO,  signals. 
This  allowed  the  use  of  off-the-shelf  integrated  cir- 
cuits with  resultant  savings  in  power  and  space. 

The  output  power  requirements  of  the  translator 
were  5 W (37  dBm)  with  30-dB  intermodulation 
distortion  products.  Intermodulation  distortion  of  this 
magnitude  at  high  power  levels  is  difficult  to  obtain 
when  weight  and  volume  are  at  a premium.  Since 
a gain  of  40  dB  (—3  to  37  dBm)  was  required, 
the  output  power  amplifier  was  divided  into  two  sec- 
tions: a low-level  section , —3  to  20  dBm,  consisting 
of  four  stages,  and  a high-level  section,  20  to  37 
dBm,  consisting  of  a hybrid  two-stage  amplifier.  It 
was  difficult  to  meet  the  intermodulation  requirement 
in  the  high-level  stages  and  some  new  development 
was  required.  The  final  device  consisted  of  three  RF 
transistors  in  a series/parallel  hybrid  configuration. 
Each  transistor  consisted  cr  a ’e-cell  pellet  from 
adjacent  areas  of  the  original  silicon  wafer.  The 
typical  power  output  from  toe  transistors  was  2.9  to 
3.2  W at  30-dB  intermodulation  distortion.  Tempera- 
ture-compensated bias  networks  provided  linearized 
class  B operation  for  the  common-base  configuration. 

To  provide  a means  for  evaluating  overall  transla- 
tor performance,  several  power  levels  and  tempera- 
tures were  telemetered  from  the  spacecraft  to  ground 
tracking  stations.  The  parameters  included  power  sup- 
ply current,  input  signal  levels,  temperatures,  and 
output  power. 

Because  of  the  space  and  weight  constraints  of  the 
satellite,  microelectronic  packaging  techniques  were 
used  to  the  fullest.  Modular  construction  consisted  of 
several  Kovar  chemical  milled  chassis.  Active  circuits 
utilized  microstrip  substrates  and  miniature  chip  com- 
ponents. Therefore,  both  electrical  isolation  and  the 
mechanical  integrity  dictated  by  the  spacecraft  launch 
environment  were  achieved. 

TRANSAT  has  passed  all  electrical  and  environ- 
mental tests  at  APL  under  vibration  and  thermal- 
vacuum  conditions.  It  was  shipped  to  Van  '•tAurg 
Air  Force  Base  for  launch  on  a Scout  D rocket  sched- 
uled for  early  October  1977. 


Subsequent  mixing  and  amplification  results  in  the 
composite  spectrum  at  the  summing  amplifier  output 
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RF  SYSTEM  FOR  THE  TIME-TRANSFER 
RECEIVER  OF  THE  NAVIGATION 
TECHNOLOGY  SATELLITE 


An  RF  system  has  been  designed  and  developed 
tor  a timing  receiver  that  NASA  will  use  to  obtain 
global  submicrosecond  timing.  The  compact,  low-cost 
receiver  passively  receives  sidetone  ranging  signals 
trom  the  Navigation  Technology  Satellite  (NTS)  and, 
under  the  control  of  a microprocessor,  is  capable  of  a 
precision  of  better  than  0.5  gs.  It  is  potentially  useful 
to  anyone  with  a need  for  an  inexpensive  time-transfer 
capability. 


BACKGROUND 

The  NASA  Goddard  Space  Flight  Center  (GSFC) 
Laser  Tracking  Network,  currently  being  implemented, 
requires  worldwide  clock  synchronization  of  ± 1 ixs. 
Recognizing  that  this  precision  is  not  achievable  by 
conventional  time-transfer  systems  such  as  Loran  C or 
by  extended  portable  clock  trips,  GSFC  initiated  the 
development  of  a single-frequency  time-transfer  re- 
ceiver to  provide  suitable  remote-dock  synchroniza- 
tion data.  Using  ranging  signals  from  the  NTS,  the 
receiver  will  provide  NASA  with  timing  data  into  the 
early  1980's.  Among  the  non-NASA  users  arc  other 
United  States  Governmental  agencies  and  several  for- 
eign countries. 

The  NTS  receiver  is  the  result  of  a joint  develop- 
ment effort  by  APL  and  the  Naval  Surface  Weapons 
Center  (NSWC),  Dahlgren,  Virginia.  NSWC  was 
responsible  for  the  digital  and  processor  portions  and 
the  mainframe  chassis.  The  RF  and  analog  portions 
were  designed  and  developed  by  APL. 


DISCUSSION 

The  NTS  receiver  (Fig.  1)  operates  automatically 
under  the  control  of  an  Intel  8080  microprocessor. 
Using  ranging  sidetones  transmitted  by  the  spacecraft, 
it  provides  an  "observed"  measurement  of  the  space- 
craft-to-receivcr  range  in  time  units,  once  each  minute 
that  the  spacecraft  is  in  view. 

To  convert  the  range  output  by  the  receiver  to 
useful  clock  error  measurements,  further  non-real-time 
processing  is  required.  Conceptually,  the  observed 
range  may  be  compared  with  a theoretical  range  based 
on  known  ground  and  spacecraft  positions  to  yield  an 
error  that  is  relatable  to  the  real  difference  between  a 
spacecraft  clock  and  a ground  station  clock.  In  actual 
practice,  however,  similar  measurements  are  made  at 
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two  sites,  a reference  site  with  a "master"  clock  and 
a remote  site  whose  relative  clock  time  is  in  question. 
If  the  two  observations  are  made  as  closely  in  time  as 
possible,  the  spacecraft  clock  may  effectively  be  elim- 
inated and  the  comparison  made  directly  between  mas- 
ter and  remote  clinks. 

The  receiver  comprises  three  basic  systems  (Fig.  2). 
The  RF  system  uses  a carrier  tracking  phase-locked 
loop  as  a reference  to  remove  carrier  doppler  fre- 
quency from  the  received  ranging  sidetones.  Signals 
generated  by  the  digital  processor  are  then  used  to 
convert  the  sequentially  received  sidetones  to  a com- 
mon 30-kHz  TTL  square  waveform.  This  signal  con- 
tains the  phase  information  carried  by  the  successive 
sidetones  and  extracted  by  the  digital  processor.  Using 
the  phase  information,  the  microprocessor  system  com- 
putes the  observed  spacecraft  range  at  a given  instant 
of  GMT,  and  this  range  is  distributed  to  a front-panel 
display,  a Teletype  output,  and  a computer-compatible 
output. 

Figure  3 is  a simplified  block  diagram  of  the  RF 
system.  The  signals  from  the  antenna  terminals  con- 
sist of  the  carrier  at  335.355  MHz,  a reference  tone 
at  3.35.325  MHz,  and  ten  sequentially  occurring  rang- 
ing tones  from  335.324  900  through  328.925  000 
MHz.  Each  signal  may  contain  a doppler  frequency 
offset.  In  the  rf  section,  the  band  of  frequencies  is 
amplified  and  converted  to  an  intermediate  frequency 
(IF)  in  which  the  tracked  component  doppler  fre- 
quency is  reduced  by  nine-tenths. 

In  the  if  section,  a 9-MHz  bandpass  filter  estab- 
lishes the  if  bandwidth,  and  automatic  gain  control 
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Fig.  2 Block  diagram  of  the  NTS  time-transfer  re- 
ceiver. 


(AGC)  is  applied  to  the  composite  signal.  The  track- 
ing loop,  which  operates  with  a 30-Hz  (double 
sided)  bandwidth,  supplies  a frequency  reference  to 
the  sidetone  mixer  and  selector,  where  carrier  doppler 
frequency  is  removed  from  the  sidetones,  which  are 
shifted  to  near-haseband  frequencies.  The  tracking 
loop  also  supplies  AGC  voltage  by  using  quadrature 
references  from  the  digital  system. 

Synthesized  sidetones  generated  by  the  digital  sys- 
tem are  used  in  a combination  of  double-  and  single- 
sideband techniques  by  the  sidetone  extractor  to  con- 
vert the  received  sidetones  to  30  kHz.  The  sidetones 
are  filtered  through  a 100-Hz  bandpass  filter  and 
detected  by  a threshold-crossing  TTI.  converter.  To 
permit  checkout  of  the  digital  processor,  a circuit  for 
generating  calibration  sidetones  has  been  included. 

The  RF  system  is  capable  of  operating  in  two 
modes:  in  one,  the  continuous  carrier  is  tracked;  in 
the  other,  the  system  tracks  the  reference  tone,  which 
is  present  only  during  the  burst  of  ranging  sidetones. 
In  either  of  these  modes,  the  tracking  loop  voltage- 
controlled  crystal  oscillator  (VCXO)  may  be  tuned 
manually  or  automatically  by  the  microprocessor  for 
rapid  acquisition. 

The  receiver  has  a system  noise  temperature  of  ap- 
proximately 600  K.  This  permits  a tracking  loop 
tracking  threshold  sensitivity  (about  + 3 dB  signal-to- 
noise  ratio)  of  about  —153  dBm.  However,  because 
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Fig.  3 RF  system  of  NTS  time-transfer  receiver. 


of  the  wider  sidetone  filter,  the  data  signal-to-noise 
ratio  is  0 dB  at  about  — 151  dBm. 

Several  months  of  data  have  been  taken  by  NTS  re- 
ceivers over  the  past  year.  The  data  shown  in  Fig.  4 
are  representative  of  the  receiver’s  performance.  The 
figure  shows  the  clock  error  measured  between  a 
master  clock  at  the  NRL  tracking  station  at  Chesa- 
peake Beach,  Maryland,  and  a remote  clock  at  the 
NASA  tracking  station  at  Kennedy  Space  Center  in 
Florida.  The  data  were  taken  between  days  2.33  and 
268  of  1977.  Data  points  show  the  determined  error 
in  the  NASA  clock  for  the  minute  during  each  day's 
pass  that  the  spacecraft  was  at  the  minimum  range 
from  the  receiver.  At  minimum  range,  the  errors 
arising  from  spacecraft  motion  and  atmospheric  effects 
are  minimum.  No  corrections  have  been  made  for 
these  effects,  nor  have  other  corrections  been  made  to 
provide  greater  accuracy.  However,  the  rms  value  of 
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Fig.  4 Representative  dock  error  data  using  NTS  time- 
transfer  receiver. 


the  data  about  a best-fit  curve  is  approximately  100  ns. 
Similar  data  taken  over  shorter  periods  processed  with 
updated  spacecraft  trajectory  predictions  have  yielded 
rms  values  about  one-half  that  amount. 
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AN  AUTOMATIC  GROUND  STATION  FOR 
COLLECTING  TELEMETRY  DATA  FROM 
THE  TRIAD  SATELLITE 


The  Tried  satellite  has  been  a means  of  collecting 
data  on  the  earth's  magnetic  field  for  more  than  five 
years.  Its  circular  polar  orbit  allows  nearly  continuous 
high-resolution  observations  to  be  made  over  both 
polar  regions  of  the  earth.  The  data  acquired  over 
these  regions  have  led  to  the  discovery  of  massive 
electric  currents  that  flow  along  the  earth's  magnetic 
field  lines  into  and  away  from  the  auroral  regions.  A 
completely  automatic  ground  station  has  been  de- 
signed, constructed,  and  tested  for  the  unattended 
co'lection  of  telemetry  data  from  Triad.  This  system 
represents  a third-generation  improvement  over  sta- 
tions now  operating  in  Resolute,  Canada,  and  McMurdo, 
Antarctica,  to  collect  Triad  data  in  the  North  and 
South  Polar  regions. 


BACKGROUND 

The  Triad  satellite  was  launched  into  a circular 
polar  orbit  at  800  km  altitude  in  September  1972  to 
test  several  techniques  for  improving  the  U.S.  Navy 
Navigation  Satellites  (Ref.  1).  A principal  component 
is  a highly  precise  magnetometer  that  measures  the 
geomagnetic  field  along  the  three  principal  axes  of  the 
spacecraft  to  determine  attitude.  After  the  navigation 
improvement  tests  were  performed,  the  satellite  was 
devoted  to  collecting  magnetic  field  data  in  the  polar 
regions  for  auroral  studies  (Ref.  2).  The  effort,  sup- 
ported by  the  Office  of  Naval  Research  and  the  Na- 
tional Science  Foundation,  has  required  the  use  of 
telemetry  collection  stations  in  remote  parts  of  the 
earth.  In  order  to  operate  this  equipment  and  still 


remain  within  the  limited  budget,  several  steps  have 
been  taken  to  make  the  data  gathering  as  simple  and 
automatic  as  possible. 

The  first  Triad  station,  established  in  1973  at  Fair- 
banks, Alaska,  was  manually  operated.  It  produced 
an  analog  tape  that  was  converted  to  digital  form  at 
APL  for  computer  processing  and  archiving.  The  sec- 
ond station,  set  up  in  1974  at  McMurdo,  Antarctica, 
included  automatic  equipment  to  record  analog  tapes 
that  likewise  were  converted  to  digital  form  at  APL. 
The  Fairbanks  station  was  backfitted  in  1975  with  the 
automatic  analog  recording  equipment  and  was  sub- 
sequently moved  to  Resolute  Bay,  Northwest  Terri- 
tories, Canada,  in  197 6.  In  1977,  a plan  was  conceived 
to  assemble  an  additional  station  in  Scandinavia  so 
that  Triad  data  might  be  correlated  with  data  collected 
by  the  European  GEOS  satellite.  This  presented  the 
opportunity  for  further  automation  of  the  collection 
process. 


DISCUSSION 

Two  costly  elements  in  the  production  of  Triad 
magnetometer  data  are  the  shipping  of  magnetic  tape 
and  the  conversion  of  data  from  analog  to  digital 
form.  The  third-generation  Triad  station  was  designed 
to  eliminate  this  conversion  by  producing  digital  tape 
directly  at  the  collection  station.  Digital  tape  has  a 
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packing  density  nearly  20  times  greater  than  that  of 
analog  tape;  thus,  shipping  costs  are  reduced.  In  order 
to  include  the  appropriate  time  and  control  data  used 
on  other  Triad  digital  data  tapes,  a considerable 
amount  of  digital  logic  is  required  at  the  collection 
station.  This  is  supplied  by  using  a microprocessor  as 
the  tape  controller.  The  presence  of  the  microproces- 
sor, in  turn,  allows  more  complex  acquisition  logic 
in  the  collection  station.  The  stations  now  at  Mc- 
Murdo  and  Resolute  accept  any  signal  on  Triad's  fre- 
quency, which  in  recent  times  includes  signals  from 
three  additional  satellites.  The  new  station,  with  its 
ability  to  identify  Triad  data,  can  eliminate  the  re- 
cording time  currently  spent  on  these  other  signal 
sources. 

The  third-generation  Triad  data-collection  station 
was  assembled  at  APL  using  a prototype  microproces- 
sor. Preliminary  tests  indicate  an  improvement  in  data 
quality  through  the  on-site  digitizing  process.  The  in- 
tegral microprocessor  opens  the  door  to  further  auto- 
mation of  the  highly  productive  Triad  data-collection 
effort. 

Figure  1 is  a block  diagram  of  the  Triad  3 station. 
The  omnidirectional  antenna  is  mounted  on  the  roof 
of  the  equipment  site,  and  all  other  equipment  is  con- 
tained within  the  building  in  two  4-ft-high  relay 
racks.  The  instruments  include  the  preamplifier,  re- 
ceiver, phase  tracking  filter,  and  synchronizer  and  are 
used  to  produce  digital  Triad  telemetry  data.  The  time 
code  generator,  beep  word  detector,  and  time  code 
reader  (in  the  tape  controller)  form  the  basic  timing 
system,  which  produces  epoch  time.  The  tape  con- 
troller (including  the  microprocessor)  and  the  tape 
recorder  form  the  data-recording  elements. 

The  processor  in  the  tape  controller  performs  three 
functions.  It  annotates  time;  inputs,  formats,  and 
qualifies  telemetry  data;  and  outputs  data  for  digital 
recording.  The  time  is  input  in  IRIG-B  serial  time 
code.  A decoder  packs  the  bits  into  words  that  are 
converted  to  binary  by  the  processor.  The  telemetry 
synchronization  pattern  is  used  to  establish  the  data 
frame  beginning  and  end.  If  the  data  are  of  good 
quality,  the  frame  is  released  for  recording  on  tape.  If 
the  data  are  from  another  source,  they  arc  rejected, 
ami  the  tracking  filter  is  returned  to  the  search  mode. 
The  recording  element  writes  data  records  consisting 
of  time  and  telemetry  data  on  the  digital  magnetic 
tape. 

Control  functions  performed  by  the  processor  make 
use  of  the  three  basic  elements  to  force  the  tracking 
filter  to  search  for  another  signal  and  to  close  data 
files  at  the  end  of  satellite  passes.  Each  time  a good 
data  frame  is  received,  a time-out  is  set  to  two  min- 
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Fig.  1 Block  diagram  of  the  Triad  3 station. 


utes.  When  this  delay  is  timed  out  at  the  end  of  a 
pass,  the  search  system  closes  the  pass  data  file  on  tape 
and  forces  a new  search.  If  a file  is  closed,  the  search 
is  set  to  delay  for  one  hour,  when  the  tracking  filter 
is  again  directed  to  search.  After  the  delay,  and  until 
more  good  data  are  received  and  recorded,  the  search 
command  is  given  every  two  minutes. 

Since  the  Triad  satellite  moves  at  a speed  of  more 
than  8 km/s,  an  accurate  time  base  is  required  to  cor- 
relate the  data  with  the  actual  satellite  position.  Local 
time  standards  are  available  in  the  McMurdo  and 
Resolute  stations,  but  none  is  available  at  the  potential 
sites  in  Scandinavia.  The  station  receiver  is  therefore 
equipped  with  an  alternate  crystal  to  allow  reception 
of  the  Navy  Navigation  Satellite  signals,  which  in- 
clude a precision  two-minute  time  mark.  The  time 
mark  is  used  to  check  and  correct  the  local  time  epoch. 
In  the  Triad  3 station,  this  will  be  done  manually  as 
part  of  a local  maintenance  check  of  the  time  code 
generator.  Considering  the  oscillator's  stability  of  I 
part  in  10',  the  drift  is  estimated  to  be  less  than  1 ms 
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per  day.  Weekly  testing  would  allow  limekeeping  on 
the  order  of  10  ms,  sufficient  to  calculate  the  satellite 
position  to  approximately  100  m. 
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COMPLETION  OF  THE  ORBIT 
DETERMINATION  PROGRAM 


An  orbit  determination  program  has  been  com- 
pleted tor  the  Navy  Navigation  Satellite  System.  This 
large  and  complicated  program  is  the  culmination  of 
seventeen  years  of  effort  in  space  technology.  Coded 
in  PL/I  computer  language,  it  is  designed  to  serve 
the  needs  of  the  Navy  Navigation  Satellite  System 
through  1990  and  beyond.  In  addition,  it  will  serve  as 
a valuable  research  tool  for  APL. 


BACKGROUND 

A major  component  of  the  Navy  Navigation  Satel- 
lite System  is  a computer  program  that  models  satel- 
lite orbit  trajectories.  The  precision  of  the  system 
depends  on  the  ability  of  this  program  to  determine 
and  predict  orbits. 

Prior  to  I960,  the  feasibility  of  doppler  tracking 
was  researched  and  demonstrated.  This  resulted  in  the 
current  orbit  program,  implemented  on  the  IBM  7094 
20  computer  at  the  Navy  Astronautics  Group  in  Cali- 


fornia since  1963.  That  program,  written  entirely  in 
assembly  language,  is  very  large  (more  than  400  000 
machine  instructions  and  more  than  20  core  overlays 
on  the  7094  computer).  Orbit  computations  for  the 
current  constellation  of  five  satellites  require  approxi- 
mately eight  hours  of  computer  time  per  day.  During 
its  14-year  life,  the  program  was  subjected  to  a num- 
ber of  revisions  and  improvements,  but  fundamental 
limitations  remain. 

In  1968,  the  design  of  a replacement  program  to 
run  on  a third -generation  computer  system,  such  as 
the  IBM  360  series,  was  started  However,  implemen- 
tation was  delayed  because  the  7094  continued  to  be 
reliable  and  because  it  was  necessary  to  obtain  a waiver 
from  the  Navy  to  use  the  pl/i  language  (in  1968 
pl/i  was  not  a "standard"  computer  language).  The 
scheduled  cessation  of  7094  maintenance  by  IBM  at 
the  end  of  1979  prompted  completion  of  the  program 
and  acquisition  of  an  IBM  360, '6S  computer,  which  is 
roughly  equivalent  to  the  7094. 


DISCUSSION 

The  new  Orbit  Determination  Program  (ODP)  is 
written  almost  entirely  in  the  PL/l  programming  lan- 
guage. This  higher-level  language  will  make  mainte- 
nance of  the  program  easier  because  it  is  more  read- 
able than  assembly  language.  Also,  since  PL/l  is  now 
a national  standard  language,  more  computers  can 
run  the  program.  ODP  is  still  a large  program,  occu- 


tinuing  improvement  in  the  ability  of  the  program  to 
determine  and  predict  orbits.  Figure  1 shows  the  im- 
provement in  prediction  precision  over  the  years.  The 
improvements  have  resulted  from  changes  in  geodesy, 
atmospheric  drag  models,  tropospheric  correction,  and 
numerical  techniques. 


pying  about  700  000  bytes  of  disk  space  and  running 
with  eight  overlays  in  a region  of  500  000  bytes  of 
core  storage.  Orbit  computations  for  the  constellation 
of  five  satellites  will  require  approximately  five  hours 
of  360/65  computer  time. 

Orbit  determination  and  prediction  involve  many 
disciplines  including  geodesy,  celestial  mechanics, 
mathematics,  and  computer  science.  The  requirement 
for  interaction  with  Navy  Navigation  Satellite  System 
hardware  results  in  a very  complicated  technical  prob- 
lem. As  a minimum  requirement,  ODP  must  predict 
orbit  trajectories  for  satellites  at  an  altitude  of  1000 
km  about  48  hours  into  the  future  with  a precision  of 
180  m.  The  prediction  is  based  on  an  orbit  determined 
from  36  hours  of  accumulated  doppler  tracking  data. 
The  program  consistently  does  better  than  this,  with 
prediction  errors  of  less  than  50  m.  The  predicted 
ephemeris  must  then  be  packed  and  formatted  accord- 
ing to  the  design  of  the  satellite  memory. 

The  largest  single  subsection  of  the  program  is  the 
computation  of  the  forces  that  act  on  the  satellite.  The 
gravitational  force  model  is  a series  expansion  in  har- 
monic functions  of  hundreds  of  terms.  Other  modeled 
forces  include  the  attraction  of  the  sun  and  moon,  the 
indirect  effect  (body  tides)  of  these  bodies  on  the 
shape  of  the  earth,  the  force  due  to  sunlight  (radia- 
tion pressure),  and  the  "drag'’  of  the  thin  upper 
atmosphere  on  the  satellite. 

The  mathematics  in  ODP  includes  the  numerical 
solution  of  a second-order  differential  equation  by 
finite  difference  techniques,  extensive  use  of  poly- 
nomial interpolation  schemes,  and  both  least-squares 
and  Chcbyshev  approximation  algorithms.  Interleaved 
with  them  arc  an  abundance  of  coordinate  transforma- 
tions, models  of  the  forces  on  the  satellite,  and  models 
of  the  signal  propagation.  Underlying  all  is  the  theory 
of  doppler  tracking,  the  basis  for  the  whole  system. 

Lise  of  the  orbit  program  as  a research  tool  has  pro- 
duced a kind  of  bootstrap  improvement.  F.xperimental 
versions  of  the  program  were  created  to  study  certain 
phenomena,  and  the  studies  produced  refined  models 
of  the  phenomena  that  were  then  included  in  the 
operational  program.  Other  refinements  resulted  from 
long-term  studies  of  program  performance  and  work  by 
other  investigators.  As  a result,  there  has  been  con- 


1963  1965  1966  1968  1974  1975 

Year 


Fig.  1 Chronology  of  prediction  improvement,  in- 
dicating years  when  program  changes  were  imple- 
mented. 


ODP  will  be  delivered  to  the  Navy  Astronautics 
Group  in  January  1978.  Comparison  runs  with  the 
7094  program  will  be  made  during  1978,  and  the 
program  is  scheduled  to  become  operational  in  1979. 
Preliminary7  comparisons  indicate  that  ODP  is  more 
precise  than  its  7094  counterpart.  It  is  also  more  effi- 
cient, requiring  roughly  half  the  time  of  the  7094 
program  (when  running  on  a comparable  computer), 
and  is  easier  to  maintain  since  it  is  coded  in  a higher- 
level  language. 
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NUMERICAL  METHOD  FOR  SOLVING 
ABLATION  PROBLEMS  AS  INVERSE 
HEAT  CONDUCTION  PROBLEMS 


It  is  cumbersome  to  use  conventional  mathemati- 
cal methods  to  account  for  ablation  (removal  of  ma- 
terial) in  the  temperature  analysis  of  a multidimen- 
sional reentry  body.  Because  of  the  need  for  accurate, 
economical  calculations  of  reentry  ablation  and  tem- 
peratures by  the  APL  Aerospace  Nuclear  Safety  Pro- 
gram. an  improved  numerical  method  has  been  de- 
vised that  eliminates  the  cumbersome  features  and 
yields  accurate  ablation  and  temperature  histories. 


BACKGROUND 

Ablation  resulting  from  melting  or  sublimation  of 
exposed  surfaces  is  an  inevitable  consequence  of  the 
severe  heating  that  occurs  when  a space  vehicle  re- 
enters the  earth’s  atmosphere.  The  ablative  process 
generally  has  a significant  effect  on  the  reentry  body's 
surface  and  interior  temperatures  because  energy  that 
would  otherwise  heat  the  body  is  used  for  phase- 
change. 

The  primary  task  of  the  Aerospace  Nuclear  Safety 
Program  is  to  perform  reentry-temperature  and  struc- 
tural analyses  of  nuclear  heat  source  assemblies  that 
power  radioisotope  thermoelectric  generators  for  such 
spacecraft  as  Pioneer  and  Voyager.  Its  object  is  to 
determine  if  the  heat  source  assembly  can  survive  re- 
entry in  the  event  of  an  abortive  launch.  For  reentry 
protection,  the  assemblies  have  a graphite  shell  that 
ablates  by  sublimation.  To  perform  the  analyses,  a 
numerical  technique  is  needed  that  calculates  the  re- 
quired temperature  histories,  taking  proper  account  of 
ablation. 

Since  the  heat-transfer  equations  and  the  system 
geometry  are  very  complex,  the  equations  are  best 
solved  by  finite-difference  techniques  that  satisfactorily 
approximate  the  procedure  in  which  a continuous 
system  is  modeled  with  a discrete  system  of  small  cells 
(see  Figs.  1 and  2).  The  cells  contain  the  points  for 
which  temperatures  are  to  be  calculated.  The  finite- 
difference  equations  are  approximations  of  the  differ- 
ential heat-transfer  equations.  Cell  definition  is  made 
to  guarantee  compatibility  between  the  difference  and 
differential  equations  so  that  as  the  cells  become 
smaller,  the  finite-difference  solution  approaches  the 
exact  solution. 

Conventional  finite-difference  algorithms  become 
cumbersome  when  ablation  is  introduced  because  the 
system  of  cells  that  represents  the  body  has  to  be  dis- 


torted continually  to  account  for  body  shape  changes 
(Ref.  1).  Figure  1 shows  a possible  two-dimensional 
model  distortion  from  an  originally  rectangular  cell 
structure.  This  distortion  is  required  in  order  to  main- 
tain compatible  cell  definition.  Three-dimensional 
ablation  problems  are  intractable  when  analyzed  using 
conventional  finite-difference  techniques.  The  elimina- 
tion of  the  cell  distortion  problem  in  the  numerical 
method  described  here  allows  two-  and  three-dimen- 
sional problems  to  be  treated  more  easily,  permits 
simpler  finite-difference  equations,  and  yields  faster 
and  more  accurate  temperature  and  ablation  predic- 
tions. Ablation  problems  that  could  not  have  been 
solved  previously  can  now  be  solved  easily  and  eco- 
nomically. 


Local  recession  distance 
Ablating  surface 
Physical  solid 


Fig.  1 Possible  finite-difference  cell  structure  for  an 
ablation  problem  posed  as  a conventional  heat-conduc- 
tion problem.  Only  physical  heating  is  applied  to  the 
exterior  surfaces  of  the  physical  solid,  which  the  finite- 
difference  cells  model  directly.  (TCP:  temperature  cal- 
culation point;  CB:  cell  boundary.) 


DISCUSSION 

An  undistorted  cell  geometry  can  be  retained  by 
embedding  the  ablation  problem  in  an  "inverse  heat- 
conduction  problem"  that  permits  no  material  re- 
moval. Figure  2 shows  how  the  cell  distortion  in  Fig. 
1 is  eliminated  by  the  introduction  of  fictitious  solid 
in  the  ablated  region  so  that  a rectangular  cell  struc- 
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Fig.  2 Finite-difference  cell  structure  for  an  ablation 

Eroblem  embedded  in  an  inverse  heat-conduction  prob- 
m.  Physical  heating  is  applied  to  the  ablating  and 
nonablating  surfaces  of  the  physical  solid.  Fictitious 
heating  is  applied  to  the  embedding  surface. 

turc  is  retained.  To  keep  the  heat  transfer  problem 
well  posed,  fictitious  heating  rates  are  applied  to  the 
embedding  surface  shown  in  Fig.  2.  The  rates  are 
determined  by  relating  them  to  the  known  ablating 
surface  temperature,  giving  the  problem  its  inverse 
character. 

In  contrast,  the  problem  in  Fig.  1 is  a conventional 
heat-conduction  problem  because  all  enclosing  bound- 
ary information  is  specified  directly.  Boley  (Ref.  2) 
first  proposed  embedding  ablation  problems  in  in- 
verse heat-conduction  problems;  the  extension  to 
finite-difference  analysis  was  first  proposed  by  the 
author  in  Refs.  1 and  3 and  further  elaborated  in  Ref. 
4.  A general-purpose  computer  program  was  written 
that  uses  embedding  (Ref.  5)  and  is  applicable  to 
one-,  two-,  and  three-dimensional  rectangular,  cylin- 
drical, and  spherical  geometries. 

The  embedding  computer  program  (designated 
SHTP-E)  borrows  heavily  from  APL's  existing  Stand- 
ard Heat  Transfer  Program  (SHTP,  Ref.  6).  It  uses 
accelerated  computational  methods  applicable  to  or- 
thogonal geometries  such  as  that  shown  in  Fig.  2 but 
not  to  distorted  ones  such  as  that  shown  in  Fig.  1 . The 
ablation  and  temperature  analysis  is  advanced  in  time 
steps  small  enough  to  guarantee  numerical  stability 
and  accuracy  but  large  enough  for  an  economical  cal- 
culation. Local  fictitious  heating  rates  are  related  to  the 
ablating  surface  temperature  and  net  heating;  this 
provides  a way  to  find  the  ablating  surface  tempera- 
ture. location,  and  rate  of  movement  as  well  as  the 
fictitious  heating  rate.  Once  the  fictitious  heating  rates 
are  found,  the  temperature  field  in  the  entire  body  of 
the  inverse  problem  is  determined,  and  the  calculation 
proceeds  to  the  next  time  level. 


References  1,  3,  and  4 compare  embedding  results 
with  those  from  more  conventional  methods.  Other 
comparisons  with  a more  exact  solution  can  be  found 
in  Refs.  1 and  3.  In  the  three  references,  the  methods 
agree  well  with  each  other,  and  the  superiority  of  the 
new  embedding  technique  for  the  analysis  of  ablation 
problems  is  demonstrated. 

Figure  3 is  a comparison  of  SHTP-E  predictions 
of  recession  distance  (defined  in  both  Figs.  1 and  2) 
for  a spherical  reentry  body  in  axisymmetric  flow  with 
those  made  by  CMA,  a widely  used  conventional  one- 
dimensional program  for  ablation  analysis  (Ref.  7). 
The  SHTP-F.  model  is  axisymmetric,  with  the  tem- 
perature depending  on  reentry  flight  time,  radius  from 
the  sphere  center,  and  angle  from  the  forward  stag- 
nation point  (the  forwardmost  point  on  the  sphere 
with  respect  to  the  airflow") . The  air  heating  of  the 
SHTP-E's  sphere  is  distributed  from  a peak  on  the 
forward  stagnation  point  to  a minimum  of  1 80° 
opposite  it.  The  CMA  model  ignores  angular  tempera- 
ture dependence  and  is  subject  only  to  stagnation 
point  heating.  For  both  models,  reentry  heating  begins 
11.06  s into  a reentry  flight  and  is  characterized  by  a 
brief  period  of  intense  heating  until  the  spheres  de- 
celerate and  the  heating  decreases.  Because  of  the 
angular  distribution  of  heat  conduction  and  external 
heating  in  the  SHTP-E  model,  these  results  differ 
from  the  CMA  results  in  an  expected  manner.  Further 
details  of  the  comparison  are  given  in  Ref.  4. 


Fig.  3 Predictions  of  recession  distances  for  a spherical 
reentry  body  in  an  axisymmetric  Aowfield.  Angles  are 
from  the  forward  point  of  the  reentry  body  facing  the 
flow.  The  SHTP-fc  finite-difference  model  is  axisym- 
metric and  is  subject  to  maximum  heating  at  0*,  de- 
creasing to  a minimum  at  180°.  The  CMA  model  ignores 
angular  heating  and  has  the  maximum  forward  point 
heating  applied  to  it. 
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SATELLITE  PATH  RAIN  ATTENUATION 
CORRELATED  WITH  PREDICTED  VALUES 
USING  RADAR  AND  DISDROMETERS 


The  use  of  radar  and  raindrop  size  distribution 
measurements  tor  determining  path  attenuation  sta- 
tistics has  been  demonstrated.  These  statistics  are 
essential  in  order  to  determine  margins  tor  power 
transmissions  and  receiver  sensitivities  in  designing  a 
satellite  configuration  with  a frequency,  f.  greater 
than  10  GHz. 


ttrements  to  obtain  case  by  case  fade  events  as  well 
as  fade  distributions  were  tested  successfully  against 
direct  measurements  with  the  COMSTAR  beacon  at 
28A6  GHz.  The  results  strongly  suggest  the  greater 
use  of  radar  for  the  acquisition  of  fade  statistics  for 
various  frequencies  and  path  configurations. 


BACKGROUND 

Because  of  the  increasing  demand  on  earth-satel- 
lite communication  systems,  frequency  hands  below 
10  GHz  are  besoming  crowded  and  therefore  prohib- 
itive. To  alleviate  this  problem,  future  systems  are 
planned  at  frequencies  above  10  GHz  where  the  atten- 
uation due  to  rain  may  represent  a severe  problem. 
To  design  highcr-frequcncy  systems,  designers  must 
know  the  attenuation  distribution,  that  is,  the  percent- 
age of  time  that  attenuations  (fades)  exceed  certain 
levels.  Such  information  is  used  in  establishing  trans- 
mitter power  margins  and  receiver  sensitivity  require- 
ments and  can  lie  acquired  directly  by  means  of  ex- 
perimental geosynchronous  satellites  with  ground 
receiving  or  transmitting  stations. 

An  indirect  way  to  determine  fade  distributions  is 
to  use  predictive  methods  on  radar  measurements, 
thus  minimizing  the  need  to  establish  expensive  and 
inefficient  satellite  communication  systems.  Predictive 
methods  using  radar  and  drop  size  distribution  mcas- 


DISCUSSION 

The  experimental  configuration  consists  of  a phase- 
locked-loop  receiving  system  operating  at  28. '6  GHz, 
an  S-hand  radar  (/  = 2.84  GHz)  located  30  m 
away,  and  a system  of  three  rain  gauges  and  two  dis- 
drometers  (for  measuring  raindrop  size  distributions) 
located  in  the  immediate  vicinity  of  the  receiving 
antenna.  The  receiving  system  has  a maximum  dy- 
namic range  of  33  dB  below  the  free  space  rw 
power  received,  which  is  — 106  dBm.  Both  the  re- 
ceiving and  radar  antennas  have  a bcamw  idth  of  0.4° 
and  look  in  the  direction  of  the  COMSTAR  satellite 
(elevation  41. 6°.  azimuth  210°). 

During  periods  of  rain,  the  radar  measures  the 
backscatter  of  rain  along  the  path  with  a range  reso- 
lution of  130  m.  A parameter  known  as  the  reflec- 
tivity factor,  7 (mm*/ms),  which  is  proportional  to 
the  backscatter  but  contains  only  the  characteristics  of 
the  rain,  is  measured  by  means  of  the  radar  (Ref.  1) 
at  each  range  bin. 
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The  disdrometer  is  an  electromechanical  sensor  that 
outputs  an  impulse  voltage  whenever  a drop  impacts 
on  a 50-cm*  cross-sectional  area  of  a Plexiglas  sensor 
(Ref.  2).  For  drops  falling  at  terminal  velocities,  the 
voltage  outputs  are  proportional  to  the  drop  momen- 
tums.  Through  calibration,  the  drop  diameters  may  be 
related  to  the  voltage  outputs,  and  the  ground  drop 
size  distributions  are  determined.  By  means  of  a given 
integration  time  interval  (e  g.,  30  s),  the  parameter  Z 
may  be  calculated  on  the  ground.  The  attenuation 
coefficient,  t (dB /km)  at  28.56  GHz,  an  additional 
parameter  calculated  using  the  measured  drop  size 
distribution,  measures  the  attenuation  per  unit  dis- 
tance of  propagation  due  to  the  sum  total  of  all  rain- 
drops. 

For  each  30-s  drop  size  distribution  measured 
with  the  disdrometer,  a pair  of  values  for  Z and  k 
may  be  calculated.  If  log  Z and  log  k are  plotted 
against  each  other,  a best-fit  linear  regression  line 
with  a relatively  high  coefficient  of  determination,  r2, 
may  be  established.  The  quantity  r 2 measures  the 
goodness''  of  fit,  1 being  perfect.  A straight  line  on 
a log  k versus  log  Z plot  implies  empirically  that 
& " uZ\  where  a and  b are  experimentally  deter- 
mined empirical  parameters. 

An  example  of  a regression  line  is  given  in  Fig.  1, 
which  corresponds  to  a 43-min  sampling  of  rain 
during  the  evening  of  14  September  1977  (20  h,  3 
min  through  20  h,  46  min  GMT).  The  vertical  and 
horizontal  axes  represent  log  k and  log  Z,  respec- 
tively. Each  of  the  43  data  points  represents  a 30-s 
integration  of  drop  size  distribution  for  which  the  rain 
rate  exceeded  2.5  mm  h.  The  straight  line  is  the  best 
fit  with  a - 3 47  X 10-’,  b - 0.733,  and  r • = 
0.95.  For  this  rain  period  as  for  others,  the  values  of 
a and  b were  inserted  into  a computer  program  that 
determines  the  radar-predicted  attenuation  as  a func- 
tion of  time.  For  the  events  considered  in  this  study, 
the  path  length  was  taken  only  up  to  the  0°  isotherm, 
above  which  ice  predominates. 

Figure  2 is  an  example  of  a corresponding  atten- 
uation event  within  the  interval  of  time  for  the  dis- 
drometer sampling  (20  h,  3 min  to  20  h,  31  min 
GMT).  The  vertical  scale  represents  the  path  atten- 
uation due  to  rain  and  the  horizontal  scale  the  time 
in  GMT.  The  solid  curve  represents  measured  atten- 
uation at  the  satellite  receiver  (sampled  every  3 s). 
The  dashed  curve  represents  the  predicted  attenuation 
using  the  radar  reflectivities  and  the  values  of  a and  b 
determined  from  the  disdrometer  measurements.  Also 
plotted  (the  dot-dashed  curve)  is  the  attenuation  pre- 
dicted using  the  radar  reflectivities  and  the  values  of 
a and  b derives!  from  the  well-known  Marshall- 
Palmer  (hereafter  referred  to  as  M-P)  distribution 


Fig.  1 Scatter  diagram  of  logio  of  the  attenuation  co- 
efficient plotted  versus  login  of  the  reflectivity  factor. 
Each  point  represents  a 30-s  sample  of  drop  size  distri- 
bution for  14  September  1977  from  20  h.  3 min  to  20 
h,  46  min  GMT  and  for  rain  rates  exceeding  2.5  mm/h. 


(Ref.  3),  namely,  a *=  2.01  X 10-’  and  b - 0.773. 
The  radar-predicted  levels  were  sampled  every  10  s. 

The  receiver  results  apparently  reached  saturation 
at  29  dB.  This  represents  a loss  of  lock  condition 
where  the  overall  receiver  dynamic  range  of  33  dB  is 
exceeded,  the  additional  4 dB  being  due  to  the  diurnal 
variation  of  the  satellite  as  well  as  to  water  on  the 
antenna  feedhorn.  The  loss  due  to  the  diurnal  varia- 
tion rarely  exceeded  2 dB  and  is  generally  known  to 
within  ±0.5  dB.  The  reflection  losses  due  to  water  on 
the  antenna  feedhorn  were  measured  by  controlled 
water  spray  tests  and  averaged  2 dB  with  an  uncer- 
tainty of  ±1  dB  when  a hydrophobic  agent  (Fusidox 
M)  was  used.  The  disdrometer-derived  results  and 
the  M-P  results  are  similar  and  show  generally  good 
correlation  below  29-dB  attenuation.  Above  this  level, 
one  to  one  correlation  is  not  possible. 

Since  the  designer  of  communications  systems  is 
interested  in  fades  from  a statistical  point  of  view, 
radar-predicted  fade  distributions  were  compared  with 
those  measured  directly.  During  the  summer  of  1977 
(June  through  September),  simultaneous  data  for 
five  rain  days  were  obtained.  Table  1 summarizes  the 
sampling  days  and  times  as  well  as  the  corresponding 
disdrometer  results.  When  all  the  drop  size  distribu- 
tions are  considered  together,  the  values  of  a and  b 
are  almost  identical  with  those  corresponding  to  the 
M-P  distribution.  This  is  attributed  to  the  fact  that 
the  M-P  distribution  represents  a large  data  base  for 
the  geographical  region  of  Washington,  D C.  (mid- 
Atlantic  east  coast).  On  an  individual  rain-period 
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Time  IGMT) 

Fig.  2 Comparison  of  attenuation  events  measured  directly  with  the 
COMSTAR  receiver  with  those  predicted  using  radar-disdrometer  measure- 
ments and  radar-M-P  distributions.  Event  of  14  September  1977. 


TABLE  1 

SUMMARY  OF  SAMPLING  DAYS  AND  TIMES  AND  CORRESPONDING  DISDROMETER  RESULTS 


Date 

(1977) 

Event 

No. 

Sampling  Time  Intervals 
(h:min) 

Disdrometer  Best-Fit  Parameters 

* 

b 

r2 

No.  of 
Samples 

Radar 

Disdrometer 

6 Ju 

1 

21  26-21:55 

21:26-22:45 

1.17X1 0* 

0.825 

0.93 

123 

2 

21:57-22:26 

* 

22:36-22:48 

9 Jun 

12:44-13:13 

12:44-14:49 

7.36X10* 

0.617 

0.75 

180 

5 

13:24-13:53 

6 

14:41-14:50 

24  Aug 

7 

14:56-15:24 

14:21-16:30 

1.23X10  = 

0.492 

0.60 

102 

8 

15:38-16:07 

9 

16:09-16:25 

24-25  Aug 

10 

23:35-23:59 

23:30-00:52 

1.45X10* 

0.806 

0.94 

107 

11 

00:05-00:34 

14  Sep 

12 

20:03-20:32 

20  03-20:46 

3.47X10* 

0.733 

0.95 

53 

1) 

20:33-20:48 

Notes:  Overall  radar  sampling  time  — 304  min  Best  fit,  b — 0.775 

Overall  disdrometer  sampling  time  - 458  min  (565  samples)  Best  fit.  r*  - 0.89 

Best  fit,  a — 1.87  X 10-*  (all  disdrometer  samples  for  all  days) 
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basis,  the  values  of  a and  b may  be  markedly  different 
from  the  M-P  distribution  (as  indicated  for  the  24 
August  period).  The  M-P  levels  of  attenuation  for 
this  period  disagreed  with  the  directly  measured  at- 
tenuation, whereas  the  disdrometer-predicted  values 
were  in  relatively  close  agreement. 

The  fade  distributions  from  the  approximately  five 
rain  days  of  sampling  (corresponding  to  13  attenua- 
tion events  and  5 h of  simultaneous  radar  data)  arc 
shown  in  Fig.  3.  The  vertical  scale  represents  the 
probability  that  the  attenuation  exceeds  the  abscissa, 
and  the  abscissa  corresponds  to  various  attenuation 
levels.  The  solid  and  open  circles  and  the  squares 
represent  the  directly  measured,  the  disdrometer-radar 
predicted,  and  the  radar-M-P  predicted  distributions, 
respectively.  On  an  absolute  basis,  the  disdrometer 
predicted  curve  is  in  good  agreement  with  the  directly 
measured  one,  whereas  the  M-P  curve  shows  sig- 
nificantly poorer  agreement. 

These  results  demonstrate  the  utility  of  using  radar 
in  conjunction  with  disdrometer  measurements  for 
predicting  fade  events  as  well  as  long-term  fade 
distributions  associated  with  satellite  communications 
through  rain. 
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Fade  depth  (d8) 


Fig.  3 Comparison  of  fade  distribution!  measured  di- 
rectly with  the  COMSTAR  receiver  with  those  pre- 
dicted using  radar-disdrometer  measurements  and  radar- 
M-P  distributions.  Curves  correspond  to  five  hours  of 
simultaneous  rain  data  covering  five  rain  days  and  13 
attenuation  events. 
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NASA  CONTROL  CENTER  DESIGN 


An  original  data  system  design  has  been  developed 
end  will  soon  be  implemented  by  NASA  at  the  God- 
dard Space  Flight  Center  (GSFC).  The  $5  million 
multiple-computer  data  system  will  satisfy  the  func- 
tional requirements  (Ref.  1)  of  a new  Network  Control 
Center  (NCC)  that  will  be  used  to  support  the  world- 
wide NASA  Spaceflight  Tracking  and  Data  Network 
during  1980-1990.  NCC  is  responsible  for  all  NASA 


spaceflight  tracking  and  data  acquisition  operations 
that  are  not  associated  with  the  Deep  Space  Network. 
It  must  ensure  that  data  from  scientific  satellites  and 
manned  spacecraft  are  constantly  supplied  to  the 
users.  It  must  also  provide  a communications  link 
between  a spacecraft  and  its  control  center  and/or 
experimenters,  who  may  be  located  at  a university  or 
other  non-NASA  location. 
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BACKGROUND 

NASA  is  planning  for  greater  use  of  automation 
in  network  management,  due  in  part  to  the  opera- 
tional Tracking  and  Data  Relay  Satellite  System 
(TDRSS).  Since  TDRSS  will  enable  a single  ground 
station  to  cover  most  satellites  and  manned  missions, 
several  overseas  ground  stations  may  be  phased  out. 
However,  the  data  system  still  must  support  approxi- 
mately eight  conventional  Ground  Spaceflight  Track- 
ing and  Data  Network  sites  along  with  the  one 
TDRSS  site,  and  typically  20  to  40  spacecraft.  In 
addition  to  the  automation  of  network  control,  the 
new  NCC  must  also  handle  much  higher  data  rates 
and  be  fault -tolerant. 

One  of  the  design  objectives  is  to  automate  as 
many  activities  as  practical  but  to  allow  for  human 
intervention  for  critical  decision-making.  To  satisfy 
this  objective,  a large,  centralized,  on-line  data  base 
is  required  for  mainraining  files  on  NCC  operations. 
A software  system  for  interactive  data  base  manage- 
ment is  also  necessary  to  allow'  network  operators  to 
monitor  and  control  the  operations. 


Input  data  rates  of  up  to  10  megabits  are  required, 
more  than  an  order  of  magnitude  greater  than  the 
present  NCC  system.  This  reflects  the  increased  data 
requirements  of  most  satellites  and  manned  space 
missions,  and  the  greater  number  of  spacecraft  that 
will  be  in  operation  in  the  future. 

NCC  must  provide  data- processing  services  24 
hours  per  day,  seven  days  per  week.  Operations  must 
be  near  real  time  in  order  to  respond  to  network 
anomalies  arising  from  link  failures,  loss  of  data,  or 
equipment  failures.  These  conditions  often  require 
the  network  operators  to  reconfigure  resources  inter- 
actively and  perform  emergency  scheduling.  At  times, 
test  and  simulation  are  required  before  this  can  be 
done.  Mandator)-  performance  requirements  were  im- 
posed on  the  data  system  design  to  ensure  sufficient 
fault-tolerance. 

Finally,  the  design  must  be  made  cost  effective  by 
minimizing  special  hardware/software  designs  and  by 
using  off-the-shelf  systems.  Also,  to  allow-  for  com- 
petitive bidding,  it  must  not  be  limited  to  a single 
system. 


Functional 

requirements 


Workload  requirements 


System  design 
considerations 


Configuration  studies  I 

— —I 


Selection  of  acceptable 
configuration 


CD  Inputs 
CD  Studies/decisions 
H Outputs,  results  of  studies 
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Fig.  t Design  methodology. 


DISCUSSION 

The  engineering  methodology  used  for  the  data 
system  design  is  illustrated  in  Fig.  1.  Various  inputs 
were  identified  and  trade  off  analyses  were  performed 
to  determine  which  functions  are  best  centralized  as 
opposed  to  those  best  distributed.  It  was  determined 


that  a configuration  for  both  centralized  and  dis- 
tributed functions  provided  the  best  solution.  The 
requirement  for  providing  a large  interactive  data 
base  that  is  accessible  and  functionally  assignable  to 
a large  number  of  terminals  was  best  satisfied  by  a 
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5 channels 
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processor 
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processor 


(message  handler 
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15  channels 


NASA  communications  network 


Fig.  2 Configuration  of  NCC  data  system. 
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centralized  system.  Trade-off  analyses  performed  on 
the  flow  of  communications  data  showed  that  these 
functions  must  be  distributed  to  accommodate  the 
number  of  linka  and  the  high  data  rates. 

Additional  trade-off  studies  were  made  on  the  use 
of  hardware  versus  software  while  minimizing  special 
designs.  A study  of  the  communication  links  revealed 
that  a considerable  amount  of  data  formatting/ 
reformatting,  packing/unpacking,  and  frame  synchro- 
nizing could  be  performed  more  efficiently  by  hard- 
ware rather  than  software.  This  resulted  in  a design 
that  required  special  communication  front-end  units 
that  could  handle  messages  or  monitor  telemetry  data. 

The  mandatory  performance  requirements  called 
for  a design  that  would  tolerate  single-point  failures 
without  losing  the  ability  to  support  time-critical 
control  functions.  The  design  must  also  provide  for 
both  scheduled  and  unscheduled  maintenance,  with- 
out disruption  of  services.  The  configuration  shown 
in  Fig.  2 provides  these  capabilities. 

The  communication  functions  are  accomplished  by 
means  of  four  identical  minicomputers  and  front -end 
units,  each  of  which  can  process  message  or  telemetry 
data.  A redundant  processor  is  required  for  backup 
and  software  development.  Analysis  of  the  network 
operations  showed  that  an  almost  equal  separation 
was  possible  between  control  and  support  functions, 
with  the  control  functions  being  time-critical.  This 
separation  allows  the  operational  functions  to  be  dis- 


tributed between  two  identical  operational  processors, 
thereby  avoiding  unnecessary  redundancy.  The  design 
also  provides  for  a "tightly-coupled"  multiprocessor 
system  with  shared  dual-port  memory  and  dual-access 
peripherals  for  auxiliary  storage.  Should  either  proces- 
sor fail,  the  high-priority  control  tasks  will  continue 
and  the  support  tasks  will  be  delayed  pending  re- 
pairs. 

Detailed  equipment  and  software  procurement 
specifications  (Ref.  2)  were  supplied  to  NASA.  They 
provided  for  extensive  potential  growth  and  the  flexi- 
bility for  the  further  distribution  of  functions  after 
the  system  becomes  operational.  Sperry  Univac  Com- 
puter Systems  was  selected  to  provide  the  total  system, 
which  will  consist  of  a Univac  1100/20  multiproces- 
sor, five  V77-613  minicomputers,  15  communication 
front  ends,  and  associated  intercomputer  interfaces 
and  system  software. 
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STIRLING  CYCLE  SATELLITE 
REFRIGERATOR  WITH 
ONE- YEAR  LIFETIME 


A Stirling  cycle  satellite  refrigerator  that  has  been 
designed,  constructed,  and  tested  will  cool  a germa- 
nium gamma-ray  detector  to  liquid  nitrogen  tempera- 
tures and  will  operate  in  the  satellite  environment  for 
at  least  one  year.  Four  of  these  refrigerators  are 
scheduled  to  be  launched  in  mid-1978  on  the  first 
satellite  experiment  designed  to  use  mechanical  re- 
frigerators to  obtain  high-resolution  gamma-ray  spectra 
In  space. 


BACKGROUND 

The  task  of  developing  a mechanical  refrigerator 
to  cool  a germanium  gamma-ray  detector  for  a satel- 
lite experiment  was  begun  by  APL  under  the  spon- 
sorship of  the  Defense  Advanced  Research  Projects 
Agency  (DARPA)  in  1973.  The  Stirling  cycle  re- 
frigerator was  selected  as  the  most  promising  candi- 
date, and  six  refrigerators  have  been  constructed. 


Previous  attempts  to  cool  radiation  detectors  in 
orbital  flight,  using  either  solid  cryogens  or  mechani- 
cal refrigerators,  are  reviewed  in  Ref.  1.  In  the  few 
instances  they  have  been  used  for  inf  rarest  detectors, 
mechanical  refrigerators  have  not  operated  more  than 
a few  hundred  hours  and  have  performed  erratically. 
The  ability  of  the  Stirling  cycle  refrigerator  to  cool 
detectors  to  cryogenic  temperatures  of  about  77  K 
with  relatively  good  efficiency  has  been  well  known 
for  some  years.  However,  it  has  been  difficult  to  ob- 
tain lifetimes  greater  than  several  hundred  hours  and 
to  cope  with  the  environmental  factors  associated 
with  satellite  launch  and  flight. 


DISCUSSION 


In  addition  to  providing  the  required  cooling  capac- 
ity, a satellite  refrigerator  must  meet  stringent  re- 
quirements for  weight,  power  consumption,  vibration 
transmitted  to  the  detector,  torques  transmitted  to 
the  satellite,  and  survival  during  launch.  The  APL/ 
DARPA  refrigerator  (Fig.  1)  has  met  the  following 
specifications: 

Output  refrigeration 
First  stage 
Second  stage 
Power  consumption 

Speeds  (three  ground - 
selectable) 

Weight 

Heat  exchange  temperature 
Vibration  of  cold  point 
Starting  torques 
Launch  vibration 


Lifetime 


140  K at  1.3  W 
63  to  73  K at  0.3  W 
30  to  33  W 

830,  1000,  1130  rpm 

13.8  lb 

0 to  43°C 
<0.00013  in. 

<0.01  in.-oz 

18.9  g rms  for  2 min 
each  axis 

7200  hours,  test 
continuing 


Five  of  the  refrigerators  have  been  delivered  to  the 
Lockheed  Palo  Alto  Research  Laboratory  (LPARL), 
which  will  fly  four  of  them  in  a satellite  experiment. 
The  remaining  refrigerator  is  being  life  tested  at 
APL.  The  cryogenic  portion  of  the  refrigerator  was 
manufactured  by  Philips  Laboratories  under  subcon- 
tract to  APL.  The  electronics  package  was  built  at 
APL;  it  contains  the  switching  circuitry  for  the 
brushless  nc  drive  motors,  the  speed-control  circuitry, 
and  the  signal-processing  circuitry  for  the  rpm,  gas 
pressure,  current,  and  temperature  measurements  of 


Fig.  1 The  APL/DARPA  refrigerator. 


the  motor,  which  are  telemetered  hack  to  earth.  All 
performance,  qualification,  acceptance,  and  life  test- 
ing has  been  performed  at  APL. 

The  mechanical  details  of  the  interior  of  the  re- 
frigerator and  the  solutions  to  the  construction  prob- 
lems are  described  in  Ref.  2. 

Unquestionably,  the  most  interesting  aspect  of  any 
satellite  refrigerator  is  its  performance  as  a function 
of  time  over  long  periods.  The  Serial  No.  1 refrigera- 
tor being  life  tested  at  APL  has  accumulated  an  un- 
precedented 7200  hours  of  operation.  Figure  2 shows 
the  first-  and  second-stage  temperatures  as  a function 
of  time  for  a continuous  run  of  360  hours.  The  tem- 
perature rise  is  about  1 K per  day  for  the  second 
stage.  When  the  refrigerator  was  stopped,  allowed  to 
warm  up,  and  then  restarted,  it  returned  to  (or  almost 
returned  to)  its  original  performance.  The  rate  of 
temperature  climb  shown  in  Fig.  2 is  related  to  the 
purity  of  the  helium  gas  within  the  refrigerator.  The 
rise  of  1 K per  day  or  less  can  only  be  achieved  with 
careful  outgassing  of  the  refrigerator  under  vacuum 
conditions  and  careful  charging  with  clean  helium. 

Superimposed  upon  the  slow  daily  temperature 
climb  is  a slower  climb  of  the  initial  start-up  tempera- 
ture (Fig.  3).  This  slow  rate  of  degradation  can  be 
caused  by  wear  of  the  displacer  seals,  gas  contami- 
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Heat  loads 

Second  stage  * 0.300  W 
First  stage  * 1 50  W 
Heat  sink  - 20'C 


Fin  2 Initial  performance  of  Serial  No.  I refrigerator 
as  a function  of  time. 


n.ints,  amt  loss  of  helium  pressure.  The  life  test  of 
the  Serial  No.  1 refrigerator  has  been  conducted  with 
no  maintenance,  disassembly,  or  gas  recharging,  and 
the  contributions  of  the  various  factors  to  degradation 
of  performance  cannot  be  assessed  until  the  test  is 
completed  The  results  shown  in  big.  .3  are  for  a re- 
frigerator speed  of  1000  rpm;  if  the  speed  is  increased 
to  ll'O  rpm,  the  second-stage  temperature  will  he 
about  75  K after  7200  hours  of  operation. 

The  intermittent  gamma- ray  experiment  to  which 
this  refrigerator  is  fitted  uses  an  intrinsic  germanium 
detector  It  does  not  require  continual  cooling  nor  is 
there  a requirement  for  continual  operation  for  pe- 
riods greater  than  200  hours.  The  resolution  and 
sensitivity  of  the  detector  do  not  degrade  at  tempera- 
tures less  than  120  K.  Thus  the  results  presented 
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Refrigerator  speed  * 1000  rpm 
Heat  loads: 

Second  stage  - 0 300  W 
First  stage  • 1.50  W 
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£ 140 
0> 

2 120 

1 100 

H 80 

60 

1000  3000  5000  7000 

Accumulated  time  (h) 


Fig.  3 Initial  temperatures  and  power  consumption 
(24  hours  after  start-up)  versus  accumulated  hours 
of  operation. 


show  that  the  API.  DARPA  refrigerator  meets  the 
requirements  of  the  experiment  for  which  it  was  de- 
signed. 
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BIOMEDICAL  SYSTEMS 


INTRODUCTION 


1 

The  Johns  Hopkins  University  School  of  Medicine,  the  School  of  Hygiene 
and  Public  Health,  and  API.  have  collaborated  in  a program  of  biomedical 
research  and  development  since  196V  An  important  objective  of  the  program 
is  to  apply  the  "know  how"  of  engineering  and  physical  sciences,  acquired 
by  APL  in  defense  and  civilian  research  and  development,  to  solve  problems 
in  medical  research  and  health  care  delivery.  The  program  has  included 
collaboration  with  nearly  all  of  the  clinical  departments  and  with  many  of 
the  basic  science  departments  of  the  medical  divisions.  More  recently,  a col- 
laborative program  has  been  established  with  the  Johns  Hopkins  Hospital. 
There  are  active  programs  in  molecular  biology,  nuclear  medicine,  ophthal- 
mology, neurosensory  research  and  instrumentation  development,  radiology, 
cardiovascular  systems,  prosthetic  systems,  computer  support  and  clinical  in- 
formation systems,  and  clinical  engineering. 

This  application  of  state-of-the-art  technology  has  contributed  to  many 
areas  of  basic  medical  research  and  to  clinical  diagnosis  and  therapy  through 
the  improvement  of  instrumentation,  techniques,  and  knowledge.  The  appli- 
cation of  systems  engineering  techniques  has  helped  solve  problems  of  the 
health  care  delivery  system.  These  activities  have  a focus  in  the  Department 
of  Biomedical  Engineering  at  the  Johns  Hopkins  Medical  Institutions.  In 
addition,  a new  educational  program  developed  in  collaboration  with  mem- 
bers of  the  Department  of  Biomedical  Engineering  leading  to  a Master's 
Degree  in  clinical  engineering  was  established  in  September  1973  at  the 
Medical  School. 

The  results  of  the  research  and  development  in  this  program  are  reported 
in  the  open  literature,  principally  in  biomedical,  biological,  and  medical 
journals.  During  the  program's  relatively  short  life,  more  than  1 so  papers 
and  several  book  chapters  have  been  published,  and  more  than  70  instru- 
ments for  research  and  clinical  application  have  been  developed.  A number 
of  API.  staff  members  have  assumed  line  responsibility  for  JHMI  in  areas 
where  technology  can  make  specific  contributions  to  patient  care  and  to  the 
health  care  delivery  system.  These  areas  include  radiation  physics,  information 


processing,  and  the  definition  of  technical  requirements  for  the  maintenance 
of  the  highly  technological  medical  equipment  used  for  patient  care. 

At  any  given  time,  there  are  40  to  60  active  projects  within  the  biomedical 
program.  They  range  all  the  way  from  fundamental  spectroscopic  studies  of 
model  biological  compounds — to  the  discovery  of  a method  for  observing  the 
blood  flow  in  usually  inaccessible  regions  of  the  human  eye— to  the  specifi- 
cation of  monitoring  equipment  and  communication  systems  for  better  patient 


I he  projects  reported  herein  include  efforts  in  fundamental  research  and 
in  the  biomedical  application  of  technology  developed  in  other  fields  of  APL 
activity.  The  first  article  demonstrates  the  utility  of  using  a mathematical 
model  to  describe  the  transport  of  ions  and  fluid  through  the  structured  bio- 
logical tissue  making  up  the  corneal  epithelium.  The  next  describes  measure- 
ments of  the  effect  of  intraocular  pressure  on  the  small  angle  scattering  of 
polarized  light  by  corneal  tissue.  The  third  describes  the  development  of 
exercise  echocardiography,  a new  technique  that  shows  promise  in  the  early 
detection  of  heart  disease.  The  fourth  deals  with  the  incorporation  of  a 
microprocessor  into  the  command  and  control  system  of  a medical  manipula- 
tor for  patients  with  high  spinal  cord  injury,  thus  vastly  reducing  the  number 
of  instructions  needed  to  control  the  manipulator. 


The  last  three  articles  describe  accomplishments  at  the  Johns  Hopkins 
Oncology  Center.  A computerized  information  system  gives  doctors  rapid 
access  to  medical  histories  for  treatment  or  research.  A program  for  testing 
and  calibrating  radiation  treatment  machines  has  been  developed  and  imple- 
mented. Finally,  a monitoring  system  using  data  telemetry  has  been  devel- 
oped for  use  with  ambulatory  patients. 


Of  course,  the  articles  that  follow  are  only  a sample  of  the  broad  range 
of  activities  within  the  collaborative  biomedical  engineering  program.  Other 
projects  will  be  highlighted  in  subsequent  Accomplishments  volumes. 
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COMPUTER  MODELING  OF  TRANSPORT 
IN  BIOLOGICAL  TISSUE 


A mathematics I model  of  transport  in  structured 
biological  tissue  has  been  constructed  and  used  to 
examine  the  transport  and  electrical  properties  of 
the  corneal  epithelium.  By  incorporating  a primitive 
structure  into  the  model,  it  becomes  possible  to  inter- 
pret experimental  data  on  entire  tissue  in  terms  of 
more  detailed  processes  ta king  place,  without  the 
errors  that  can  result  when  structure  is  ignored. 


BACKGROUND 

The  corneal  epithelium  is  a multiple  layer  of  cells 
at  the  very  front  of  the  eye,  bathed  by  the  tears  on 
the  surface  and  supported  by  the  corneal  stroma,  a 
layer  of  gel-like  connective  tissue  that  makes  up  most 
of  the  thickness  of  the  cornea  (Fig.  1 ).  Because  of  its 
gel  like  character,  the  stroma  has  a tendency  to  swell, 
and  stromal  swelling  is  known  to  decrease  corneal 
transparency.  Since  the  epithelium  is  the  only  barrier 
preventing  stromal  imbibition  of  the  tears,  it  is  im- 
portant to  understand  its  physiology  and  role  in  the 
maintenance  of  normal  corneal  thickness.  The  com- 
plex set  of  transport  pathways  through  the  epithelium 
cannot  lie  characterized  fully  by  available  experimental 
techniques;  the  purpose  of  this  study  was  to  gain  use- 
ful information  by  analyzing  experimental  data  in 
the  context  of  a simpler  structural  model  of  this  tissue. 


Fig.  I Three-dimensional  drawing  of  the  corneal  epi- 
(helium. 


DISCUSSION 

The  structure  of  the  epithelial  model  used  in  this 
work  (big.  2)  consists  of  a transcellular  pathway  in 
parallel  with  a shunt  between  the  cells.  This  structure 
is  much  simpler  than  that  of  the  real  epithelium  (com- 
pare Fig.  1 ) but  is  sufficient  to  explain  many  experi- 
mental data  iiii  this  tissue.  Also,  because  of  the  sim- 
plicity. the  parameters  of  the  model  (e.g.,  the  ionic 
permeabilities  of  the  shunt  and  cell  membranes  and 
the  rates  of  mctabolically  coupled  ion  transport  into 
and  out  of  the  cell)  arc  few  enough  to  permit  them  to 
be  evaluated  Irom  clcctrophysiologic  data  on  entire 
tissue  and  microelectrode  punctures  of  the  electrically 
coupled  epithelial  cells. 


Fig.  2 Structure  of  the  model  epithelium. 


In  the  model,  the  tears  and  stromal  fluid  can  con- 
tain any  number  of  permeant  monovalent  ions  (typi- 
cally Na\  K\  and  C.T)  and  charged  impermeants. 
Fluxes  across  the  shunt  and  the  cell  membranes  arc 
given  by  the  very  general  Ncrnst-Planck  equations 
that  relate  ion  flux  to  concentration  and  potential  dif- 


TABLE  1 

PREDICTED  AND  EXPERIMENTAL  THINNING  RATES 


Condition 

Thinning  Rate  (gm/h) 

Experimental. 
Refs.  2.  3 

Calculated 

R,  - 36  kSI-cm*  * 

R.  - 9 kil-cnC  * 

Resting 

- 1.0  ± 0.2  •* 

0 08 

0.2 

Stimulated  f 

1.3  ± 0.1 

0.7 

1.9 

Chloride-free 

bath,  resting  { 

3.9  * 0.6 

32 

5.6 

Chloride-free 

bath,  stimulated 

6.2  ± 0.6 

4.6 

7.6 

* These  bracket  the  experimental  variation  of  the  shunt  resistance,  R*  which 
is  variable  because  or  its  sensitivity  to  the  traumatic  effects  of  handling 


**  Minus  sign  denotes  swelling. 

t On  the  basis  of  other  data,  stimulation  is  modeled  as  a fourfold  increase  in 
the  permeability  of  the  outer  membrane  to  chloride  ion. 

* Chloride  in  tear-side  solution  replaced  by  impermeant  sulfate. 


ference  driving  forces.  Kinetic  expressions  for  the 
dependence  of  the  active  ion  transport  rates  on  sub- 
strate concentration  are  also  included.  The  model  is 
described  in  much  fuller  detail  in  Ref.  1. 

The  parameters  of  the  model  (e  g.,  the  permeability 
of  each  barrier  to  each  ion)  were  evaluated  largely 
from  intracellular  and  transepithclial  potential  meas- 
urements. The  parameterized  model  has  been  used  to 
predict  the  fluxes  across  the  epithelium  and  the  re- 
sponse of  these  fluxes  and  of  the  intracellular  and 
transepithelial  potentials  to  changes  in  ambient  com- 
position and  passive  permeabilities.  Of  more  general 
interest,  mode]  calculations  have  demonstrated  that 
several  experimental  techniques  used  to  measure  the 
transport  properties  of  biological  tissue  can  give 
misleading  or  erroneous  results  when  the  tissue  has  a 
structure  even  as  primitive  as  that  used  here. 

To  illustrate  the  predictive  capabilities  of  the 
model,  consider  an  experiment  (Refs.  2 and  3)  in 
which  a cornea  is  removed  from  the  eye  and  fluid 
exchange  at  the  back  surface  of  the  stroma  is  in- 
hibited by  a layer  of  oil.  Because  of  the  metabolic 
ion  pumps  in  the  epithelium,  there  will  be  a net  salt 
Hux  into  or  out  of  the  stroma.  This  will  induce  a 


corresponding  water  flow  by  osmosis,  and  the  stroma 
will  swell  or  thin.  The  rate  of  swelling  or  thinning 
can  be  influenced  by  altering  the  composition  of  the 
medium  bathing  the  front  surface  of  the  cornea 
(changing  the  driving  force  for  salt  flow)  or  by 
agents  such  as  adrenaline  or  theophylline,  which  alter 
epithelial  membrane  permeability.  Experimental  data 
on  the  stromal  thinning  rate  for  two  ambient  com- 
positions, at  rest  and  theophylline-stimulated,  are 
compared  with  predictions  in  Table  1 . Tile  agreement 
is  quite  good. 
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EFFECTS  OF  PRESSURE  ON  CORNEAL 
SMALL-ANGLE  LIGHT  SCATTERING 
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The  scattering  of  polarized  laser  light  by  the  cornea 
is  investigated  as  a function  of  transcorneal  pressure. 
Marked  changes  in  the  light  scattering  occur  as  the 
pressure  increases.  Also,  the  standard  procedures 
used  to  prepare  corneal  tissue  for  electron  micro s- 
copy  have  been  modified  so  that  a transcorneal 
pressure  is  maintained  during  the  fixation  step.  This 
modification  produces  structural  alterations  that  could 
explain  the  observed  changes  in  light  scattering  with 
applied  pressure. 

BACKGROUND 

The  cornea  is  the  clear,  tough  window  in  the  eye 
that  admits  light  to  the  photoreceptors  in  the  retina. 
The  cornea's  mechanical  strength  is  provided  hy  its 
constituent  rod  like  collagen  fibrils.  The  fibrils  scatter 
the  incoming  light,  and  modern  theories  recognize 
that  corneal  transparency  is  the  result  of  destructive 
interference  among  the  light  waves  scattered  in  all 
except  the  forward  direction  (Refs.  1 to  3).  These 
interference  effects  arise  because  the  spatul  arrange- 
ment of  the  fibrils  is  partly  ordered.  We  have  recog- 
nized that  light  scattering  provides  a noninvasivc 
pnibc  of  structure  anil  that  it  can  be  used  to  test  the 
predictions  of  various  structural  models  that  have 
been  proposed.  lrp  to  now.  by  comparing  light-scat- 
tering measurements  with  theoretical  calculations 
based  on  structures  shown  in  electron  micrographs 
(EM),  we  have  been  able  to  show  that  the  wave- 
length dependencies  of  the  observed  light  scattering 
arc  consistent  with  the  structures  depicted  in  F.M. 

Other  researchers  have  been  attempting  to  elucidate 
corneal  structure  by  using  the  small-angle  light- 
scattering (SALS)  methods  that  hive  been  a power- 
ful tool  in  understanding  polymer  morphology  (Refs. 
4 and  3).  SALS  differs  from  our  previous  methods 
in  that  polarization  effects  arc  investigated;  moreover, 
SALS  probes  larger  structures.  However,  comparisons 
of  the  observed  scattering  with  model  structures  for 
which  the  SALS  is  known  and  with  corneal  struc- 
tures known  histologically  have  not  been  satisfactory. 

Despite  these  difficulties,  our  success  with  conven- 
tional scattering  methods  led  us  to  believe  that  a 
similar  approach  should  work  with  corneal  SALS.  In 
this  regard,  we  noted  corrugations  in  the  lamellar 
structure  shown  in  EM  of  corneas  fixed  conventionally 
in  the  absence  of  pressure  (cf.  Fig.  la).  These  cor- 
rugations have  a fairly  regular  period  of  several  mi- 
crons and.  although  no  calculations  have  been  made, 
we  felt  that  the  waves  could  influence  corneal  SALS. 


(a) 


Fig.  1 Low-power  EM  of  corneal  stroma  of  pair  eyes 
from  same  rabbit,  (a)  Fixation  in  the  usual  manner  with 
no  applied  pressure.  Note  the  waviness  of  the  lamellae, 
(b)  Fixed  with  applied  pressure  of  18  mmHg.  Note 
that  the  lamellae  are  straight. 

Moreover,  J.  Francois  (private  communication)  had 
suggested  that  the  waviness  was  reduced  in  EM 
where  the  normal  intraocular  pressure  (IOP)  was 
maintained  during  fixation.  Previous  SALS  measure- 
ments were  made  without  maintaining  a transcorneal 
pressure;  consequently,  we  set  out  to  investigate  the 
effect  of  IOP  on  corneal  small-angle  light  scattering 
and  to  confirm  Francois'  suggestion  (Ref.  6). 

DISCUSSION 

In  the  SALS  technique,  the  cornea  is  illuminated 
with  a collimated  beam  of  polarized  light  from  a 
laser.  Light  scattered  with  polarization  either  parallel 
or  perpendicular  to  that  of  the  incident  beam  is  ob- 


served  at  scattering  angles  of  less  than  S°.  These 
scattered  light  intensities  form  patterns  that  are  called 
the  / .m.i  I patterns,  respectively.  The  /■  scattering 
pattern  results  primarily  from  isotropic  optical  density 
fluctuations  within  the  scattercr,  whereas  the  /,  scat- 
tering pattern  is  caused  by  anisotropic  optical  density 
fluctuations  The  / pattern  shows  the  most  significant 
effects  of  transcorneal  pressu  e 

Figure  2 shows  the  effect  of  increasing  IOP  on 
the  /.  scattering  from  the  central  region  of  a rabbit 
cornea  (Ref.  6).  The  effect  of  corneal  birefringence 
on  these  measurements  was  minimized  by  standard 
means  (Refs.  4 and  * ) . The  zero-pressure  pattern  is 
consistent  with  those  in  Ref.  5.  As  the  pressure  is  in- 
creased. the  zero-pressure  pattern  persists  but  its 
intensity  decreases  until  it  is  barely  discernible  at 
the  normal  IOP  (18  mmHg)  Above  zero  pressure, 
a new,  previously  unobserved  pattern  becomes  visible. 

We  suggested  (Ref.  6)  that  the  zero-pressure  pat- 
tern arises  from  the  waviness  in  the  corneal  lamellae 
that  is  seen  in  conventionally  fixed  EM  of  corneas 
(cf.  Fig.  la).  Subsequently,  we  devised  a relatively 
simple  means  for  EM  fixation  while  maintaining  a 
transcorneal  pressure.  When  fixed  at  18  mmHg,  the 
lamellae  are  almost  perfectly  straight  (cf.  Fig.  lb), 
which  confirms  Francois'  suggestion  and  provides 
strong  histological  support  to  our  suggestion  for  the 
origin  of  the  zero-pressure  scattering  pattern. 

Subsequent  to  our  original  experiments  reported  in 
Ref.  6,  we  have  been  attempting  to  explain  the  new- 
inner  pattern.  We  have  shown  that,  in  fact,  it  is  not 
a SALS  effect.  The  inner  pattern  is  actually  a projec- 
tion of  a cross  that  appears  over  the  whole  cornea 
when  it  is  illuminated  with  polarized  light  and  ob- 
served through  a crossed  polarizer  (Ref.  7).  In  our 
original  experiments,  we  did  not  eliminate  the  diffuse 
polarized  light  that  fell  over  the  whole  cornea  be- 
cause we  felt  that  it  was  insignificant  compared  to 
the  intensity  of  the  laser  beam  through  the  center. 
If  this  stray  polarized  light  is  carefully  removed  by- 
means  of  apertures,  both  the  corneal  cross  and  the 
inner  pattern  can  be  eliminated.  We  note,  however, 
that  the  origin  of  the  corneal  cross  has  never  been 
satisfactorily  explained,  and  we  arc  presently  attempt- 
ing to  learn  its  structural  basis. 
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QUANTIFYING  HEART-WALL  MOTION 
AND  THICKENING  BY  NONINVASIVE 
TWO-DIMENSIONAL  ECHOCARDIOGRAPHY 


Newly  developed  ultrasonic  imaging  instruments 
that  use  two-dimensional  electronic  scanning  have 
been  applied  to  the  visualization  of  individual  planes 
through  the  heart  and  to  the  quantification  of  the 
resultant  video  data.  These  data,  taken  while  the  sub- 
ject is  exercising,  can  differentiate  on  a regional  basis 
between  normal  motion  of  the  heart  wall  on  the  one 
hand  and  abnormal  motion  and  thickening  of  the 
heart  wall  on  the  other.  The  technique  shows  promise 
for  the  early  detection  of  heart  disease  and,  being 
noninvasive.  is  a candidate  tor  screening  people  in 
the  general  population  who  are  believed,  because  of 
other  factors,  to  risk  such  disease. 


BACKGROUND 

For  years,  the  two  tests  used  most  often  to  deter- 
mine the  presence  of  coronary  artery  disease  have 
been  cardiac  arterial  angiography  and  the  electrocardi- 
ogram (fkc)  taken  during  exercise.  Although  the 
FKG  method  is  noninvasive,  it  is  difficult  to  interpret, 
results  in  an  undesirable  number  of  false  positives, 
and  is  not  sufficiently  sensitive  to  mild  ischemia  (in- 
adequate oxygen  supply).  The  angiography  technique 
requires  that  a catheter  be  inserted  into  the  aorta  to 
deliver  opaque  contrast  material  to  the  coronary  ar- 
teries. which  are  then  viewed  by  X ray.  This  method 
is  highly  invasive  and  cannot  detect  mild  ischemia. 


The  echocardiographic  technique  is  totally  noninvasive 
and  shows  promise  of  detecting  changes  in  heart  wall 
motion  and  thickening  resulting  from  changes  in  the 
elastic  constants  of  the  heart  wall  caused  by  mild 
ischemia. 


DISCUSSION 

Figure  la  shows  one  video  frame  (taken  at  the 
rate  of  30  per  second)  of  a cross  section  of  the  left 
ventricle  of  the  heart  (white  represents  tissue,  black 
the  blood-filled  interior).  Several  frames  taken  at 
this  section  record  the  dynamics  of  the  ventricle  as 
the  heart  beats.  In  Fig.  lb,  computer-generated  con- 
tours are  superimposed  on  the  video  frame  of  Fig. 
la  and  are  modified  by  an  operator  who  moves  cir- 
cumferential points  (usually  16)  equally  spaced  in  an 
angle  about  the  center  of  the  area.  The  contours 
are  smoothed,  remodified  as  necessary,  and  stored. 
The  next  frame  is  processed  similarly  using  the  previ- 
ous contours  as  the  beginning  approximation. 

The  contours  arc  adjusted  by  the  operator  by  trans- 
lation, size  change,  and  radial  motion  of  a prese- 
lected number  of  points  After  contouring  ,V  frames 
covering  a full  cardiac  cycle,  the  N sets  ol  16  points 
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Fig.  1 (a)  Normal  two-dimensional  echo  cross  section  at  end  diastole  at 

the  level  of  the  papillary  muscles,  (h)  Same  section  with  inner  and  outer 
wall  contours  superimposed.  Papillary  muscles  are  excluded  from  this  con- 
tour. 


are  Fourier  transformed,  filtered,  and  inverse  trans- 
formed to  the  desired  number  of  time  samples.  Such 
data  from  planar  sections  can  be  combined  to  provide 
three-dimensional  volume  and  regional  wall  measure- 
ments. 

Figure  2 indicates  the  shortening  of  16  segments 
around  the  left  ventricle  at  the  end  of  ventricular 
contraction  (end  systole ) obtained  by  subtracting  the 
segmental  radius  at  the  beginning  of  contraction  (end 
diastole)  from  that  at  end  systole  Figure  2a,  for  a 
typical  heartbeat  of  a normal  control  subject,  shows 
relative  systolic  radial  shortening  around  the  whole 
left  ventricle  In  contrast.  Fig  2b  shows  the  same 
regional  motion  in  a patient  who  has  suffered  con- 
siderable damage  in  the  frontal  portion  of  his  left 
ventricle  from  a previous  heart  attack  (anterior  myo- 
cardial infarction).  One  can  see  the  segments  that 
exhibit  reasonable  motion,  those  that  show  relatively 
little  motion  (hypokinesis  or  akinesis),  and  those  that 
show  paradoxical  lengthening  (dyskinesis) . 

These  and  similar  results  demonstrate  the  ability 
of  the  system  to  obtain  quantitative  measurements  of 
regional  myocardial  performance.  The  major  efforts 
during  the  coming  year  will  be  to  validate  the  ultra- 
sonic results  by  comparing  them  with  results  from 
current  techniques  (eg.,  angiography  and  gamma 
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Fig.  2 Maximal  segmental  systolic  shortening  of  the  in- 
ternal left  ventricular  radius  (end-diastolic/end-systolic 
radius). 


camera  studies),  and  to  explore  techniques  for  auto- 
mating the  contouring  system. 
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MICROPROCESSOR  CONTROL  AND 
COMMAND  FOR  MEDICAL  MANIPULATOR 


The  application  of  a microprocessor  to  the  con- 
trol and  command  system  of  an  electrically  powered 
medical  manipulator,  formerly  implemented  with  dis- 
crete digital  and  analog  circuitry,  improved  the  manip- 
ulator s dynamic  performance  and  versatility.  Easily 
reprogrammable  sequences  command  the  manipula- 
tor through  a variety  of  motion  sequences  to  aid  the 
quadriplegic  patient  with  certain  everyday  tasks. 


BACKGROUND 

Quadriplegics  are  in  the  unfortunate  position  of 
having  little  or  no  control  of  their  trunk  and  limbs 
but  often  retain  the  ability  to  move  their  head  or 


their  facial  muscles.  These  muscles  may  be  used  to 
command  a powered  manipulator,  thereby  offering 
the  quadriplegic  patient  some  degree  of  independence 
in  performing  everyday  tasks  such  as  reading,  typing, 
eating,  and  telephoning. 

A six-degree-of-freedom  electrically  powered  ma- 
nipulator designed  by  APL  (Ref.  1)  gave  the  patient 
the  ability  to  execute  any  desired  sequence  of  motions 
with  the  aid  of  a chin  commander  (Fig.  1).  Physical 
control  consists  of  a linear  and  a discrete  command, 
the  linear  command  being  proportional  to  the  chin 
commander's  vertical  position  and  the  discrete  com- 
mand emanating  from  a microswitch  located  at  the 
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Manipulator 


far  cm!  of  the  chin  loop.  Motion  in  each  degree  of 
freedom  required  selection  and  termination  of  that 
degree  of  freedom  with  the  discrete  command  and 
positioning  with  the  linear  command  by  means  of 
servocontrol. 

Complex  sequences  of  motions  thus  required  many 
mentally  demanding  command  sequences.  The  need 
was  obvious  for  preprogrammed  command  sequences 
for  the  unvarying  motions  associated  with  many  tasks. 


DISCUSSION 

The  decision  to  use  a microprocessor  to  provide 
the  manipulator  with  the  capability  of  executing  pre- 
programmed command  sequences  also  provided  the 
opportunity  to  replace  much  of  the  discrete  logic  and 
analog  control  circuitry  with  microprocessor-executed 
programmed  logic. 

The  microprocessor  was  designed  to: 

1.  Scan  the  six  degrees  of  freedom  (or  program 
menu)  on  the  display  to  provide  the  patient 
an  opportunity  for  selection, 

2.  Monitor  the  position  states  of  the  manipulator 
for  display  and  control. 


3.  Update  the  display  when  necessary, 

4.  Monitor  the  linear  and  discrete  commands  for 
motion  control  and  program  sequencing, 

v.  Control  the  motion  in  each  degree  of  freedom 
so  that  it  is  smooth  and  safe,  and 

6.  Sequence  the  manipulator  through  the  com- 
mand sequences  of  the  selected  program. 

The  block  diagram  of  the  manipulator  (Fig.  2) 
shows  that  three  motors  driven  by  one  driver  move 
the  manipulator  in  six  selectable  degrees  of  free- 
dom. All  motions  except  track  displacement  are  exe- 
cuted with  cable-driven  spring-return  mechanisms. 
Motion  in  a specific  degree  of  freedom  begins  with 
the  establishment  of  a torque  balance  between  motor 
and  spring  so  that  a solenoid  lock  can  be  unlatched 
to  allow  motion  in  that  degree  of  freedom.  Linear 
commands  from  the  chin  commander  are  processed  to 
determine  the  direction  and  velocity  of  motion  sub- 
ject to  constraints  of  maximum  acceleration,  maxi- 
mum velocity,  and  position  limits.  Termination  of 
motion  involves  deceleration  to  a minimum  velocity 
to  allow  a smooth  latching  of  the  motion  at  the 
nearest  latch  point. 

Position  is  determined  with  an  8-bit  analog-to- 
digital  (A/D)  converter  and  high-quality  potentiom- 
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Fig.  2 Block  diagram  of  microprocessor-controlled  manipulator. 


cters.  The  digital  servomechanism  extracts  rate  infor- 
mation from  the  position  to  preclude  stalling  under 
varying  load  conditions.  The  display  and  digital  ma- 
nipulator commands  are  serially  updated  at  200  Hz 
when  the  servomechanism  is  being  used,  as  is  the 
linear  command  to  the  motor  driver  that  comes  from 
an  8-hit  digital-to-analog  (D/A)  converter. 

The  preprogrammed  command  sequences  are  exe- 
cuted by  an  interpreter  that  was  designed  to  imple- 
ment the  following  instructions: 

1.  Move  in  a selected  degree  of  freedom,  D,  to 
position  P. 

2.  Make  the  transition  from  programmed  mode 
to  manual  mode.  If  D is  selected  in  manual 


mode  and  the  motion  proceeds  past  position 
T,  return  to  the  programmed  mode  and  con- 
tinue the  motion  to  P. 

3.  Branch  to  any  interpreter  instruction. 

4.  Halt  program  execution  until  commanded  to 
continue  by  the  discrete  command. 

5.  Exit  program  to  return  to  manual  mode. 

The  F-8  microprocessor  used  to  implement  this 
particular  design  is  a highly  integrated,  low-power, 
8-bit  controller.  It  consists  of  a central  processing 
unit  (CPU),  a static  memory  interface  (SMI),  and  a 
parallel  input/output  chip  (PIO).  Two  interrupts, 
two  clocks,  32  hits  of  input  output,  addressability  to 
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65  556  bytes,  and  an  instruction  execution  time  of  2 
fis  are  provided. 

The  microcomputer  performs  its  control  task  in 
1500  bytes  of  memory.  An  additional  500  bytes  are 
used  for  interpreter  code  for  12  preprogrammed 
command  sequences. 
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ONCOLOGY  CLINICAL  INFORMATION 
SYSTEM  (PROVISIONAL  VERSION) 


The  use  of  complex  multimodality  treatment  plans 
and  the  necessity  for  the  detailed,  precisely  timed 
collection  of  data  quantifying  tumor  response,  drug 
toxicity,  biochemical  markers,  and  other  specific  pa- 
rameters of  therapeutic  evaluation  have  created  the 
need  tor  a clinically  oriented  data  system  than  can 
assist  the  oncologist  in  organizing , collecting,  and 
analyzing  data  related  to  research  and  treatment.  To 
serve  this  need,  a prototype  computer-based  system 
has  been  developed  and  is  in  routine  use  at  the 
Johns  Hopkins  Oncology  Center  (JHOC). 


BACKGROUND 

JHOC  is  one  of  18  Comprehensive  Cancer  Centers 
established  throughout  the  country  as  part  of  the 
National  Cancer  Plan  initiated  by  the  National  Can- 
cer Act  of  1971.  The  Center  has  major  programs  in 
laboratory  and  clinical  research,  education  at  all  levels 
of  training,  and  collaborative  activities  with  com- 
munity physicians.  As  part  of  its  clinical  activities,  it 
provides  cancer  patents  with  specialized  care  using 
advanced  techniques  in  chemotherapy,  radiation  ther- 
apy, surgery,  and  combined  modality  treatment. 

The  care  of  adult  patients  is  centered  in  a 56-bed 
facility  that  also  can  serve  500  outpatients  per  week. 
Inpatients  and  outpatients  are  usually  managed  under 
the  direction  of  similar,  well-defined  therapy  plans 
that  require  the  precisely  timed  administration  of 
drugs  and  the  timely  collection  of  clinical  and  lab- 
oratory data  quantifying  tumor  response,  drug  toxic- 
ity, biochemical  markers,  and  other  parameters  of 
therapeutic  evaluation.  Because  of  the  relative  com- 
plexity of  these  plans  and  the  need  to  maintain  rec- 
ords to  evaluate  an  individual's  response  compared  to 


that  of  other  patients  receiving  similar  treatment,  the 
Oncology  Center  has  developed  a prototype  informa- 
tion system  that  is  run  on  the  APL  computer. 


DISCUSSION 

The  Oncology  Clinical  Information  System  (OC1S) 
has  been  defined  as  a network  of  three  subsystems 
that  support  the  functions  of  center  management, 
patient  care,  and  radiation  therapy.  Because  funding 
for  a dedicated  computer  was  not  immediately  avail- 
able, it  was  decided  to  develop  a provisional  system 
using  the  APL  computer.  The  primary  objectives  of 
this  provisional  system  are  to  provide  useful  support 
in  patient  management,  to  give  designers  and  users 
an  opportunity  to  refine  system  requirements,  and  to 
develop  a data  base  for  later  versions  of  the  system. 

OCIS  operates  as  an  adjunct  of  the  Johns  Hopkins 
University  Medical  Information  System,  which  pro- 
vides on-line  support  in  the  areas  of  admissions,  out- 
patient registration,  inpatient  drug  distribution,  and 
retrieval  cf  selected  clinical  data  such  as  radiology 
reports.  The  system  also  performs  all  billing  func- 
tions and  produces  a variety  of  management  reports. 
Thus,  the  capabilities  of  the  prototype  system  are 
limited  to  supporting  functions  unique  to  the  On- 
cology Center.  The  organization  of  the  patient  data 
base,  the  type  of  data  displays  and  reports  available, 
and  potential  uses  of  the  system  are  discussed  briefly 
below. 

The  patient  data  base  is  organized  into  three  levels. 
The  highest  level  contains  summary  information  about 
all  patients  treated  for  cancer  at  Johns  Hopkins 


since  196-1.  Data  for  more  than  20  000  patients  are 
stored  in  records  that  contain  identification,  coded 
diagnoses,  dates,  and  follow-up  information.  Because 
the  data  base  may  he  listed  in  several  categories,  there 
is  a search  and  tabulation  capability.  The  primary  pur- 
poses of  this  file  are  to  serve  as  a master  index  to 
the  paper  records,  to  assist  in  follow-up  of  cancer 
patients,  and  to  identify  subsets  of  the  data  base  for 
further  analysis. 

The  second  data  level  contains  a complete  abstract 
tor  each  new  Johns  Hopkins  Hospital  cancer  patient. 
A data  base  of  4000  abstracts  has  been  accumulated 
since  the  initiation  of  OCIS.  Each  abstract  gives  all 
the  data  required  by  the  Hospital  Tumor  Registry 
and  serves  as  a nucleus  for  the  Oncology  Center 
patient  records.  The  file  combines  fixed  formatted  in- 
formation, coded  information,  and  free  text.  A por- 
tion of  a sample  abstract  is  shown  in  Fig.  1. 

The  third  level  contains  clinical  data  such  as  test 
results,  medications,  vital  signs,  and  blood  products. 
Such  data  are  collected  for  all  patients  treated  at 
[HOC!  and  arc  available  for  almost  1000  patients. 

The  main  purpose  in  collecting  the  clinical  param- 
eters is  to  produce  effective  data  displays  in  order  to 
understand  the  therapy  prixess  better.  The  provisional 
system  prepares  plots  and  flow  sheets  to  display  key 
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Fife.  1 Portion  of  sample  abstract. 

data  interactions.  Figure  2 is  a semilog  plot  of  plate- 
let count,  P,  and  white  cell  count,  II”,  versus  time. 
The  plot  also  shows  the  interrelationships  among 
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Fig.  2 Sample  plot  of  white  cell  count  and  platelet  data  versus  time. 
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platelet  anil  white  hint'd  cell  transfusions,  chemo- 
therapy, temperature  (decrees  above  37 °C),  and 
antibiotics  Alternate  plot  formats  create  histograms 
or  scatter  plots. 

In  addition  to  generating  plots,  the  system  also  pro- 
duces flow  sheets  that  contain  either  a general  sum- 
mary or  a detailed  survey  for  a specific  organ  system. 
Depending  on  their  use,  flow  sheets  can  be  printed 
horizontally  (Fig.  3)  or  vertically. 
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Fig.  3 Sample  flow  sheet. 
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Fig.  4 Sample  bacteriology  report. 


agement,  platelet  donor  matching,  and  other  func- 
tions. 


Finally,  a special  set  of  reports  is  prepared  from 
the  microbiology  data.  Because  time  is  required  to 
grow  a culture  and  identify  an  organism,  results  fre 
qucntly  are  delivered  several  weeks  after  the  speci- 
men has  been  collected.  To  make  effective  use  of  this 
information,  the  reports  are  sorted  by  date,  by  source 
of  specimen  (e.g.,  blood),  and  by  organism.  The 
sorted  reports  allow  the  physician  to  trace  an  organism 
from  site  to  site  and  to  note  changes  in  antibiotic 
sensitivity.  A sample  report  is  shown  in  Fig.  4. 

The  provisional  version  of  the  system  maintains 
an  appointment  system  that  schedules  outpatient  visits 
and  prepares  data  summaries  automatically  prior  to 
patient  visits.  There  is  also  a census  system  to  main- 
tain records  of  all  patients  treated  at  JHOC,  and 
there  arc  special  programs  to  support  facility  man- 


CONCLUSION 

The  provisional  system  was  put  into  operation 
with  the  opening  of  the  new  JHOC  facilities  in  Octo- 
ber 1976.  In  its  first  year,  it  satisfied  all  of  the  princi- 
pal objectives.  It  is  being  replaced  by  a second  ver- 
sion that  uses  a dedicated  minicomputer  located  at 
JHOC.  The  medical  staff  is  familiar  with  the  system 
and  its  outputs;  in  many  cases,  disease-specific  out- 
puts have  been  designed  by  the  JHOC  physicians. 
There  is  a large  and  accurate  data  base  to  support 
patient  management  and  research.  Moreover,  the  re- 
quirements for  the  next-generation  OCIS  are  well 
understood  and  should  be  relatively  easy  to  imple- 
ment. 

Author:  B.  I.  Blum 

Support:  Johns  Hopkins  Oncology  Center 
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CALIBRATION  OF  RADIATION  ONCOLOGY 
TREATMENT  MACHINES 


Emphasis  on  excellence  of  treatment  in  radiation 
oncology  requires  a major  physics  and  engineering 
effort  in  operating  and  calibrating  large  linear  accel- 
erators used  for  therapy.  During  the  preceding  year, 
all  major  therapy  machines  in  the  Johns  Hopkins 
Oncology  Center  were  calibrated  and  are  now  in  use. 
In  addition,  the  elements  of  a quality-control  program 
f or  accurate  dose  delivery  were  formulated  and  ini- 
tiated. Future  efforts  will  entail  consolidation  of  the 
various  activities  into  a unified  program  that  includes 
machine  maintenance  and  calibration,  treatment  plan- 
ning, clinical  dosimetry,  and  periodic  dose  verification 
to  assure  the  high  degree  of  excellence  demanded  in 
a modern  radiation  oncology  facility. 


BACKGROUND 

The  new  radiation  oncology  wing  of  the  Johns 
Hopkins  Oncology  Center  opened  in  1977.  The  fa- 
cility has  three  major  high-energy  therapy  machines. 
A Siemens  “"Co  teletherapy  unit  with  a nominal  pho- 
ton energy  of  1.25  MeV  is  the  facility's  only  isotopic 
gamma-ray  source.  In  addition,  two  linear  accelerators 
have  been  installed:  a Varian  Clinac-4  linear  accel- 
erator that  produces  X rays  at  4 MeV  and  a Clinac-18 
linear  accelerator  that  produces  X rays  at  10  MeV 
and  electrons  at  6,  9,  12,  15,  and  18  MeV.  The  three 
machines  and  a diagnostic  simulator  are  now  opera- 
tional. More  than  100  patients  per  day  are  being 
treated  in  the  Center. 


DISCUSSION 

Before  the  machines  were  used,  it  was  necessary 
to  measure  the  dosage  supplied  by  the  beams.  Ini- 
tially, the  dose  rate  of  each  machine  was  calibrated 
absolutely  by  measuring  the  beam  ionization  produced 
in  a water  phantom  (equivalent  to  tissue)  under 
standardized  conditions.  Baldwin-Farmer  ionization 
probes  and  electrometers  calibrated  at  the  National 
Bureau  of  Standards  were  used  to  convert  the  meas- 
ured ionization  to  absorbed  dose. 

In  addition  to  the  absolute  calibration,  relative 
dose  measurements  were  made  on  each  machine  to 
define  the  treatment  beam  parameters.  These  data 
were  obtained  by  driving  a scanning  probe  (under 
computer  control)  through  a water  phantom  that  was 
being  irradiated.  This  established  the  relative  dose 
versus  depth  and  also  the  relative  dose  over  the  lat- 
eral beam  geometry.  The  data  have  been  computerized 


and  are  used  for  the  detailed  planning  of  each  pa- 
tient’s treatment  by  means  of  computer  modeling. 
This  work  required  approximately  eight  months  of 
intensive  effort. 

A major  program  has  been  initiated  to  ensure  dose 
accuracy  and  the  repeatability  of  treatment.  The  major 
elements  of  the  program  are: 

1.  An  aggressive  preventive  maintenance  pro- 
gram consisting  of  daily  checks  of  specific 
parameters  and  quarterly  checks  by  the  vendor 
of  the  therapy  machine; 

2.  Daily  checks  for  constancy  and  symmetry  of 
the  output  level  of  the  beam,  by  means  of  di- 
odes, plus  weekly  checks  by  means  of  ioniza- 
tion chamber  measurements; 

3.  A comprehensive  weekly  review  of  all  patient 
charts  for  dose  verification;  and 

4.  An  extensive  program  in  clinical  dosimetry’ 
(i.e.,  measurement  on  the  patient)  using  di- 
odes and  thermoluminescent  dosimeters  (LiF). 

Many  of  these  elements  are  being  incorporated 
into  a computerized  dose-quality-control  program 
that  will  be  the  basis  of  a model  system.  The  new 
system  will  include  computerized  chart  checks  and  an 
automatic  dose-measuring  system. 

In  addition,  all  treatment  planning  and  dose  cal- 
culation algorthims  are  being  reviewed  in  a continu- 
ing effort  to  maintain  a library  of  only  the  most  up- 
to-date  computational  techniques. 

An  investigation  of  the  dosimetry  of  total-body- 
irradiation  (TBI)  patients  is  under  way.  Isodose  con- 
tours of  dose  levels  through  various  body  sections 
have  been  obtained  for  adult-  and  child-size  phan- 
toms. The  responses  of  hundreds  of  dosimeters  placed 
in  the  life-size  phantoms  are  processed  and  the  data 
(plane  by  plane)  are  interpolated  to  establish  the 
isodose  pattern. 

As  a part  of  the  overall  program,  a physics  course 
appropriate  for  technologists  and  residents  has  been 
initiated.  It  is  presented  as  part  of  the  accredited 
Technologist  School. 


Authors:  D.  G.  Grant  (APL)  and  W.  C.  Lam  (JHOC) 
Support:  Johns  Hopkins  Hospital 
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THE  PATIENT  MONITORING  SYSTEM 
AT  THE  JOHNS  HOPKINS 
ONCOLOGY  CENTER 


A complete  physiological  monitoring  system  has 
been  installed  and  is  operating  at  the  Johns  Hopkins 
Oncology  Center.  It  provides  information  on  basic 
vital  signs  at  a bedside  panel  and  also  in  central 
nurses'  stations  for  all  rooms  in  the  56-bed  inpatient 
facility  of  the  Johns  Hopkins  Oncology  Center.  The 
system  meets  physician/patient  requirements  and  has 
improved  patient  care  significantly. 


BACKGROUND 

The  physiologic  monitoring  of  patients  covers  an 
extreme  range  of  physician/patient  requirements.  On 
the  one  hand,  physicians  require  fundamental  data 
such  as  electrocardiogram  (ecg),  respiration  rate, 
heart  rate,  temperature,  and  blood  pressure  in  a for- 
mat that  provides  predictive  trends  and  also  has  an 
alarm  mode.  On  the  other  hand,  the  patient  should 
he  encumbered  as  little  as  possible  in  the  gathering 
of  such  data  and,  in  particular,  should  remain  ambu- 
latory if  possible. 


DISCUSSION 

The  Johns  Hopkins  Oncology  Center  is  a 56-bed 
facility  situated  on  two  floors  and  serviced  by  four 
nurses'  stations.  Each  nurses'  station,  the  focal  point 
of  the  patient  monitoring  system,  serves  14  to  16 
rooms.  All  rooms  function  identically  (Fig.  1).  Pa- 
tient monitoring  panels  behind  the  bed  relay  some  sig- 
nals (temperature  and  infusion  pump  alarm)  directly 
to  the  nurses'  station  and  relay  others  (primarily  ECG 
and  blood  pressure)  to  a second  wall  panel  located 
in  a position  somewhat  removed  from  the  bed.  This 
second  relay  panel  is  connected  to  a bedside  cart 
that  contains  a modular  amplifier  for  physiologic 
signals.  The  amplified  signals  are  transmitted  through 
the  same  panel  to  the  central  nurses'  station.  Each 
room  has  a pair  of  36-conductor  cables  leading  back 
to  the  nurses'  station. 

The  Oncology  Center  monitoring  system  ha  pro- 
vided for  16  hardwired  ecg  channels,  9 telemetry 
FCG  channels,  15  respiration-rate  channels,  16  heart- 
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Fig.  1 Schematic  of  patient  monitoring  system. 


rate  channels,  9 pressure-monitoring  channels,  2 
sphygomomanometer  channels,  and  1 oximeter  chan- 
nel (two  additional  channels  have  been  ordered), 
all  of  which  can  be  distributed  flexibly  among  all  the 
patient  rooms.  Also,  up  to  three  simultaneous  tem- 
perature channels  are  available  from  each  room. 

Installation  of  the  patient  monitoring  system  is 
95%  complete.  Two  cart  formats,  acute  and  intensive, 
are  available,  depending  upon  the  severity  of  the 
illness.  Signals  for  ecg  (both  hardwire  and  telem- 
etry), heart  rate,  respiration  rate,  arterial  and  venous 


blood  pressure,  and  patient  temperatures  are  avail- 
able. All  vital  physiologic  signals  have  variable  up- 
per and  lower  alarm  thresholds.  Visual  scope  dis- 
plays and  chart  recorders  facilitate  data  interpreta- 
tion. Future  plans  include  more  extensive  applica- 
tions of  data  telemetry. 

Authors:  D.  G.  Grant  (APL)  and  E.  Hertz  and  D.  M. 
Riker  ( JHOC ) 
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URBAN  TECHNOLOGY 


INTRODUCTION 


The  recognition  that  the  great  cities  continue  to  play  a key  role  in  main- 
taining the  vitality  of  the  nation  has  spurred  renewed  interest  in  the  institu- 
tional structures  that  support  urban  life.  Through  a variety  of  special  projects, 
APL  has  participated  in  and  contributed  to  programs  directed  at  developing 
and  applying  modern  technology  to  a variety  of  civil  problems  directly  rele- 
vant to  current  urban  issues.  These  programs  have  included  transportation, 
fire  research,  siting  of  pow’er  plant  facilities  near  urban  centers,  health  care 
delivery  systems,  and,  more  recently,  location  of  leaks  in  buried  natural-gas 
distribution  lines.  While  not  all  of  these  program  areas  are  included  in  this 
section,  several  are  brought  together  here  in  an  effort  to  provide  a unifying 
perspective  on  their  scope  and  their  relationship  to  urban  life. 

The  articles  in  this  section  relate  specifically  to  public  transportation  and 
to  questions  of  public  safety.  The  article  on  short-headway  (hence,  poten- 
tially high  volume)  automated  guideway  transit  systems  represents  APL’s 
continuing  interest  in  improving  urban  transportation  via  systems  analysis 
and  technological  innovation. 

The  article  on  leak  detection  in  buried  natural-gas  distribution  lines  con- 
cerns safety  and  resource  conservation.  The  explosion  hazard  from  leaking 
buried  gas  distribution  lines  has  been  well  documented,  and  several  unsuc- 
cessful engineering  attempts  have  been  made  to  develop  methods  of  pin- 
pointing leaks. 

However,  the  need  for  an  effective  method  continues  to  increase  since 
throughout  the  United  States  the  gas  distribution  systems  have  reached 
the  age  where  sharp  increases  are  projected  in  the  mean  annual  leak  rate  per 
mile.  The  APL  program  is  directed  at  analyzing  the  active  acoustic  method 
as  the  only  viable  current  method  of  leak  detection.  The  project  is  both 
theoretical  and  experimental  but,  unlike  prior  attempts,  has  a research  orien- 
tation and  includes  studies  of  acoustic  propagation  in  soils,  instrumentation 
development,  and  theoretical  analysis. 

The  last  two  articles  describe  different  aspects  of  the  APL  Fire  Program. 
The  first  concerns  the  combustion  characteristics  of  materials  commonly  in- 
volved in  fire,  including  the  quantitative  determination  of  surface  ignition 
temperatures  and  combustion  products.  The  second  concerns  the  clinical  and 
technical  evaluation  of  fire-related  deaths.  This  work  represents  the  broader 
Laboratory  involvement  in  fire  technology  and  fire-fighting  tactics.  The  devel- 
opment of  new  materials  for  construction,  including  some  that  produce  toxic 
agents  on  burning,  has  increased  the  threat  of  fire  to  both  occupants  and  fire 
fighters.  The  mortality  aspect  of  the  APL  program  attempts  for  the  first  time 
to  quantify  specific  contributors  to  deaths.  A third  component  of  this  program 
concerning  documentation  and  dissemination  of  information  related  to  fire 
appeared  in  previous  Accomplishments. 


PRECEDING  FAGS 


DATA  RATE  AND  QUANTIZATION 
REQUIREMENTS  FOR  A SHORT-HEADWAY 
VEHICLE-FOLLOWER  AGT  SYSTEM 
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An  important  consideration  in  designing  Automated 
Guideway  Transit  (AGT)  Systems  is  the  communica- 
tion requirements  tor  longitudinally  controlling  ve- 
hicles. Thus,  trade  off  studies  are  required  to  establish 
the  communication  data  rates  and  accuracies  needed 
to  maintain  an  acceptable  level  of  vehicle  perfor- 
mance. 


BACKGROUND 

Vehicle  following  in  an  AGT  system  is  a longi- 
tudinal control  scheme  in  which  the  state  of  a vehicle 
is  determined  by  the  state  of  the  immediately  pre- 
ceding vehicle.  It  has  been  shown  (Ref.  1)  that  at 
headways  (the  time  between  successive  vehicles  pass- 
ing a fixed  point  on  the  guideway)  of  less  than  3 s, 
a fundamental  kinematic  constraint  creates  bounds  on 
vehicle  motion.  This  constraint  dictates  a minimum 
allowable  spacing  between  vehicles  that  is  a function 
of  trailing  vehicle  state,  preceding  vehicle  state,  and 
the  future  maneuvering  capability  of  each  vehicle. 

The  difficulties  encountered  in  controlling  a ve- 
hicle subject  to  the  kinematic  constraint  have  been 
resolved  by  the  use  of  a controller  that  is  a nonlinear 
function  of  states  (Ref.  2).  The  approach  presented 
in  Ref.  2 explicitly  incorporates  the  kinematic  con- 
straint into  the  control  law  and  consequently  pro- 
vides near-optimal  performance  while  satisfying 
safety  requirements  under  all  operating  conditions. 
The  control  law  thus  provides  a baseline  system  for 
studying  the  effects  upon  vehicle  performance  of 
sampling,  quantization,  and  time  delays  associated 
with  intervehicle  communication  and  control  compu- 
tation. 

To  implement  a longitudinal  control  law,  there  are 
many  ways  to  configure  the  overall  system,  that  is,  to 
measure  vehicle  states,  to  transmit  this  information 
between  vehicles,  and  subsequently  to  control  the  in- 
dividual vehicles.  As  a result,  the  problem  of  the 
required  data  rate  creates  various  trade-offs  related  to 
vehicle/wayside  control  allocation  and  on-board  digi- 
tal/analog (D/A)  interfacing.  Control  allocation  de- 
termines whether  control  loops  should  be  closed  on 
board  the  vehicle  (a  "smart"  vehicle)  or  through 
the  wayside  (a  "dumb"  vehicle).  The  specification 
of  the  many  parameters  involved  in  the  problem  is 
approached  through  analytically  derived  results  using 
simplified  linear  models  and  then  verification  and  ex- 


pansion of  these  results  with  a detailed  computer 
simulation. 


DISCUSSION 

The  general  system  configuration  that  will  be  con- 
sidered for  this  study  is  shown  in  Fig.  1.  An  on- 
board digital  computer  receives  information  from  the 
wayside  that  is  sampled,  quantized,  and  delayed.  The 
on-board  digital/analog  interface  determines  how 
much  of  the  on-board  control  computation  is  per- 
formed digitally.  The  figure  reflects  a smart  vehicle 
configuration  where  all  control  computation  is  per- 
formed on  board  the  vehicle  and  where  preceding 
vehicle  information  is  received  from  the  wayside  in 
the  form  of  spacing  and  its  derivatives  (relative  ve- 
locity and  acceleration).  A dumb  vehicle  configura- 
tion would  allocate  control  computation  to  a wayside 
computer.  The  parameters  considered  in  this  study 
that  affect  sampling  rate  and  quantization  require- 
ments arc: 

1.  Headway  (0.3  and  3 s); 

2.  Smart  or  dumb  vehicle; 

3.  On-board  sampling  interval; 

4.  Wayside  sampling  interval; 

3.  Bit  number  to  represent  wayside  information; 

6.  Bit  number  for  on-board  computation; 

7.  Location  of  D /A  conversion  on  board; 

8.  Bit  number  for  D/A  conversion;  and 

9.  Inclusion  of  analog  filtering  following  D/A 
conversion. 

The  first  phase  of  the  study  was  to  consider  closed- 
loop  stability  as  a function  of  headway,  wayside  sam- 
pling rate,  and  control  allocation.  A z-transform 
analysis  in  conjunction  with  root  locus  plots  is  used 
for  this  purpose.  From  a stability  viewpoint,  the 
effects  of  sampling  arc  more  strongly  related  to  the 
jerk  and  acceleration  limits  on  the  vehicle  rather  than 
to  the  particular  operating  headway.  This  is  because 
the  control  law  derivation  is  based  on  the  kinematic 
constraint  and  thus  controller  bandwidth  is  largely 
determined  by  the  gains  that  reflect  this  constraint, 
independent  of  headway. 
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Fig.  1 General  system  configuration. 


Another  result  of  the  analysis  is  that  a smart  ve- 
hicle offers  an  advantage  over  a dumb  vehicle  in 
that  the  required  sampling  rate  from  the  wayside 
is  lowered.  For  a service  jerk  limit  (nonemergency) 
of  2 m/sa  and  a service  acceleration  limit  of  1.5 
m/s*,  a dumb  configuration  requires  a sampling  in- 
terval of  0.35  s to  maintain  stability  at  a 0.5-s  head- 
way (Ref.  3 V On  the  other  hand,  for  a smart  vehicle 
this  requirement  is  relaxed  to  a sampling  interval  of 
0.6  s.  At  a 3-s  headway,  the  requirements  are  1 s for 
a dumb  vehicle  and  2 s for  a smart  vehicle.  Conse- 
quently, at  either  headway  there  is  approximately  a 
twofold  reduction  in  the  required  sampling  rate  when 
a smart  vehicle  configuration  is  used. 

The  question  of  setting  quantization  requirements 
becomes  more  complex  because  of  the  many  factors 
that  influence  the  resulting  noise  levels.  To  determine 
the  number  of  bits  required  for  on-board  computation 
and  wayside  information,  a statistical  analysis  is  used 
with  quantization  noise  added  at  selected  points  in 
the  linear  model.  A covariance  propagation  equation 
is  derived,  thus  providing  the  variance  in  vehicle 
acceleration  as  the  performance  measure.  The  param- 
eters considered  are  the  on-board  bit  number,  bit 
number  for  D/A  conversion,  wayside  bit  number, 
wayside  sampling  interval,  and  on-board  sampling 
interval. 

A linearized  model  of  the  closed-loop  system  with 
discrete  vehicle  dynamics  may  be  represented  as  in 
Fig.  2.  The  index,  /,  represents  updates  at  the  on- 
board sampling  interval,  Tv.  The  discrete  time  state 
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Fig.  2 General  form  of  cloaed-loop  system. 

equations  are  then  given  by 

*(«  + 1)  - *x(i)  + r,(»)  + B*( «) 
and 

■CO  - + mO) 

where 

♦ = n X n system  matrix, 
x = n-dimensional  state  vector, 
r — n X m matrix, 

t|  = ws-dimensiontl  vector  of  state  noise, 
u - scalar  input, 

M =■  noise  representing  A/D  conversion  noise 
on  wayside  inputs, 

H - 1 X n vector, 

B»  »X  1 vector, 


d 
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and  the  index,  *,  represents  updates  at  the  wayside 
sampling  interval,  Tw.  The  input,  *(»),  is  constant 
over  N samples  of  the  on-board  computer;  Tw 
NTP.  By  evaluating  the  state  at  wayside  sampling 
instants,  an  equation  for  propagation  of  state  covari- 
ance at  wayside  sampling  instants  is  (Ref.  3) 

P(/  + N)  = FPCOF-  + G (l) 


PC  0 = £[*(0*r(0] 

F = * v + Y VBH 

/-  o 

G = + &BM  (&By] 

;-o 

Q = £[iGV(0] 

M = £[/i(»rCO]- 

Consequently,  for  given  Af,  and  sampling  intervals 
Tt,  and  Ttc,  we  first  compute  the  matrices  F and  G. 
For  an  initial  covariance  P(0)  — 0,  the  state  co- 


variance may  be  propagated  using  Eq.  1 to  determine 
a steady-state  noise  level  as  a function  of  various 
quantization  levels  and  sampling  rates  of  the  on-board 
and  wayside  computers. 

With  the  above  technique,  the  trade-offs  between 
various  digital/analog  interfaces  and  sampling  inter- 
vals are  determined.  Based  on  the  results  in  Ref.  3, 
the  selected  configuration  requires  8 bits  for  D/A 
conversion  (placed  at  the  acceleration  command),  8 
bits  for  on-board  computation,  and  a total  of  23  bits 
from  the  wayside.  At  a headway  of  0.3  s,  the  on- 
board sampling  interval  is  set  at  0.02  s;  for  a 3-s 
headway,  it  is  0.03  s. 

The  final  system  problem  considered  is  the  effect 
of  time-delayed  wayside  information.  A modified  z- 
transform  analysis  with  n>ot  locus  plots  is  used  to 
demonstrate  the  rapid  increase  in  system  damping 
with  increased  time  delay.  From  'he  simulation  it  is 
found  that  at  a 0.3-s  headway,  a 0 0'  -s  time  delay 
may  be  tolerated  with  a 0 1-s  wayside  sampling  in- 
terval. At  a 3-s  headway,  a 0.3-s  delay  with  a 1-s 
wayside  sampling  interval  is  acceptable. 
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Fig.  3 Simulation  example. 


All  the  results  have  been  verified  with  a detailed 
computer  simulation  of  the  system.  The  simulation 
includes  quantization  of  wayside  information,  D/A 
conversion  on  board,  and  quantization  at  various 
stages  of  the  on-board  control  computation.  A sepa- 
rately excited  w motor  model  is  used  for  the  pro- 
pulsion system,  including  nonlinear  effects  due  to 
aerodynamic  drag  An  example  of  the  system  response 
with  sampling  and  quantization  is  shown  in  Fig.  ). 
A trailing  vehicle  (dashed  line')  traveling  at  15  m s 
is  overtaking  a preceding  vehicle  (solid  line)  travel- 
ing at  10  m/s  with  an  initial  separation  (tail  to  nose) 
of  6?  m.  The  dashed  line  on  the  separation  plot  rep- 
resents the  kinematically  required  separation  (which 
must  always  be  less  than  the  actual  separation).  At 
20  s,  the  preceding  vehicle  decelerates  on  jerk  and 
acceleration  limits  to  a speed  of  5 m/s.  The  result- 
ing noise  level  is  most  pronounced  on  the  plot  of 
vehicle  jerk  and  is  within  acceptable  bounds. 

Continuing  study  in  this  area  includes  further  trade- 
offs in  terms  of  vehicle- way  side  control  allocation. 
All  allocation  computations  considered  in  this  study 
were  performed  either  on  board  or  at  wayside.  How- 


ever, there  are  alternatives  in  which  some  informa- 
tion may  be  combined  at  wayside  before  being  trans- 
mitted to  the  vehicle.  Another  major  assumption  in 
this  work  was  the  existence  of  ideal  measurements 
of  position,  velocity,  and  acceleration.  Hence,  the 
effects  of  measurement  noise  will  be  studied  to  de- 
termine, for  example,  whether  states  may  be  esti- 
mated from  position  measurements  alone. 
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ACTIVE  ACOUSTIC  DETECTION  OF 
LEAKS  IN  UNDERGROUND  NATURAL-GAS 
DISTRIBUTION  LINES 


The  most  promising  technique  now  being  investi- 
gated for  the  detection  of  leaks  in  buried  natural-gas 
pipelines  is  the  active  acoustic  method.  In  support 
of  an  APL  research  program  to  provide  fundamental 
information  required  in  developing  an  optimal  system 
of  active  acoustic  leak  detection,  a theoretical  and 
experimental  investigation  was  made  of  the  radiation 
accompanying  the  propagation  of  acoustic  waves 
within  an  underground  pipe.  The  investigation  has 
shown  that  the  radiated  field  consists  of  both  shear 
and  compressional  waves.  This  result  is  significant 
because  the  coherent  radiation  from  the  pipe  walls 
represents  the  largest  obstacle  to  reliable  leak  detec- 
tion by  the  active  acoustic  method.  In  a related  devel- 
opment, a sensitive,  non-earth-loading  laser  interferom- 
eter that  can  be  used  In  the  field  has  been  devel- 
oped to  measure  earth  vibrations. 


BACKGROUND 

The  detection  of  leaks  in  residential  natural-gas 
distribution  lines  is  a matter  of  concern  to  both  in- 
dustry and  federal  regulatory  agencies.  A research 
effort  directed  toward  understanding  the  fundamentals 
of  the  active  acoustic  detection  of  leaks  has  been  con- 
ducted at  APL  (Ref.  1).  In  this  promising  method, 
sound  is  injected  into  the  pipe  and  the  leak  is  detected 
by  measuring  displacements  of  the  earth's  surface 
caused  by  radiation  emanating  from  the  leak. 

The  program  at  APL  consists  of:  (a)  an  experi- 
mental field  test  that  makes  use  of  an  abandoned 
water  line  on  the  APL  property  to  study  radiation 
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from  (he  walls  of  nonleaking  and  leaking  pipes.  (b) 
a theoretK.il  investigation  of  elastic  waves  radiated 
from  underground  piping  due  to  vibration  generated 
by  coupling  of  the  pipe  walls  to  the  internal  acoustic 
pressure  variations  and  from  leak-radiated  signals, 
and  (c)  the  development  of  improved  detection  in- 
struments for  measuring  vibrations  of  the  earth's 
surface. 

In  this  article,  we  describe  the  results  of  an  experi- 
mental and  theoretical  study  of  radiation  from  the 
vibrating  walls  of  a buried,  nonleaking  pipe  (Ref.  2). 
This  problem  is  important  to  the  overall  effort  because 
background  radiation  from  the  walls  of  a nonleaking 
pipe  was  the  major  reason  previous  investigators  failed 
to  develop  an  active  acoustic  leak-detection  system  A 
portable  field  instrument  for  detecting  vibrations  of 
the  earth’s  surface  by  optical  means  is  also  described. 


DISCUSSION 

Surface  acceleration  accompanying  the  propagation 
of  plane  progressive  sound  waves  within  an  air-filled 
pipe  was  measured  as  a function  of  frequency  by 
means  of  surface-mounted  accelerometers.  The  pipe, 
a 20-year-old  abandoned  water  line  approximately 
1000  ft  long  with  an  inside  diameter  of  2 in  and  a 
wall  thickness  of  3/16  in.,  is  buried  30  to  36  in. 
deep  m clay  earth.  Before  the  pipe  was  instrumented, 
it  was  dried  by  pumping  and  was  tested  for  leaks; 
none  were  found.  The  acoustic  source  was  a standard 
electromagnetic  unit  of  the  type  used  to  drive  horns 
in  public  address  systems  It  was  impedance- matched 
to  one  end  of  the  pipe.  Surface  accelerations  were 
measures!  by  means  of  a Briiel  and  Kjacr  type  ^698 
accelerometer  located  122  ft  from  the  driver  and 
directly  above  the  pipe.  Checks  were  made  through- 
out the  course  of  the  experiments  to  ensure  that  the 
disturbance  measured  at  the  surface  of  the  earth  was 
due  to  the  transmission  of  acoustic  energy  in  the  gas 
within  the  pipe  and  not  from  the  steel  walls. 

The  results  of  a typical  scries  of  measurements  are 
shown  in  Fig.  1 for  pipe  depths  of  30  and  12  in. 
The  accelerometer  output,  f: ,,  in  microvolts  is  plotted 
against  frequency.  During  these  measurements,  the 
sound  pressure  level  at  the  driver  was  held  constant 
at  approximately  140  dB  relative  to  2 X 10~<  dyne/ 
cm’.  The  drop-off  in  signal  with  frequency  is  caused 
mainly  by  attenuation  in  the  gas  within  the  pipe, 
attenuation  in  the  earth  medium,  variation  in  ac- 
celerometer sensitivity,  and  the  mismatch  of  acceler- 
ometer impedance  relative  to  that  of  the  earth.  The 
last  two  factors  make  it  desirable  to  develop  an  al- 
ternative method  (other  than  accelerometers’)  to  meas- 


Fig. 1 Accelerometer  output  plotted  as  a function  of 
frequency. 


ure  earth  displacements.  The  important  thing  to  note 
in  this  figure  in  the  case  of  the  30-  in.  pipe  depth  is 
the  scries  of  maxima  separated  by  about  400  Hz. 
indicative  of  an  interference  phenomenon. 

The  case  of  an  infinitely  long  steel  pipe  buried  in 
earth  modeled  as  an  isotropic  homogeneous  clastic 
half-space  was  considered  theoretically.  The  gas  within 
the  pipe  conducts  a plane  pressure  wave  that  induces 
radial  and  longitudinal  displacements  of  the  wall. 
These  displacements  excite  compression.il  and  shear 
waves  in  the  earth  medium  that  satisfy  the  wave  equa- 
tions (Ref.  3) 
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giving  the  comprcssional  and  sheer  potentials  of  <3 
and  <fi.  respectively.  Here,  n is  the  comprcssional  wave 
velocity  and  /?  the  sheer  wave  velocity;  t is  time.  These 
velocities  are  related  to  the  elastic  constants  of  the 
earth  medium 
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The  problem  is  solved  by  first  determining  the 
amplitudes  of  the  radiated  cylindrical  waves  satisfying 
Eqs  t and  2 by  matching  the  boundary  conditions 
at  the  pipe  earth  interface.  The  waves  radiate  to  the 
surface  where,  directly  over  the  pipe,  they  are  reflected 
as  bulk  shear  and  compressional  waves,  the  magni- 
tudes of  which  are  determined  by  the  existing 
Ik nindary  conditions  (the  vanishing  of  normal  and 
tangential  stresses).  The  radial  and  longitudinal 
displacements  at  the  earth's  surface  directly  over  the 
pipe  are  given  in  terms  of  the  incident  and  reflected 
compressional  and  shear  potentials. 

Theoretical  vertical  displacement  is  shown  in  Fig. 
2 where  the  displacement,  V,  is  plotted  versus  fre- 
quency for  the  case  of  a pipe  buried  30  in.  deep. 
The  experimental  vertical  displacements  also  plotted 
in  Fig.  2 were  obtained  by  correcting  the  data  in 
Fig.  1 for  the  aforementioned  attenuation  and  ac- 
celerometer sensitivity  and  mismatch  factors  (Ref.  2). 
It  is  noted  from  Fig.  2 that  the  level  of  wall  radia- 
tion from  a buried  pipe  can  be  accounted  for,  within 
a factor  of  two,  by  means  of  theory,  i.e.,  axisymmetric 
vibration  of  the  pipe  walls.  This  result  stands  in 
marked  contrast  to  similar  experiments  with  pipes  in 
an  air  medium.  There,  radiation  levels  exceeded  those 
due  to  axisymmetric  excitation  by  more  than  an  order 
of  magnitude  (Ref.  4)  and  have  been  found  to  be 
generated  by  coupling  into  the  transverse  bending 
modes.  It  is  concluded  that  bending  vibrations  of  a 
buried  pipe  must  be  strongly  damped. 


Fig.  2 Comparison  of  calculated  and  experimental 
values  of  earth  surface  displacement. 


Another  important  result  is  that  the  amplitudes  of 
the  compressional  and  shear  waves  generated  by  the 
pipe  wall  vibration  are  comparable,  leading  to  the 
interference  evidenced  by  the  experimental  data.  Val- 
ues of  the  compressional  and  shear  wave  velocities,  « 
and  p,  were  chosen  to  give  the  best  fit  between  theory 
and  experiment.  To  the  authors'  knowledge,  values 
of  n and  ft  in  situ  were  never  previously  determined 
simultaneously. 

In  order  to  eliminate  problems  associated  with 
surface-mounted  accelerometers,  a Fizeau  laser  inter- 
ferometer (Ref.  1)  was  constructed  and  mounted  on 
a tripod  (Fig.  3).  By  means  of  electronic  stabiliza- 
tion, the  system  was  freed  of  many  of  the  environ- 
mental instabilities  that  commonly  trouble  interferom- 
eters; hence  it  was  suitable  for  field  use.  Extremely 
small  displacements  could  be  detected  with  this  in- 
strument. If  the  displacement  of  the  earth.  A,  is 
represented  by  8 cos  <,>,/.  where  u>,  is  the  angular  fre- 
quency of  sound,  (be  minimum  displacement,  8. 
detected  by  this  system  was  approximately  10  ,J  m 
(0.001  nm),  a dimension  considerably  smaller  than 
an  atom. 


Fig.  3 Field-portable  iaser  interferometer  for  detect- 
ing displacements  of  the  earth's  surface. 
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The  defector  w.ts  used  to  study  surface  or  Rayleigh 
waves  (Ref.  1 ) in  sand  and  earth  in  the  laboratory. 
The  waves  were  generated  by  pulses  from  a piezo- 
electric driver  in  contact  with  the  soil  surface.  The 
Rayleigh  wave  velocity  in  sand,  FV*’,  was  about  *>2 
m s.  In  earth,  Rayleigh  velocity  ranged  from  66  to 
14'  m s,  depending  on  soil  compaction.  The  sand 
measurements  agree  with  results  obtained  by  other 
worker*  ( Ref.  ' ) . 
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STEADY  TEMPERATURE  IN  A 


A group  of  boundary  value  problems  (Ret.  1)  Is 
solved  that  describes  the  steady  state  temperature 
distribution  In  a circular  cylinder  of  infinite  length, 
moving  In  the  direction  of  its  axis.  The  cylinder  is 
heated  by  a hot  gas  flowing  perpendicularly  to  the 
cylinder.  The  theory  supports  experiments  on  the 
steady-state  combustion  of  polymers  carried  out  by 
Hunter,  Fristrom,  and  Grunfelder  (Ref.  2). 


BACKGROUND 

Polymers  arc  commonly  occurring  solids  whose 
molecules  arc  made  up  of  repeating  units;  fabrics 
and  plastics  are  examples.  In  the  experiments  of  Ref. 
2,  a polymer  is  coated  on  a supporting  wire  (or  fiber- 
glass strand)  and  drawn  at  a constant  speed  in  the 
direction  of  its  axis  across  the  wake  of  a (lame,  which 
serves  as  a heat  source.  The  coated  wire  then  reaches 
a steady  state  at  each  position,  fixed  in  laboratory 
coordinates.  These  experiments  have  provided  in- 
sights into  the  fundamental  physics  and  chemistry 
of  polymer  combustion.  The  present  theory  has  aided 
go  in  the  interpretation  of  the  data 
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MOVING  CYLINDER 


DISCUSSION 

The  chemical  reactions  in  polymer  combustion  arc 
controlled  mainly  by  the  temperature  of  the  solid 
phase.  First,  combustion  requires  a mixture  of  flam- 
mable gas  and  oxygen.  The  flammable  gas  evolves 
from  the  surface  as  a result  of  chemical  reactions  in 
which  the  polymer  breaks  down  into  simpler  species, 
some  of  which  are  gaseous  Which  reactions  lake 
place  and  how  fast  they  occur  arc  controlled  mainly 
by  the  temperature  of  the  solid  phase.  In  particular, 
the  onset  of  the  evolution  of  flammable  gas  occurs 
at  a well  defined  surface  temperature  that  is  char- 
acteristic of  the  composition  of  the  polymer.  The 
ignition  delay  time  is  simply  the  warm-up  time  re- 
quired for  the  surface  to  reach  the  flammable  gas 
temperature. 

Reference  I calculates  the  temperature  of  the  mov- 
ing polymer-coated  wire  in  terms  of  experimental 
variables,  including  the  wire  speed  and  diameter, 
the  exposure  time,  the  thermal  properties  of  the 
polymer,  and  the  intensity  of  the  heat  source  as  a 
function  of  longitudinal  position. 


Figure  1 trc.its  .1  simple  case  in  which  the  polymer 
coating  is  very  thin.  The  plot  determines  surface 
temperature.  T,  at  a given  (dimensionless')  exposure 
time  /*.  and  a given  (dimensionless)  heat  transfer 
coefficient,  N:  T,„.  is  the  angle  averaged  temperature 
in  the  gas  near  the  wire,  and  T„  is  the  room  tem- 
perature The  figure  can  he  applied  to  the  determina- 
tion of  the  polymer's  flammable  gas  temperature.  Hie 
value  of  1*  at  ignition  is  measured  when  T„,  and  N 
are  known.  The  flammable  gas  temperature  is  then 
the  corresponding  value  of  T read  from  the  chart. 
This  would  be  difficult  to  measure  directly,  although 
it  is  possible  in  a few  cases  In  those  eases,  the  agree- 
ment between  theory  (refined  somewhat)  and  experi- 
ment is  satisfactory. 

The  inverse  of  big  I,  in  which  exposure  time  is 
plotted  against  surface  temperature,  determines  an 
ignition  delay  time  when  the  flammable  gas  tempera- 
ture is  known  That  kind  of  illustration  has  provided 
correlations  of  data  obtained  in  different  heal  sources 
and  with  different  w ire  sizes 

Currently,  the  results  in  Ref  1 are  being  applied  to 
studies  of  polymer  combustion  at  low  pressure  The 
evolving  gases  are  sampled  for  chemical  analysis,  and 
the  temperature  of  the  surface  where  the  sample  is 
taken  is  then  calculated  from  Ref  I . 
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HEART  AND  CYANIDE  STUDIES 


The  term  "smoke  Inhalation"  has  been  used  to  de- 
scribe tbe  cause  of  death  of  fire  victims  for  many 
years.  The  term  is  imprecise  since  many  physiological 
and  toxicological  Interactions  could  have  caused  the 
death.  The  API  Fire  Casualty  Study  was  initiated  In 
1971  to  clarity  the  machanlsms  involved  whan  a par- 
son dies  as  a result  ot  Insult  to  the  respiratory  sys 


t* 

Fig.  I Surface  temperatures  of  a polymer-coated  wire 
moving  at  constant  speed  in  the  direction  of  its  axis 
through  a cross  flow  of  hot  gas. 
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tern  caused  by  the  toxic  atmosphere  generated  in 
residential  fires. 

Two  segments  of  fhe  ongoing  study  that  have 
provided  an  insight  Info  the  fire  fatality  problem  are 
the  special  heart  and  cyanide  studies.  The  work 
through  the  middle  ot  Calendar  Year  1977  is  dis- 
cussed here,  but  the  studies  are  continuing. 


BACKGROUND 

Smoke  generated  by  a fire  is  complex,  comprising 
both  gaseous  anil  particulate  matter.  Many  constituents 
are  toxic  irritants  to  the  human  body.  This  view  is  at 
variance  with  previous  general  concepts,  which  con- 
sidered carbon  monoxide  produced  by  combustion  to 
be  the  only  toxic  part  of  smoke.  No  serious  concern 
had  been  given  to  the  possibility  of  other  causes  of 
the  toxic  atmosphere  in  a fire.  Even  now,  this  concept 
is  held  so  solidly  that  hospital  emergency  rooms  do 
not  think  of  any  other  possible  problems  when  treat- 
ing smoke  inhalation  victims. 

The  APL  program  includes  detailed  autopsies  of 
fire  victims,  including  blood  toxicology.  The  results 
of  the  initial  work  indicate  that  carbon  monoxide  was 
the  lethal  agent  in  M) % of  the  cases,  but  some  other 
reason  had  to  be  found  for  the  cause  of  death  in  the 
rest.  With  i'ni;  fact  in  mind,  the  study  effort  was  ex- 
panded in  an  attempt  to  determine  whether  preexist- 
ing heart  disease  or  other  toxic  agents  were  involved. 

One  of  the  toxic  agents  considered  was  gaseous 
hydrogen  cyanide.  It  can  be  produced  by  the  combus- 
tion of  nitrogen-containing  materials  such  as  wool, 
nylon,  polyurethane,  and  acrylonitrile,  and  a special 
effort  was  made  to  determine  if  the  fire  victims  had 
been  exposed  to  it. 


DISCUSSION 
Heart  Study 

This  study  (Ref.  1),  initiated  and  carried  out  by 
the  Maryland  State  Office  of  the  Chief  Medical  Exam- 
iner, explored  the  relationship  between  preexisting 
heart  disease,  blood  carbon  monoxide  levels,  and  the 
incidence  of  fire  mortality.  Persons  with  preexisting 
heart  disease  might  be  expected  to  succumb  to  lower 
levels  of  carbon  monoxide,  that  is,  survive  for  a 
shorter  period  of  time  than  healthy  individuals  in  an 
atmosphere  contaminated  with  carbon  monoxide.  This 
premise  assumes  that,  in  an  individual  with  heart 
disease,  the  myocardial  oxygen  supply  is  in  jeopardy 
even  in  a normal  atmosphere. 

The  heart  study  population  of  the  Maryland  Medi- 
cal Examiner's  Office  consisted  of  adult  fire  victims 
(20  years  of  age  and  older)  who  died  during  the 
course  of  a fire  or  within  six  hours  thereafter.  Because 
of  the  diversity  of  the  group  satisfying  these  criteria, 
comparative  analysis  required  the  group  to  be  nar- 
rowed down  to  a more  homogeneous  population.  Also 
in  an  effort  to  achieve  some  degree  of  uniformity  re- 
garding the  terminal  environment  and  the  circum- 


stances, only  those  who  died  as  a result  of  a "typical" 
house  fire  were  analyzed  further.  This  excluded  all 
victims  of  flash  fires,  automobile  collision  fires,  etc. 

The  only  readily  measurable  and  quantifiable 
parameter  of  preexisting  heart  disease  is  coronary 
stenosis,  which  is  the  most  important  and  widespread 
form  of  heart  disease.  The  hearts  of  all  victims  in  the 
study  were  examined.  The  coronary  arterial  tree  was 
dissected  from  the  heart  and  decalcified,  and  the 
arteries  were  serially  sectioned  throughout  their 
lengths.  The  degree  of  luminal  stenosis  is  expressed 
as  the  percentage  of  narrowing  with  respect  to  the 
original  lumen.  The  arteries  were  divided  into  three 
segments:  proximal,  middle,  and  distal,  based  on  the 
total  length  of  each  artery.  The  maximum  stenosis  for 
each  segment  of  each  major  coronary  artery  was  used 
in  conjunction  with  the  other  autopsy  findings  to 
quantify  the  extent  of  heart  disease  in  the  victims. 

One  of  the  early  and  striking  findings  was  that  the 
amount  of  heart  disease,  expressed  in  terms  of  coro- 
nary artery  stenosis,  is  higher  in  the  fire  victims  than 
would  be  expected  in  the  general  population,  particu- 
larly in  the  younger  group.  Of  the  -II  victims  in  the 
young  group  (20  to  39  years  of  age),  eight  had 
greater  than  90%  stenosis  in  at  least  one  segment  of 
the  major  coronary  arteries.  In  addition.  40%  of  the 
overall  total  had  at  least  one  fonis  of  stenosis  of  90% 
or  more,  and  more  than  half  had  narrowing  greater 
than  75%. 

Despite  the  fact  that  other  toxic  gases  and  sub- 
stances may  be  generated  during  a fire,  carbon  mon- 
oxide poisoning  remains  the  primary  killer.  The  gen- 
erally accepted  minimal  lethal  level  for  carbon 
monoxide  intoxication  is  50%  blood  saturation. 
Whenever  death  occurs  with  a lower  level,  other  con- 
tributory factors  must  be  sought.  The  original  analysis 
focused  on  individuals  who  succumbed  with  less  than 
50%  blood  carbon  monoxide  saturation. 

If  the  presence  of  significant  heart  disease  leads  to 
a greater  susceptibility  to  carbon  monoxide  poisoning, 
one  would  expect  a considerable  number  of  these  vic- 
tims to  die  with  a low  level  of  blood  carbon  mon- 
oxide. The  relationship  between  blood  carbon  mon- 
oxide level  and  coronary  artery  stenosis  was  examined. 
No  constant  relationship  between  these  two  parameters 
was  found.  Of  the  victims  with  significant  heart 
disease,  a preponderance  achieved  a blood  saturation 
level  of  greater  than  50%  before  dying. 

Why  did  such  a substantial  number  not  succumb 
to  low  carbon  monoxide  levels?  It  could  be  postulated 
that  preexisting  heart  disease  makes  escape  from  the 
fire  difficult  but  does  not  necessarily  contribute  to  an 


early  death.  The  incapacitated  group  would  then  in- 
hale variable  amounts  of  carbon  monoxide,  accounting 
for  the  wide  range  of  lethal  levels  found.  This 
cardiac-based  incapacity  would  also  explain  the  high 
incidence  of  heart  disease  in  our  series  since  the  heart 
disease  would  ''select”  fatalities.  However,  before 
accepting  this  explanation,  another  variable  had  to  be 
considered. 

Alcohol  has  long  been  known  to  be  a factor  in 
house  fires,  and  various  theories  have  been  pro- 
pounded regarding  the  interaction  of  carbon  mon- 
oxide and  alcohol  as  a cause  of  death.  Blood  alcohol 
determinations  are  made  for  all  fire  victims,  and  this 
parameter  could  readily  be  considered.  Of  the  110 
victims,  92  had  measurable  amounts  of  alcohol  in  their 
blood  at  the  time  of  death,  including  67  of  the  85 
males.  There  is  no  evidence  to  suggest  that  the  pres- 
ence of  alcohol  per  re  has  an  effect  on  the  level  of 
carbon  monoxide  associated  with  death.  In  addition, 
although  both  are  central  nervous  system  depressants, 
there  is  no  apparent  additive  effect. 

Because  the  group  was  relatively  small,  only  the 
males  were  sufficient  in  number  to  warrant  further 
analysis.  They  were  divided  into  "sober"  and  "in- 
ebriated” groups,  the  inebriated  having  a blood  alco- 
hol content  of  0.1  5%  or  higher.  Of  the  sober  individ- 
uals, approximately  half  of  those  with  significant  heart 
disease  died  with  a low  level  of  carbon  monoxide;  the 
remainder  achieved  a level  of  greater  than  50%  be- 
fore dying.  In  the  inebriates!  group,  18  of  21  individ- 
uals with  significant  heart  disease  had  a carbon  mon- 
oxide level  of  greater  than  50%.  When  these  values 
arc  averages!  for  the  two  groups,  we  find  that  sober 
individuals  with  heart  disease  had  a mean  carbon 
monoxide  level  of  51%,  while  the  inebriated  victims 
had  a mean  level  of  64%.  Those  without  significant 
heart  disease  in  Kith  groups  had  similar  levels,  6}  and 
62%  in  the  sober  and  inebriated  groups,  respectively. 

Apparently  people  with  significant  heart  disease 
live  longer  in  a fire  when  unsler  the  influence  of  alco- 
hol than  when  sober.  There  are  several  possible  ex- 
planations, the  simplest  being  that  the  test  group  was 
too  small  for  the  results  to  be  valid.  Possibly  alcohol 
so  impairs  judgment  that  the  inebriated  victims  are 
unaware  of  the  danger  in  the  situation  and  do  not 
panic,  whereas  some  of  the  sober  individuals  die  pre- 
maturely because  of  stress.  However,  if  this  were  the 
case,  one  would  expect  the  inebriated  individuals  with- 
out heart  disease  to  have  higher  levels  of  carbon  mon- 
oxide than  the  sober  ones  and  this  is  not  so.  A more 
interesting  speculation  is  that  the  inebriated  individ- 
uals with  heart  disease  survive  longer  because  alcohol 
acts  as  a dilator  of  atherosclerotic  coronary  arteries. 


Hydrogen  Cyanide  Study 

The  potential  for  the  involvement  of  cyanide  in 
fires  has  increased  with  the  increasing  use  of  poly- 
meric materials  in  residences.  An  extensive  study  was 
undertaken  by  the  Maryland  State  Office  of  the  Chief 
Medical  Examiner  and  APL  to  determine  the  possible 
role  that  cyanide  might  have  in  causing  fire  casualties. 

An  analytical  method  that  would  be  accurate,  re- 
peatable, and  acceptable  to  the  toxicological  commu- 
nity had  to  be  defined.  Suggested  methodologies  were 
evaluated  (Ref.  2),  and  a modification  of  those  sug- 
gested in  Ref.  3 was  selected.  It  is  specific  and  sensi- 
tive to  0.01  gg/ml  of  cyanide. 

The  cyanide  molecule  is  elusive  in  biological  ma- 
terials. It  can  be  lost  and  it  can  be  generated  after 
death.  The  stability  of  cyanide  in  blood  stored  under 
various  conditions  after  death  was  studied  extensively 
in  fire  death  victims  and  in  simulated  blood  samples 
(Ref.  4).  To  minimize  difficulties,  the  blood  for 
analysis  should  be  obtained  within  24  hours  after 
death. 

A set  of  normal  values  was  obtained  by  measuring 
blood  cyanide  levels  in  32  living  subjects,  64  random 
nonfire  victims,  and  22  airplane  crash  victims  who 
died  of  impact  injuries.  The  range  of  readings  for 
the  118  cases  was  0.00  to  0.22  gg/ml;  therefore, 
values  ranging  from  0.00  to  0.25  /ag/ml  were  con- 
sidered normal. 

Cyanide  concentrations  were  determined  in  256  fire 
fatalities  that  occurred  in  Maryland  during  the  period 
1 January  1973  to  28  June  1977  (Refs.  2 and  4).  The 
distribution  of  the  results  of  these  analyses  is  shown 
in  Table  1 . Normal  concentrations  were  found  in 
30%  of  the  cases.  Subtoxic  concentrations  (0.26  to 
1.00  /ig/ml)  were  found  in  33%  and  would  not  be 
expected  to  elicit  biological  effects,  although  concen- 
trations exceeding  0.26  gg/ml  are  indicative  of  the 
presence  of  cyanide  precursors  in  building  materials. 
Possibly  toxic  concentrations  (1. 01  to  2.00  gg/ml) 

TABLE  1 

DISTRIBUTION  OF  BLOOD  CYANIDE 
CONCENTRATIONS  IN  256  FIRE 
FATALITIES  DURING  1 JANUARY  1975  TO 
28  JUNE  1977 


Cyanide 

Concentration 

Kg/ml) 

Number 

Percent 

0.00-0.25 

77 

30 

0.26-1.00 

89 

35 

1.01-2.00 

66 

26 

2 01  and  greater 

24 

9 

256 

100 

T 


were  found  in  26%  of  the  cases,  and  probably  toxic 
concentrations  (2.01  /tg/ml  and  greater)  in  9%. 
These  concentrations  are  indicative  of  exposure  and 
toxicity  and  may  result  in  interaction  with  other  toxic 
substances  and  possible  incapacitation. 

The  detection  of  cyanide  in  concentrations  exceed- 
ing 0.26  fifi/ml  in  70%  of  the  fire  Vitalities  is  a clear 
indication  of  the  widespread  distribution  of  nitrogen- 
containing  polymeric  materials  in  the  average  building 
or  vehicle. 

Blood  cyanide  concentrations  correlated  significantly 
with  carbon  monoxide  saturation  levels.  When  cyanide 
was  high,  carbon  monoxide  was  high;  when  cyanide 
was  low,  carbon  monoxide  was  low.  Apparently, 
cyanide  did  not  explain  deaths  not  otherwise  ex- 
plained. However,  its  role  as  an  intermediary  agent 
(causing  incapacitation  and  subjecting  a victim  to  a 
prolonged  exposure  to  carbon  monoxide)  has  not  been 
fully  elucidated. 

An  automobile  fire  fatality  was  studied  in  which 
the  vehicle  contained  a padded  dashboard  and  plastic- 
covered  seats  filled  with  urethane  foam.  The  blood 
cyanide  concentration  was  3.13  fig/ml,  and  the  blood 
was  18%  saturated  with  carbon  monoxide.  The  blood 
alcohol  concentration  was  0.22%  and  the  blood 


methadone  W’as  0.03  gg/ml;  these  factors  were  con- 
tributory to  the  accident.  Limited  chest  and  bone  in- 
juries and  the  carbon  monoxide  were  insufficient  to 
cause  death,  but  the  cyanide  concentration  was  toxic 
and  hence  was  the  primary  cause.  The  study  of  future 
cases  will  probably  reveal  similar  situations  that  point 
clearly  to  the  toxic  role  of  cyanide.  A similar  case  was 
recently  observed  in  Virginia;  the  saturation  level  of 
blood  carbon  monoxide  was  6%  and  the  cyanide  con- 
centration was  3.1  pg/ml.  The  fire  was  of  the 
smoldering  type  and  no  other  factors  were  evident. 
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OCEAN  SCIENCE  AND  TECHNOLOGY 


INTRODUCTION 


An  understanding  of  the  physics  of  the  ocean  environment  is  essential  to  APL 
in  performing  its  mission  in  support  of  the  Fleet.  This  area  of  physics  includes 
propagation  of  electromagnetic  and  acoustic  radiation;  absorption  and  scat- 
tering both  within  and  at  the  surface  of  the  ocean;  hydrodynamic  phenomena 
such  as  surface  waves,  internal  waves,  and  currents;  and  the  effects  of  physical 
variables  such  as  water  temperature,  pressure,  salinity,  and  density  on  ocean 
dynamic  phenomena.  Central  to  APL's  effort  is  the  development  of  sensing 
and  data-processing  systems  to  detect  and  characterize  signals  propagated 
through  or  at  the  surface  of  the  ocean.  While  analysis,  simulation,  and  labora- 
tory research  are  extremely  useful,  many  phenomena  of  interest  must  be 
studied  at  sea.  Accordingly,  APL  has  developed  the  capability  to  conduct 
large-scale  scientific  studies  in  the  open  ocean. 

Much  of  the  ocean  research  is  not  reported  here  for  reasons  of  national 
security,  but  a small  sample  has  been  selected  to  illustrate  different  areas  of 
interest.  Theoretical  studies  are  exemplified  by  a mathematical  analysis  of  the 
transient  effects  produced  by  the  impulsive  injection  of  fluid  into  a medium 
having  a vertical  density  gradient,  such  as  the  ocean.  The  Navigation  Infor- 
mation and  Computing  System,  described  in  the  next  article,  provides,  for 
the  first  time,  an  automatic  recording  of  ship  position  reference  data  and  a 
real-time  readout  of  ship  position.  As  oceanographic  research  becomes  more 
sophisticated,  greater  precision  and  accuracy  are  required  in  determining  ship 
position. 

Also  reported  are  the  first  successful  measurements  of  the  earth's  gravita- 
tional field  from  a fixed-wing  aircraft.  The  final  two  articles  describe  work 
related  to  SEASAT-A,  a NASA-sponsored  satellite  and  the  first  dedicated  to 
oceanographic  research.  One  article  describes  a wideband  data  link  for  the 
transmission  of  synthetic  aperture  radar  signals  to  the  ground.  The  other  re- 
ports the  development  of  a device  to  simulate  the  radar  return  signal  for  the 
radar  altimeter  aboard  SEASAT-A.  The  simulator  was  used  to  test  the  altim- 
eter’s capability  to  measure  ocean  wave  height  and  satellite  altitude. 

Although  APL’s  increasing  expertise  in  ocean  physics  to  date  has  been  ap- 
plied to  Navy  problems,  it  should  be  available  in  the  future  to  enhance  the 
specific  understanding  of  the  ocean  environment  and  to  assist  in  its  peaceful 
exploitation.  Site  surveys  for  determining  optimum  locations  for  ocean  thermal 
energy  plant-ships  is  one  example  of  many  possible  applications. 


SINGULAR  SOLUTIONS  IN 
BUOYANT  FLOWS 


The  fundamental  solution  for  a system  of  equations 
that  describes  a model  of  buoyant  flow  represents  an 
impulsive  injection  of  fluid  at  a point.  Examination  of 
the  velocity  and  displacement  fields  associated  with 
this  solution  suggests  that  transient  effects  in  general 
are  likely  to  remain  strongly  localized  near  the  region 
where  they  are  generated  and  to  last  a long  time. 


BACKGROUND 

In  the  flow  of  a fluid  of  variable  density,  such  as 
the  ocean,  buoyant  effects  induced  by  gravity  acting 
on  the  fluid  can  play  a significant  role.  Among  these 
effects  are  the  well  known  slug  flows  at  very  low 
velocities,  the  pronounced  oscillatory  flows  induced 
near  the  buoyancy  frequency,  and  the  steady-state 
flows  caused  by  bodies  moving  at  constant  velocity. 

Recently  we  have  been  able  to  find  the  funda- 
mental solution  for  a system  of  equations  that  de- 
scribes the  transient  flow  of  such  a fluid  in  three 
dimensions.  Our  study  (Ref.  1)  of  the  fundamental 
solution,  which  corresponds  physically  to  an  impul- 
sive injection  of  a unit  mass  of  fluid  at  a point, 
graphically  illustrates  the  fact  that  in  this  model  of 
fluid  behavior  the  transient  motion  remains  largely 
localized  in  the  vicinity  of  the  original  disturbance, 
does  not  propagate  a way,  and  lasts  a long  time. 

This  suggests  that  it  would  be  unwise  to  use  steady- 
state  results  for  fluids  that  might  be  modeled  similarly 
if  transient  effects  are  likely  to  be  at  all  significant, 
even  though  in  the  real  world  the  amplitudes  would 
eventually  dissipate  because  of  viscous  and  nonlinear 
effects. 
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DISCUSSION 

The  usual  conservation  laws  for  an  inviscid  incom- 
pressible fluid  with  an  initial  stable  unperturbed  den- 
sity, p„,  that  depends  only  on  the  vertical  coordinate, 
x3,  yield,  when  linearized  about  the  static  solution, 
a system  of  equations  given  by 


d2w  j , «’  dp 

— — -f  n2w  H — = 0 

dt2  Po  5X3 


and 


d2<p  dV  , d 1 dip  dw 

dxi2  dxj!  dx»  po  dxj  P°  dxj 


to  be  satisfied  at  every  point  away  from  regions 
where  there  are  sources  or  body  forces  other  than 
gravity.  The  quantity  n2  stands  for  — po)dpo! dx3. 
The  velocities,  which  arc  the  most  important 
quantities  associated  with  a solution  of  these 
equations,  are  given  for  the  horizontal  components, 
U\  and  «2,  by 


1 dip 
Po  OXa 


a = 1,  2 


and  for  the  vertical  component  by 

, 1 dip  . 

ui  = w -\ — 

Po  dxj 


When  n and  po  are  treated  as  constants  in  these 
equations  (a  frequently  made  approximation),  the 
fundamental  solution  discussed  in  Ref.  2 can  be 
obtained  explicitly.  Analytical  details  of  the  solution 
are  given  in  Refs.  2 and  3.  Here  we  note  that  the 
velocities  are  purely  radial  into  and  out  from  the 
disturbance,  and  that  they  can  be  integrated  exactly 
to  obtain  the  radial  displacement  field. 

The  impulsive  source  for  which  we  show  velocity 
and  displacement  graphs  emits  a unit  sphere  of  fluid 
at  t = 0.  Initially  the  fluid  is  pushed  out  radially  to 
accommodate  the  injected  fluid.  The  subsequent  flow 
develops  slowly.  Above  and  below  the  disturbance 
the  fluid  velocity  is  vertical  and  oscillates  sinusoidally 
at  the  buoyancy  frequency.  Horizontally,  the  fluid  ve- 
locity is  steadily  outward  from  the  origin  and  tends 
in  time  to  the  value  l/(3r2),  where  r is  the  hori- 
zontal distance  from  the  origin. 

The  first  series  of  graphs  (Fig.  1)  shows  the 
velocity  field  development.  The  flow  is  axially  sym- 
metric about  the  vertical  axis  and  symmetric  with 
respect  to  the  horizontal  plane  through  the  origin; 
hence  one  quadrant  suffices  to  show  everything.  Since 
the  velocities  are  radial,  they  are  plotted  radially;  the 
scale  is  such  that  the  maximum  velocity  above  the 
origin  is  of  magnitude  one.  For  visual  reference,  the 
0 and  ±5  values  of  the  velocities  are  shown. 

The  six  graphs  of  Fig.  1 show  the  development 
of  the  velocity  field  at  quarter-period  intervals  through 
the  first  one  and  one-half  buoyancy  periods.  Its 
further  development  is  shown  at  each  period  from 
the  second  through  the  tenth  in  Fig.  2.  Significant 
in  these  graphs  are  the  steady  oscillation  at  the  buoy- 
ancy frequency  above  the  disturbance,  the  develop- 
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Fig.  1 Impulse  source  velocities  at  quarter-period  in- 
tervals from  0.25  to  1.5  periods. 
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Fig.  2 Impulse  source  velocities  at  one-period  inter- 
vals from  2 to  10  periods. 


mcnt  of  convolutions  or  waves  in  the  velocity  as 
time  passes,  and  the  rapid  approach  to  a steady  out- 
ward flow  in  the  horizontal  plane.  It  is  also  evident 
that  the  vertical  oscillations  above  and  below  the 
source  propagate  toward  the  horizontal  plane  to  give 
a convoluted  and  wavy  velocity  field  as  time  passes. 

The  next  series  of  graphs  shows  the  displacements 
that  result  from  the  same  impulsive  source.  Figure  3 


shows  what  happens  to  spherical  shells  at  quarter- 
period  intervals  to  one  and  one-half  periods.  The 
shells  are  originally  at  the  positions  shown  in  the 
first  graph  of  Fig.  3.  The  second  graph  of  Fig.  3 
shows  the  radial  displacements  of  these  shells  imme- 
diately following  the  impulse.  Figure  1 shows  the 
subsequent  development  of  the  displacements  at  each 
period  through  the  tenth.  The  displacements,  like 


the  velocities,  also  oscillate  steadily  above  the  source, 
become  convoluted  with  waves  that  move  toward 
the  horizontal,  grow  steadily  outward  in  the  hori- 
zontal plane,  but  do  not  propagate  away  from  the 
source. 
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NAVIGATION  INFORMATION  AND 
COMPUTING  SYSTEM 


The  Navigation  Information  and  Computing  System 
(NICS)  was  developed  to  provide  a real  time  display 
and  a permanent  record  of  a submarine's  position  as 
well  as  to  predict  the  submarine's  location  at  the  time 
of  missile  launch.  NICS  can  perform  these  functions 
using  one  or  several  position  measurement  systems, 
selecting  or  combining  the  outputs  of  each  system  as 
directed  by  the  operator. 


BACKGROUND 

Evaluation  of  the  performance  of  submarine  navi- 
gation systems  requires  accurate  measurement  of  the 
submarine's  position  while  at  sea.  Recent  advances 
in  digital  electronics  and  signal  processing,  and  an 
increasing  demand  for  the  precise  determination  of 
ship  position  for  such  purposes  as  offshore  oil  ex- 
ploration, have  contributed  to  the  development  and 
marketing  of  several  new  position  reference  systems. 
To  determine  the  suitability  of  these  systems  for  the 
evaluation  of  navigation  systems,  the  Navy's  Strategic 
Systems  Projects  Office  requested  APL  in  1975  to 
conduct  a field  demonstration,  comparison,  and  per- 
formance evaluation  of  several  position  reference 
systems  (Refs.  1 and  2).  A special-purpose  interface 
device  was  designed  and  built  by  APL  to  provide 
time  coordination  and  simultaneous  data  recording 
for  as  many  as  six  systems  (Fig.  1).  The  reliable 
performance  and  convenient  operation  of  the  inter- 
face and  control  unit  (IFCU)  motivated  further  de- 
velopment and  the  subsequent  production  of  two 
NICS  units. 


Fig.  1 Interface  and  control  unit  configuration  for 
demonstration  of  position  reference.  Six  position  refer- 
ence systems  were  interfaced  to  the  tape  recorder  by 
the  predecessor  of  the  NICS. 


DISCUSSION 

NICS  (Fig.  2)  comprises  several  commercially 
available  hardware  components  as  well  as  the  spe- 
cially designed  IFCU  (Ref.  3).  Included  are  a 
Hewlett-Packard  9825  programmable  calculator  with 
thermal  printer  and  X-Y  plotter  peripherals,  a 
cartridge-type  tape  recorder,  and  a graphics  unit  with 
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a cathode  ray  tube  (CRT)  display.  The  IFCU  func- 
tions as  the  interface  between  the  several  position 
reference  systems  and  the  calculating,  displaying,  and 
recording  equipment  (see  Fig.  3).  It  also  contains  the 
system  clock  and  an  intervalometer,  which  controls  the 
transfer  of  data  to  the  tape  recorder. 


Fig.  2 The  Navigation  Information  and  Computing 
System,  shown  with  Lorac  hyperbolic  equipment. 


The  data-recording  function  of  NICS  is  independ- 
ent of  calculator  operation.  Once  each  second,  the 
output  of  each  position  measurement  system  is  sam- 
pled and  stored  in  a buffer.  When  the  buffer  has  been 
filled,  the  calculator  input  output  status  is  checked. 
If  the  calculator  is  ready  to  accept  data,  the  con- 
tents of  the  buffer  are  transferred  to  the  IFCU  data 
bus  under  calculator  control.  If  the  calculator  is  not 
ready  for  input  or  is  not  present,  its  control  lines 
are  disabled  and  a built-in  counter  transfers  the 
buffer  contents  to  the  IFCU  bus  If  the  intervalom- 
eter has  counted  down  to  zero,  a buffer  for  the  tape 
recorder  is  enabled  and  accepts  the  data  on  the  IFCU 
bus  for  subsequent  transfer  to  the  tape  recorder. 

Data  are  transferred  from  the  calculator  to  the 
thermal  printer  and  the  X-Y  plotter  via  the  manu- 
facturer's interface.  API.  devised  a means  of  display- 
ing and  plotting  processed  data  on  the  CRT  using 
the  Megatek  BP-741  graphics  display  interface. 

The  primary  application  of  NICS  (Ref.  4)  is  in 
Polaris  and  Poseidon  submarine  Demonstration  and 
Shakedown  Operations  conducted  near  Cape  Canav- 
eral, Florida.  The  position  reference  system  currently 
used  is  Lorac  Hyperbolic.  Lorac  data  in  the  form 
of  Red  Lane  count  and  Green  Lane  count  are  trans- 
ferred to  the  calculator  every  30  seconds.  These  data 
are  converted  to  latitude  and  longitude  by  the  cal- 
culator and  are  output  to  the  printer. 

The  CRT  displays,  located  in  the  submarine  navi- 
gation center  and  in  the  control  room,  show  a rep- 
resentation of  the  nearby  Florida  coast,  operating 
area  boundaries,  and  the  submarine's  position.  The 
calculator  differentiates  the  calculated  position  and 
maintains  an  estimate  of  speed  over  ground  and 
course.  This  velocity  estimate  is  then  used  to  calcu- 
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Fig.  .3  Block  diagram  of  the  Navigation  Information  and  Computing  System. 


late  the  effects  of  (Kean  current  on  ship  movement 
and  to  predict  the  coordinates  of  the  submarine  posi- 
tion at  the  time  of  missile  launch.  Lorac  lane  count 
and  time  are  recorded  on  magnetic  tape  and  are  con- 
verted to  geodetic  position  for  post-mission  evalua- 
tion of  the  navigation  systems. 
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AIRBORNE  GRAVITY  MEASUREMENTS 


Gravity  measurements  made  during  Project  Magnet 
flight  tests  proved  the  feasibility  of  airborne  gravity 
surveys.  Direct  data  processing  provided  gravity  pro- 
files accurate  to  five  parts  per  million. 


BACKGROUND 

In  November  1972,  APL  was  funded  to  study 
the  feasibility  of  using  the  Project  Magnet  airplane 
to  conduct  gravity  surveys.  The  APL  Gravity  Study 
(Ref.  1)  concluded  that  data  from  precision  navi- 
gation and  altimetry  instruments  would  be  needed 
to  process  data  from  a gravity  meter.  A precision 
navigation  system  was  installed  in  the  aircraft  in 
1974.  Precision  altimetry  instruments  and  a gravity 
meter  were  installed  in  1976.  During  a test  flight 
on  I November  1976,  collected  data  proved  that 
precision  gravity  measurements  could  be  made  from 
a fixed  wing  aircraft  (Ref.  2). 


DISCUSSION 

Gravity  Measurfments  The  APL  Gravity  Study 
concluded  that  airborne  gravity  surveys  could  pro- 
vide accuracies  better  than  ten  parts  per  million.  The 
cgs  unit  of  acceleration,  the  gal,  is  1 cm/s*.  The 
acceleration  of  the  earth's  gravity  is  approximately 
980  gal  This  value  wrics  by  ±S  gal  over  the  sur- 


face of  the  earth.  Thus,  a gravity  survey  accurate  to 
ten  parts  per  million  means  that  the  processed  data 
are  accurate  to  ± 10  mgal. 

The  LaCoste  and  Romberg  gravity  meter  measures 
the  force  required  to  balance  a mass  against  the 
force  of  gravity.  Two  levels  of  force  are  established: 
a constant  force  of  approximately  980  gal  to  cancel 
the  bulk  of  the  earth  s gravity  force,  and  a dynamic 
force  to  balance  variations  in  the  ambient  gravity 
field.  In  an  airborne  survey,  motions  of  the  airplane 
cause  variations  in  the  force  required  to  balance  the 
mass. 

The  final  product  of  a gravity  survey  is  an  area 
map  depicting  the  gravity  anomaly.  The  anomaly  is 
the  difference  between  measured  and  computed  grav- 
ity values.  Thr  computed  values  are  derived  from  a 
mathematical  model  of  the  earth's  gravity  field.  In 
order  to  make  an  accurate  comparison  between  the 
two  gravity  values,  the  exact  position  of  the  airplane 
must  be  known.  For  instance,  in  an  area  where  the 
gravity  changes  10  mgal/mi,  a navigation  error  of 
1 mi  could  cause  a 10-mgal  error. 

Gravity  surveys  are  conducted  at  constant  altitude; 
an  airplane  flying  at  constant  altitude  over  the  surface 
of  the  rotating  earth  experiences  centrifugal  accelera- 
tions. The  acceleration  due  to  an  east  or  west  velocity’ 
component  is  approximately  7 mgal  /kt. 


The  altitude  of  the  airplane  affects  the  gravity  data 
in  two  ways.  At  survey  altitudes,  the  acceleration  of 
gravity  decreases  approximately  1 mgal  for  each  10 
ft  of  altitude.  Also,  vertical  motions  of  the  airplane 
cause  forces  on  the  mass  in  the  gravity  meter,  hut 
the  gravity  meter  cannot  detect  the  difference  be- 
tween the  force  due  to  gravity  and  the  forces  due 
to  airplane  motion. 

Navigation  System.  The  Honeywell  electrically 
suspended  gyro  (ESG1  navigation  system  enables  the 
Magnet  airplane  to  navigate  surveys  with  the  accuracy 
required  for  gravity  measurements.  The  accuracy  of 
the  ESG  is  specified  to  be  0.1  nmi/h.  This  means 
that  the  ESG  navigation  data  will  be  accurate  to 
within  1 nmi  during  a 10-h  survey. 

A conventional  "local  vertical"  system  attempts  to 
maintain  a level  platform  in  order  to  align  the  hori- 
zontal accelerometers  on  the  platform  perpendicular 
to  the  force  of  gravity.  Thus,  the  vertical  accelerom- 
eter is  aligned  parallel  with  the  local  force  of  grav- 
ity. The  gyros  must  be  torqued  constantly  to  track 
the  local  vertical. 

The  ESG  is  a space-stable  inertial  system.  It  is 
aligned  initially  in  the  local  vertical  mode.  However, 
once  aligned,  the  gyros  are  left  to  spin  in  the  initial 
alignment.  Thus,  the  accelerometers  on  the  ESG  plat- 
form measure  the  accelerations  of  the  airplane  in 
space-stable  coordinates.  A computer  transforms  the 
ESG  gyro  data  to  values  computed  for  local  vertical 
coordinates.  The  absence  of  perturbing  torques  re- 
duces system  errors  by  a factor  of  ten  or  better. 

Altimetry  Systems.  The  Honeywell  precision 
radar  altimeter  enables  the  Magnet  system  to  deter- 
mine the  altitude  of  the  airplane  over  water.  This 
instrument  provides  altitude  data  that  are  accurate  to 
10  ft.  In  operation,  a burst  of  radar  pulses  is  trans- 
mitted by  the  altimeter,  which  then  pauses  to  moni- 
tor the  returning  signals.  When  "locked  on,”  the 
radar  circuits  detect  and  track  the  first  signals  return- 
ing from  each  radar  burst.  Timing  circuitry  computes 
the  transit  time  between  transmission  and  return,  and 
this  transit  time  is  used  to  determine  the  altitude  of 
the  airplane. 

The  Rosemount  pressure  port  calibrator  compares 
the  atmospheric  pressure  with  a reference  pressure. 
At  survey  altitude,  a test  valve  is  opened  to  fill  a 
reference  chamber  at  ambient  pressure  and  is  then 
closed  to  trap  the  reference  pressure  in  the  chamber 
During  a survey,  changes  of  less  than  1 ft  in  alti- 
tude can  be  detected. 

Gravity  Survey.  On  . November  1976,  a gravity 
survey  was  conducted  at  16  000  ft  to  escape  unusual 


turbulence  at  lower  altitudes.  The  airsneed  was  ap- 
proximately 2S0  kt.  The  gravity  data  were  recorded 
on  magnetic  tape  by  the  LaCoste  and  Romberg  Air- 
Sea  Gravity  System,  and  the  navigation  data  were 
recorded  on  magnetic  tape  by  the  Magnet  system. 

The  results  of  the  gravity  survey  were  analyzed,  and 
gravity  measurements  accurate  to  5 mgal  were  deter- 
mined. Computations  were  made  from  data  samples 
taken  at  10-s  intervals.  Figure  1 (from  Ref.  1) 
shows  plots  of  the  gravity  anomaly  for  one  track  of 
the  survey.  This  track  was  flown  in  an  area  "or  which 
the  gravity  anomaly  is  well  known.  One  curve  shows 
the  gravity  profile  observed  during  a shipboard  sur- 
vey of  the  area.  The  second  curve  shows  the  profile 
extrapolated  to  depict  the  profile  computed  for  an 
altitude  of  16  000  ft.  The  third  curve  shows  the 
profile  obtained  after  subjecting  the  extrapolated  data 
to  a low-pass  filter.  This  curve  was  used  to  evaluate 
the  results  of  the  Project  Magnet  gravity  survey 
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Fig.  1 Effect  of  attenuation  and  filtering  on  a gravity 
profile. 


Figure  2 (from  Ref.  2)  shows  a comparison  of 
the  observes!  air  gravity  profile  with  the  reference 
profile.  Di.  LaCoste  made  the  following  comments 
on  this  figure: 

"The  filtered  attenuated  sea  gravity  profile  is  the 
smoother  of  the  two  curves.  However  the  observed 
air  profile  follows  it  closely.  Since  we  have  made  no 
arbitrary  systematic  corrections  to  the  observed  air 
profile,  its  negligible  average  offset  from  the  sea 
profile  is  significant.  It  indicates  good  performance  of 
all  the  system  components  over  the  longer  periods. 
The  rest  of  the  air  gravity-  results  for  which  we  have 
comparison  sea  gravity  are  about  the  same  as  shown 
in  [this  figure].  Unfortunately,  we  do  not  yet  have 
comparison  sea  gravity  for  about  two-thirds  of  the 
data  we  have  processed.  Such  data  will  be  available 


Fig.  2 Comparison  of  the  gravity  obtained  by  air  with  the  gravity  ob- 
tained at  sea  using  the  same  filtering. 


shortly.  For  the  air  gravity  for  which  we  have  com- 
parison data,  the  noise  or  roughness  in  the  air  gravity 
amounts  to  only  about  4 or  5 milligals,  which  is  as 
good  as  I had  hoped  for  when  I began  analysing 
the  data.  Obviously  the  roughness  can  be  reduced  by 
more  filtering  but  only  with  loss  of  detail.” 

Figure  3 compares  readings  from  the  gravity 
meter  with  data  from  three  altitude  sensors  in  the 
Magnet  airplane.  The  pressure  port  calibrator  gives 
the  best  correlation. 


Six  20-s  filter  stages 


Fig.  3 Comparison  of  gravity  meter  readings  with  ac- 
celerations obtained  from  three  types  of  measurement. 


CONCLUSION 

The  Project  Magnet  gravity  tests  demonstrated  the 
feasibility  of  airborne  gravity  surveys.  Direct  data 
processing  provided  gravity  profiles  accurate  to  5 
mgal.  Estimates  of  system  errors  could  be  used  to 
enhance  these  results. 

The  accuracy  of  the  navigation  data  could  be  im- 
proved by  applying  Kalman  filtering,  which  com- 
bines independent  navigation  measurements  to  pro- 
vide a statistical  estimate  of  system  errors.  Gravity 
surveys  accurate  to  1 mgal  may  be  possible  with  this 
technique. 
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DATA  LINK  FOR  SEASAT-A  SYNTHETIC 
APERTURE  RADAR 


A wideband  anal og  data  link  that  has  been  devel- 
oped by  APL  for  NASA  allows  the  transmission  of  both 
radar-return  and  timing-reference  signals  from  the 
SEASAT-A  synthetic  aperture  radar  (SAR).  Significant 
features  of  the  data  link  include  the  realization  of  a 
linear  (versus  frequency)  20-MHz  analog  channel  and 
the  transfer  of  time  from  the  spacecraft  with  a preci- 
sion of  better  than  1 ns. 


BACKGROUND 

The  SEASAT-A  spacecraft  will  contain  an  experi- 
mental SAR  developed  by  the  Jet  Propulsion  Labora- 
tory of  the  California  Institute  of  Technology.  The 
SAR  will  provide  all-weather  imaging  of  both  land 
and  sea  surfaces  with  a resolution  of  25  m.  As  orig- 
inally conceived,  the  wideband  radar-return  signal 
was  to  be  digitized  at  the  spacecraft  and  the  data 
were  to  be  transmitted  using  one  or  more  standard 
NASA  S-band  telemetry  channels.  However,  a short 
study  by  APL  revealed  that,  for  full  resolution 
imagery,  the  on-board  digitizing  approach  consumed 
more  transmission  bandwidth  than  was  readily  avail- 
able. As  an  alternative,  APL  proposed  an  analog  data 
link  that  uses  simple  frequency  translation  of  the 
I. -band  echo  to  S band,  thereby  minimizing  the  rf 
allocation  problem. 

In  September  1975,  the  SEASAT  project  office 
directed  APL  to  commence  development  of  the  data 
link.  During  the  past  year,  the  work  was  completed 


with  the  delivery  of  i flight  electronics  unit  and  three 
ground  demodulators  for  satellite  integration. 


DISCUSSION 

The  SAR  data  link  (Fig.  1)  consists  of  a space 
segment  and  a ground  segment.  Inputs  to  the  modu- 
lator/transmitter on  the  spacecraft  consist  of  the  radar 
echo  obtained  from  the  SAR  receiver,  the  radar  pulse- 
repetition- frequency  (prf)  timing  signal,  and  a 
sample  of  the  radar  oscillator  signal  called  the 
STALO.  Within  the  modulator/transmitter,  the 
STAI.O  signal  is  used  to  generate  a coherent  beacon 
signal  centered  in  the  middle  of  the  radar  return 
spectrum.  It  also  generates  the  conversion  frequency 
required  to  translate  the  L-band  (1275  MHz)  radar 
output  to  S band  (2265  MHz) . 

The  beacon  signal  is  phase-modulated  by  a 4096- 
bit  pseudorandom  noise  (PRN)  binary  code  that  is 
synchronized  with  the  radar  prf  input;  the  code 
subdivides  the  prf  interval,  nominally  600  gs,  into 
precise  subintervals  of  about  160  ns  and  enables  re- 
covery on  the  ground  to  the  subnanosecond  level. 
The  modulated  signal  is  RF-summed  with  the  beacon 
and  radar  echo  signals,  and  the  composite  is  trans- 
lated to  S band,  amplified  in  a linear  5-W  amplifier, 
and  transmitted  through  a right-hand  circular  an- 
tenna with  (nominally)  2.5  dB  of  gain.  The  data 
link  flight  electronics  package  is  shown  in  Fig.  2. 
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Fig.  1 SEASAT-A  SAR  data  link. 


A 


Fig.  2 Data  link  flight  electronics  (without  cover). 


The  ground  segment  of  the  data  link  comprises 
receiving  equipment  that  is  already  available  in  the 
NASA  Spacecraft  Tracking  and  Data  Network 
(STDN)  stations  and  a demodulator  unit  unique  to 
the  SAR,  developed  as  part  of  the  data  link  by  API.. 
The  composite  SEASAT  signal  is  received  by  the 
STDN  30— ft  dish  antenna  and  then  input  to  the 
receiving  equipment.  The  STDN  receiver  phase  locks 
to  the  data  link  beacon  component  and  thereby  re- 
moves both  the  amplitude  and  phase  variations  caused 
by  the  downlink  slant-range  dynamics.  The  output 
of  the  STDN  receiver  consists  of  the  composite  sig- 
nal centered  at  110  MHz  and  normalized  in  ampli- 
tude; the  only  amplitude  variations  left  on  this  signal 
are  due  to  target  cross  section  and  the  geometry  of  the 
target  link. 

The  function  of  the  SAR  demodulator  unit  is  to 
decompose  the  composite  signal  from  the  STDN 
receiver  and  essentially  regenerate  the  original  input 
signals  to  drive  the  ground-based  recording  equip- 
ment. The  PRF  and  STALO  signals  are  recovered  by 
a delay-locked  loop  that  performs  a real-time  corre- 
lation between  the  received  data  link  prn  code  and 
a stored  replica,  producing  a PRF  output  signal  that 
matches  the  input.  The  delay-locked  loop  also  recov- 
ers the  STAI.O  signal  by  using  a special  reference- 
signal  input  from  the  STDN  receiver.  This  signal 
contains  the  real-time  downlink  doppler  frequency 
required  to  preserve  -.ign.il  coherency. 

The  radar  echo  component  is  recovered  by  mixing 
the  receiver  output  with  a reference  signal  that  has 
also  been  dopplcr-corrected  to  preserve  coherency. 
The  video-echo  output  is  centered  at  11.4  MHz  and 
is  equal  in  amplitude  to  the  input  level  at  the  space- 
craft. The  demodulator  hardware  is  shown  in  Fig.  3. 


Fig.  3 Data  link  demodulator. 


Because  of  its  analog  nature,  the  data  link  can 
be  thought  of  as  an  extension  of  the  SAR  on-board 
receiver.  The  perfect  data  link  would  provide  a unity 
transfer  function  (i.c.,  what  comes  out  is  identical 
to  what  goes  in);  the  practical  data  link  will  exhibit 
errors  due  to  such  effects  as  bandwidth  limitations, 
electronic  nonlinearities,  system  noise,  and  oscillator 
instabilities.  Specifications  for  the  allowable  data  link 
errors  were  developed  early  in  the  program  and  can 
be  grouped  into  two  categories:  those  relating  to  the 
transfer  function,  which  are  measured  by  comparing 
the  video  output  signal  with  the  L-band  input,  and 
those  relating  to  the  recovery  of  the  timing  signals. 
The  latter  ran  be  verified  by  comparing  the  input  PRF 
and  STAI.O  signals  with  the  recovered  versions  in 
the  presence  of  a simulated  radar  return  signal,  which 
is  the  predominant  source  of  timing  jitter. 

Table  1 compares  the  measured  performance  of 
the  data  link  flight  unit,  the  STDN  receiver,  and  the 
SAR  demodulator  connected  as  an  end-to-end  sys- 
tem. Design  specifications  are  also  indicated.  During 
these  tests,  the  antennas  were  replaced  by  a hard- 
wired connection. 

The  first  six  entries  in  Table  1 show  the  perform- 
ance of  the  transfer  function.  The  gain  variation 
occurs  over  a temperature  range  of  —10  to  + 40°C, 
a much  greater  range  than  will  occur  in  orbit.  Over 
the  more  limited  temperature  range,  the  gain  is  well 
within  ±0.5  dB.  Normally,  a linear  system  driven 
beyond  its  linear  range  will  exhibit  a reduction  in 
gain.  This  is  not  the  case  with  the  data  link  because 
the  nonlinearity  of  the  S-band  power  amplifier  causes 
suppression  of  the  data  link  beacon  signal  in  the 
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TABLE  1 

PERFORMANCE  MEASUREMENTS  OF  SAR  DATA  LINK 


Specification 

Requirements 

Measurements 

Gain  at  center  frequency 

0 ± O S dB 

-O  S to  +1.3  dB 

Output  noise  at 

T - 88  K 

< — 17  dBm 

< - IS  dBm  at  T - 420  K 

Intermodualtion  (two 
input  tones  at  — 3 
dBm  each) 

< -20  dBc 

< -23  dBc 

Input  level  for  1 dB 
output  expansion 

> + 5 dBm 

> ' + 3 dBm 

Amplitude  versus  frequency 
± 9.3  MHz 

< 2 dB  peak  quadratic 
error 

< 0.8  dB 

Phase  versus  frequency 
± 9.5  MHz 

< 26°  peak  cubic 
error 

< 10° 

Phase  jitter  of  offset 
video  at  input  of 

0 dBm 

< 10°  rms  for 

T “ 5 ms 

< 7° 

Timing  jitter  of  prf 

< 1 ns  for 

< 0.8  ns  at  0 dBm 

at  input  of  0 dBm 

r - 2 s 

— 1 ns  at  +4  dBm 

0.7  ns  for  time  gate 
on  at  all  input  levels 

leakage 

Fig.  4 Simulated  SAR  return  signal. 


presence  of  a larger  radar  echo;  the  STDN  receiver 
reacts  by  increasing  its  gain  to  overcome  the  sup- 
pression, hut  the  gain  increment  is  also  applied  to 
the  radar  echo  and  leads  to  the  expansion  effect  shown 
in  the  table. 


To  evaluate  the  phase  and  time-recovery  perform- 
ance, a SAR  return  simulator  was  used.  It  makes  use 
of  a white  noise  source,  an  active  chirp  generator, 
and  RF-gating  circuits  to  generate  a repetitive  pulscd- 
RF  signal  (Fig.  4).  This  signal  is  a fairly  good  ap- 
proximation of  the  input  to  the  data  link  in  actual 
operation.  As  the  table  shows,  the  prf  time-recovery 
jitter  exceeds  I ns  for  input  echo  levels  greater  than 
+ 4 dBm;  however,  the  demodulator  has  a special 
gating  feature  that  eliminates  the  leakage  and  echo 
energy  from  the  delay-locked  loop  and  makes  the 
time  jitter  independent  of  the  radar  echo  level. 
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SIMULATION  OF  A RADAR  ALTIMETER 
OCEAN-RETURN  SIGNAL 


Under  the  sponsorship  ot  the  National  Aeronautics 
and  Space  Administration,  API  has  designed  and 
developed  a return  signal  simulator  (RSS).  It  generates 
a signal  with  the  expected  characteristics  ot  the  ocean 
return  for  the  SEASAT-A  radar  altimeter  transmit 
waveform.  With  this  equipment  (Fig.  I),  a means  is 
provided  tor  testing  and  calibrating  the  altimeter  over 
the  entire  range  ot  its  specified  performance  capa- 
bilities. 


BACKGROUND 

The  SFASAT  A radar  altimeter  is  designed  to 
measure  the  distance  between  the  satellite  and  the 
ocean  surface  along  the  nadir-pointing  axis.  For  this 
purpose,  transmitted  microwave  pulses  are  modu- 
lated with  a linear  fm  waveform.  The  processing  of 
the  ocean-surface  reflections  of  these  pulses  yields  pre- 
cision altitude  data  and  provides  a sensitive  measure 
of  the  surface  roughness  or  "sea  state." 

To  appreciate  the  implementation  of  the  return- 
signal  simulator,  a basic  understanding  of  the  nature 
of  the  return  signal  from  the  ocean  and  the  process- 
ing performed  by  the  altimeter  on  this  signal  is 
required.  The  altimeter's  function  is  to  measure  the 
time  between  the  transmission  of  a pulse  and  its 
received  reflection.  The  region  of  linear  rise  of  the 
return  signal  (the  leading  edge)  defines  its  time  of 
arrival.  The  altimeter  uses  a linear  fm  transmit  wave- 
form and  full  deramp  processing  of  the  return  signal. 
Figure  2 depicts  the  full  deramp  process. 

The  altimeter  transmits  a 3.2-ps  pulse  that  has 
a linear  320  MHz  excursion  centered  at  13.3  GHz. 
Reflections  from  the  ocean  surface  over  a wide  range 
arrive  at  the  altimeter  as  a packet  of  time-dispersed 
replicas  of  the  transmit  pulse.  The  altimeter  gen- 
erates a local  oscillator  (I  O)  pulse  that  is  timed  to 
match  the  arrival  of  the  return  signal.  The  charac- 
teristics of  the  1-0  arc  identical  to  those  of  the  trans- 
mit pulse  except  that  they  are  offset  in  frequency  by 
an  amount  equal  to  the  intermediate  frequency  of 
the  altimeter  receiver.  Mixing  of  the  1.0  pulse  with 
the  return  signal  results  in  a translation  of  the  arrival 
of  the  return  signal  from  the  time  to  the  frequency 
domain  (Fig.  2).  A spectral  analysis  of  the  mixer 
output  and  averaging  of  a sufficient  number  of  re- 
flected signals  develop  the  average  return.  The  power 
level  of  a spectral  line  indicates  the  surface  area 
illuminated  at  a given  range  or  time  delay. 


Fig.  1 The  return  signal  simulator. 


Local 


Fig.  2 Full  Jerantp  processing. 
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DISCUSSION 


Instead  of  a packet  of  time-dispersed  replicas  of 
the  transmit  waveform  (as  is  characteristic  of  the 
ocean  reflections  of  the  altimeter  transmit  signal), 
the  return  signal  synthesized  by  the  RSS  is  a packet 
of  frequency -dispersed  pulses  that  are  generated 
simultaneously.  Because  of  the  equivalence  of  fre- 
quency and  time  dispersion  with  the  altimeter’s  full 
deramp  processing  technique,  there  is  no  essential 
difference  between  the  two  types  of  representation. 

The  development  of  a composite  return  signal 
starts  with  the  simulation  of  the  average  spectrum 
for  the  sea  state  condition.  This  is  the  reverse  of  the 
altimeter  process,  which  develops  the  average  spec- 
trum from  a composite  ocean  return.  For  each  simu- 
lated return,  the  RSS  adds  independent  phase  and 
amiplitude  noise  to  the  spectral  components  to  sim- 
ulate the  random  character  of  the  distribution  of 
ocean  wave  heights.  The  resulting  spectrum  modu- 
lates a single  linear  fm  pulse  to  produce  the  packet 
of  frequency-dispersed  pulses  that  is  received  and 
processed  by  the  altimeter. 

The  average  spectrum  used  to  generate  the  return 
signal  is  derived  from  a mathematical  model  de- 
veloped by  Miller  and  Brown  (Ref.  1).  The  model 
provides  the  normalized  average  power  of  the  return 
signal  versus  time  or  frequency,  when  the  effective 
transmit  pulse  width,  the  rms  sea  state,  the  altim- 
eter antenna  heamwidth,  and  the  nominal  altitude 
are  known. 

Synchronization  of  the  altimeter  with  RSS  timing 
is  accomplished  by  closed-loop  control  of  an  80  MHz 


voltage  controlled  oscillator  in  the  RSS.  The  RSS 
receives  the  altimeter's  linear  fm  transmit  pulse  and 
mixes  it  with  an  internally  generated  1.0  pulse  with 
the  same  characteristics  but  offset  in  frequency.  A 
frequency  discriminator  centered  at  the  offset  fre- 
quency detects  misalignment  between  the  two  pulses 
and  adjusts  the  RSS  oscillator  to  null  the  error. 

The  RSS  is  a major  component  of  a complement 
of  equipment  used  for  testing  and  evaluating  the 
SFASAT  A radar  altimeter.  This  test  system  is  called 
the  Radar  Altimeter  System  Evaluator  (RASE).  It 
provides  complete,  automatic,  computer-based  testing 
and  control  of  the  altimeter.  Options  for  partial  or 
complete  manual  operation  of  the  altimeter  are  also 
provided.  Reference  2 gives  a detailed  description  of 
the  facility. 

The  RSS  operation  can  be  controlled  and  moni- 
tored by  the  RASE  computer  or  by  an  operator  at 
the  front  panel.  Control  consists  of  establishing  the 
parameters  that  determine  the  sea  state  to  be  simu- 
lated, setting  the  altitude,  including  rate  or  accelera- 
tion, and  adjusting  the  return  signal  power  level. 
Once  these  basic  parameters  have  been  entered,  the 
RSS  generates  the  appropriate  return  signal,  properly 
timed,  solely  in  response  to  the  altimeter  transmis- 
sions. 

In  addition  to  the  hardware  that  generates  the  re- 
turn signal,  a complement  of  commercial  test  equip 
merit  is  included  with  the  RSS.  This  provides  a means 
of  verifying  the  RSS  operation  and  measuring  the 
characteristics  of  the  altimeter  transmissions.  Test 
points  for  these  purposes  arc  provided  on  the  front 
panel  of  the  RSS. 
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Also  for  test  purposes,  a portion  of  the  re- 
turn signal  is  diverted  to  the  receiver  section  of  the 
RSS  where  full  deramp  processing  takes  place.  Here 
the  linear  fm  carrier  signal  is  used  for  the  LO. 
Spectrum  analysis  and  averaging  of  the  receiver  if 
output  allow  verification  of  the  desired  characteristics 
of  the  return  signal. 

Figure  3 is  a plot  of  the  average  spectrum  ob- 
served at  the  RSS  receiver  if  output  for  the  ocean 
return  simulation.  Expected  values  determined  from 
the  mathematical  model  of  Ref.  1 are  presented  on 
the  plot  for  comparison. 
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INTRODUCTION 


The  development  of  energy  systems  and  the  assessment  of  the  environ- 
mental impact  of  energy  facilities  have  become  significant  program  areas  for 
APL.  Work  on  these  important  civil  challenges  is  increasing. 

API.  began  studying  the  environmental  impact  of  energy  facilities  in  1971 
in  collaboration  with  The  Johns  Hopkins  University  Department  of  Geog- 
raphy and  Environmental  Engineering  and  the  Chesapeake  Bay  Institute.  Since 
then,  APL  has  had  the  responsibility  to  evaluate  for  the  State  of  Maryland 
the  potential  environmental  and  social  impact  of  all  power-generating  facilities 
proposed  for  development  in  the  State.  This  work  has  led  to  projects  and 
contributions  in  a wide  variety  of  technical  areas  related  to  the  impact  of 
energy  facilities,  impact  mitigation,  and  facility  location. 

One  such  project  involves  the  prediction  of  probable  locations  for  future 
power  plants  within  a multistate  region  under  various  scenarios  for  the 
growth  of  power  demand.  1 he  probable  locations  predicted  by  a computerized 
optimization  model  eventually  will  be  used  in  analyzing  the  environmental 
impact  of  alternative  future  energy  supply  strategies.  The  first  stage  of  this 
work  is  reported  in  the  first  article. 

APL  is  contributing  to  the  development  of  a number  of  new'  energy  sys- 
tems. By  1973,  it  had  analyzed  energy  technologies  and  concluded  that  the 
use  of  low-grade  heat  will  provide  increasingly  important  alternatives  to 
scarce  fossil  fuels.  These  sources  include  solar  thermal,  geothermal,  and  waste 
heat.  Programs  were  initiated  in-house  in  the  areas  of  ocean  thermal  energy 
(solar  energy  stored  in  tropical  oceans),  geothermal  applications,  and  annual 
storage  systems  that  save  summer  heat  for  winter  use. 

Subsequently,  APL  has  become  a leader  in  developing  the  concept  of  utiliz- 
ing  ocean  thermal  gradients  to  produce  ammonia.  Domestic  ammonia,  used 
in  nitrogen  fertilizers,  is  now  made  from  natural  gas.  The  APL  version  of 
Oiean  Thermal  Energy  Conversion  (OTEC)  could  have  a substantial  impact 
on  natural  gas  consumption  in  the  United  States.  Efficient  heat  exchange  at 
small  temperature  differences  is  important  to  the  success  of  the  OTEC  con- 
cept, and  APL  s efforts  have  included  heat  exchanger  design  and  evaluation. 
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Biofouling  could  degrade  the  performance  of  heat  exchangers.  In  the  second 
article,  experiments  to  measure  the  rate  of  biofouling  and  its  effect  on  heat 
transfer  are  reported.  Preliminary  studies  of  cleaning  mechanism  are  discussed. 
Work  is  also  proceeding  on  issues  relating  to  the  commercialization  of  OTEC. 
Studies  are  being  initiated  on  the  design  and  response  to  stress  of  the  massive 
concrete  pipe  that  would  be  required  to  bring  cooling  water  from  the  deep 
(Kean  (3000  ft)  to  the  surface. 


The  work  reported  in  the  third  article  on  the  Community  Annual  Storage 
Energy  System  (CASES)  seeks  to  explain  how  summer  heat  can  be  eco- 
nomically recycled  for  use  in  winter  by  means  of  annual  storage  on  a commu- 
nity basis.  It  is  anticipated  that  subsequent  phases  of  the  effort  will  address 
site-specific  questions  and  lead  to  a demonstration  and  evaluation  program. 

A number  of  other  energy-related  programs  are  continuing.  Pioneering 
work  has  been  done  in  the  development  of  efficient  flyu  heels  from  filjmentar y 
materials.  The  next  phase  of  this  work  will  be  to  demonstrate  a complete 
low-cost  flywheel  energy-storage  system  based  on  110  V,  60  Hz  power. 


APL  is  responsible  for  studies  of  geothermal  resources  in  the  Department 
of  Energy/Di  vision  of  Geothermal  Energy  Region  5.  Scenarios  are  being 
developed  that  outline  the  potential  magnitude  of  identified  resources  with 
development  and  use  schedules,  and  problems  that  must  be  surmounted  before 
the  resource  could  be  used  are  being  identified.  Site-specific  studies  of  poten- 
tial geothermal  energy  applications  are  being  initiated  at  three  locations. 

In  addition,  institutional  impediments  to  the  development  of  low-head 
hydroelectric  power  are  being  studied.  Research  is  being  conducted  on  the 
systematic  development  of  organic  conductors  and  on  the  active  acoustic  detec- 
tion of  leaks  in  underground  natural  gas  lines. 


The  articles  included  herein  are  only  a sample  of  the  energy  work  being 
done  at  APL.  Many  areas  indicated  in  this  introduction  are  not  reported. 
These  and  other  projects  will  be  reported  in  subsequent  Accomplishments 
volumes. 


REGIONAL  MODEL  FOR  LOCATING 
POWER  PLANTS 
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Assessing  the  environmental  impact  of  various 
courses  of  action  to  meet  the  Nation's  needs  for  elec- 
tric energy  is  impossible  without  some  idea  of  the 
locations  of  future  power  plants.  APL  has  developed  a 
computer  model  that  predicts  probable  locations  for 
large  numbers  of  future  power  plants  within  a multi- 
state region. 


BACKGROUND 

In  the  past  few  years,  "energy  scenarios"  have 
gained  increasing  favor  among  planners  for  the  use 
of  electric  power.  An  energy  scenario  is  a hypothetical 
course  of  events  that  might  occur  as  the  result  of  a 
given  set  of  energy  supply  decisions,  usually  on  a 
national  level.  Recently,  planners  have  come  to  realize 
that  projections  of  energy  variables  on  a national 
level  are  of  little  help  in  estimating  the  chances  of 
achieving  a given  energy  supply/demand  state  and 
in  assessing  the  possible  ways  to  achieve  it.  Analysis 
at  a finer  geographic  level  is  required.  The  concept 
of  a multistate  region  is  in  favor  because  physical 
resources  (e.g.,  river  basins)  often  cover  several 
states  and  because  electric  power  production  and  dis- 
tribution systems  often  extend  over  several  states. 

Current  computer  models  produce  schedules  of 
power  plants  needed  to  implement  a given  electric 
energy  scenario,  on  both  a national  and  a regional 
basis.  The  regional  versions  simply  provide  a list  of 
specified  power  plants  that  must  be  brought  into 
service  in  a multistate  region  within  a specified  time. 
In  order  to  deal  with  environmental  impact  questions, 
however,  a finer  geographic  resolution  is  required. 
APL’s  computer  model  (Ref.  1)  provides  this  finer 
spatial  resolution  by  predicting  probable  locations 
within  a multistate  region  for  large  numbers  of 
power  plants  associated  with  the  given  energy  sup- 
ply scenarios.  The  locations  are  not  precise  enough 
for  actual  siting  efforts;  however,  they  are  sufficiently 
accurate  for  use  in  studies  that  attempt  to  evaluate 
the  cumulative  environmental  impact  of  the  scenarios. 


DISCUSSION 

The  model  has  a grid  of  cells,  each  240  sq  mi, 
which  corresponds  exactly  to  the  15-minute  quad- 
rangle map  system  of  the  United  States  Geological 
Survey.  The  cell  is  the  smallest  element  of  locational 


resolution.  After  manual  screening  to  eliminate  cells 
with  obvious  negative  features  (such  as  lack  of  ade- 
quate cooling  water),  the  remaining  cells  are  desig- 
nated as  potential  power  plant  locations.  Figure  1 a 
shows  the  cell  grid  and  potential  location  cells  super- 
imposed on  a map  of  the  mid-Atlantic  test  region. 

The  regional  electric  power  demand  for  some  fu- 
ture time  (the  year  2000  was  used  in  initial  tests)  is 
input  to  the  model,  which  allocates  the  demand  to 
load  centers  across  the  region  on  the  basis  of  popu- 
lation projections.  The  model  then  uses  linear  pro- 
gramming to  locate  nuclear  and  coal-fired  power 
plants  to  satisfy  the  demand  while  optimizing  a speci- 
fied objective  and  while  ensuring  that  a set  of  con- 
straints is  not  violated.  The  objectives  and  constraints 
are  intended  to  represent  real-world  goals  and  re- 
strictions associated  with  power  plant  siting;  a great 
deal  of  input  from  regional  electric  utility  planners 
and  state  regulatory  officials  was  needed  in  their  for- 
mulation. 

Constraints  currently  embodied  in  the  model  in- 
clude limitations  on  consumption  and  withdrawal  of 
cooling  water,  restrictions  on  the  proximity  of  pop- 
ulations to  nuclear  plants,  limits  on  coal  availability, 
and  standards  for  ambient  air  quality.  The  cooling 
water  constraints  are  of  special  interest,  since  they 
are  formulated  for  the  entire  region,  which  encom- 
passes all  or  part  of  six  major  river  basins.  The  model 
shows  the  effect  on  all  downstream  users  of  all  con- 
sumptive withdrawals  at  any  point  in  a basin  and 
thus  could  be  used  as  a planning  tool  for  water 
resources. 

The  objectives  chosen  for  the  initial  development 
of  the  model  include  the  minimization  of  required 
electric  power  transmission,  of  required  coal  trans- 
portation, of  required  reservoir  capacity  (region-wide) 
necessary  to  provide  low-flow  augmentation,  and  of 
aggregate  population  proximity  to  nuclear  plants.  The 
first  three  objectives  serve  as  surrogates  for  costs  that 
depend  on  location,  while  the  population  objective 
reflects  the  component  of  regulatory  and  public 
opinion  that  favors  siting  away  from  population  cen- 
ters. The  objectives  can  be  used  in  weighted  linear 
combinations,  the  weights  being  proportional  to  the 
supposed  importance  of  the  corresponding  objectives. 
Interaction  with  State  and  Federal  officials,  industry 
planners,  and  the  public  is  necessary  to  determine  the 
relative  importance  of  each  objective. 
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Note: 

These  locations  are  not 
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locations  in  an  energy 
scenario  modeling 
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Fig.  I Mid-Atlantic  test  region,  (a)  Eligible  cells  after  screening;  (b) 
location  pattern  that  minimizes  coal  transportation;  (c)  location  pattern 
that  most  closely  resembles  utility  siting  plans. 


For  initial  testing,  the  model  was  exercised  by 
requiring  it  to  locate  25  000  MWe  of  nuclear  capac- 
ity and  25  000  MWe  of  coal-fired  capacity,  an 
amount  corresponding  to  a 2.7%  annual  regional 
growth  rate  from  1977  to  2000.  Figure  lb  shows  the 
resultant  locations  for  an  extreme  case  in  which  the 
minimization  of  coal  transportation  was  given  the 
greatest  weight.  The  location  pattern  that  most  closely 
resembles  utility  siting  plans  is  shown  in  Fig.  lc.  This 


pattern  arose  when  minimization  of  transmission  and 
of  flow  augmentation  were  given  equal  and  primary 
importance,  coal  transportation  was  made  much  less 
important,  and  population  proximity  was  not  con- 
sidered at  all  (beyond  meeting  legal  regulations). 


Future  efforts  will  be  directed  toward  determining 
definitive  objective  weights  and  extending  the  model 
to  other  regions  of  the  United  States. 
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REFERENCE 

1.  Probable  Pninbutiou  of  F.fiutnt  Sources  from  Energy  Supply 
und  Comer \ton  APL/JHU  Final  Report  (to  be  published)- 


BIOFOULING  AND  CLEANING 
EXPERIMENTS  FOR  OTEC  HEAT 
EXCHANGERS 


In  the  APL  Ocean  Thermal  Energy  Conversion 
(OTEC)  plant-ship  concept,  large-diameter  (3-in.) 
aluminum  tubes  with  ammonia  flowing  inside  are  being 
proposed  for  use  as  low-cost  multipass  heat  ex- 
changers. Experiments  to  measure  the  rate  of  biofoul- 
ing and  its  effect  on  exterior  heat  transfer  have  been 
conducted  by  APL  at  the  Natural  Energy  Laboratory  of 
Hawaii.  Ke-Ahole  Point,  Hawaii.  The  results  are  directly 
applicable  to  the  design  of  the  APL  heat  exchanger. 


BACKGROUND 

In  an  OTEC  plant,  warm  surface  seawater  will  be 
used  in  evaporators  to  vaporize  ammonia.  The  am- 
monia vapor  will  drive  a power  turbine  and  then 
will  be  condensed  by  means  of  heat  exchange  with 
cold  seawater  drawn  from  a depth  of  2500  ft.  The 
electric  power  produced  by  the  turbine  will  be  used 
on  board  to  produce  ammonia,  aluminum,  or  other 
energy-intensive  products. 

The  plants  will  succeed  commercially  only  if  the 
overall  system  cost  is  low.  A promising  low-cost  heat 
exchanger  concept  (Ref.  1)  has  arrays  of  large-diam- 
eter (3  in.)  multipass  aluminum  tubes  for  both  the 
evaporator  and  the  condenser,  with  ammonia  as  the 
working  fluid  flowing  inside  the  tubes.  Seawater  is 
pumped  to  head  ponds  above  the  tubes  and  flows 
vertically  downward  over  them  by  gravity. 

The  use  of  large-diameter  tubes  results  in  a con- 
figuration that  has  sufficient  space  in  the  vertical 
planes  between  rows  of  tubes  to  permit  a low-cost 
system  of  cleaning  heads  to  remove  the  biofouling 
from  the  outside  of  the  tubes. 

Since  biofouling  is  a critical  factor  that  affects 
the  long-term  performance  of  the  heat  exchanger, 
experiments  were  needed  to  measure  its  effect  and  to 
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develop  low -cost  cleaning  methods  in  an  environment 
similar  to  that  of  actual  OTEC  plant-ship  operation. 
An  apparatus  was  designed,  built,  and  deployed  in 
Hawaii,  a potential  OTEC  site,  and  experiments  have 
been  conducted  since  mid-1977. 


DISCUSSION 

Figure  1 is  a schematic  of  the  biofouling  test  ap- 
paratus. It  simulates  an  array  of  heat  exchanger  tubes 
arranged  in  one  of  the  geometries  under  considera- 
tion. The  array  is  located  in  a box  through  which 
seawater  is  passed.  There  are  five  lateral  rows  of 
tubes  and  10  tubes  in  each  longitudinal  column,  with 
3 separation  of  1*4  in.  between  adjacent  vertical  tube 
hanks. 

The  first,  fourth,  seventh,  and  tenth  tubes  (tubes 
A,  B,  C and  D,  respectively)  of  the  center  row  are 
instrumented  as  described  below  to  measure  heat 
transfer.  The  remaining  tubes  are  made  of  Alclad 
3004  aluminum  and  are  vented  to  overcome  buoy- 
ancy. These  dummy  tubes  are  held  in  position  by 
0.75-in. -thick  spacers  epoxied  in  place  at  the  ends. 

Each  instrumented  tube  (Fig.  2)  consists  of  a 
6061  aluminum  heat  sink  with  an  exposed  length  of 
5I/2  in.  One  end  of  the  tube  extends  through  one 
side  of  the  test  box.  The  heat  sink  is  heated  by  a 
5t/,-in.-long,  400- W (nominal)  quartz  lamp  that  is 
inserted  into  a 1 V^-in. -diameter  axial  hole.  The  heat 
sink  is  divided  into  four  quadrants  and  instrumented 
with  thermocouples  so  that  heat  transfer  is  measured 
on  the  top,  bottom,  and  both  sides.  After  the  lamp 
is  turned  on,  the  heat  sink  reaches  a steady-state 
temperature  on  each  quadrant  in  approximately  five 
minutes.  The  heat  transfer  coefficients  are  computed 


from  the  magntiude  of  these  temperatures,  and  the 
week-to-week  changes  in  the  coefficients  determine 
the  biofouling  coefficients.  (The  biofouiing  resistance 
coefficient  is  defined  as  R,  = l/k,oulrd  — \/hcUan, 
where  h is  the  measured  heat  transfer  coefficient  for 
the  surface.) 


All  dimensions  are  in  inches  (drawing  is  not  to  scale) 

• Instrumented  tube  (4  in.  dia.) 

O Dummy  tube  (4  in.  dia.) 


I* 78 1 


— ' — 3 1 


58 


Fig.  I Test  setup  showing  the  staggered  tube  arrange- 
ment simulating  the  full-scale  geometry  of  the  OTEC 
plant-ship  heat  exchanger. 


As  of  September  1977,  two  series  of  tests  have 
been  run.  During  series  I (}0  May  to  15  July  1977), 
the  flow  rate  through  the  assembly  was  set  at  approxi- 
mately 1775  gal/min,  which  gives  a maximum  flow 
velocity  of  3 ft/s  in  the  minimum  gap  of  the  tube 
array.  The  top  of  the  test  tank  was  left  uncovered, 
and  data  were  taken  at  approximately  weekly  inter- 
vals. Series  II  (10  August  to  21  September  1977) 
duplicated  the  conditions  of  the  first  except  that  the 
top  of  the  test  tank  was  covered  to  exclude  light  from 
the  upper  parts  of  the  test  assembly. 

Results  from  these  tests  are  illustrated  in  Fig.  3, 
in  which  the  average  fouling  coefficient  for  each  tube 
is  plotted  against  time;  details  of  the  fouling  growth 
on  each  quadrant  may  be  found  in  Ref.  2.  An  APL- 
designed  scrubber  mechanism  using  Scotchbrite  pads 
that  fitted  into  the  spaces  of  the  tube  bank  was  used 
at  the  end  of  Series  II.  The  following  observations 
were  made: 

1.  Bifouling  accumulates  at  a rate  to  produce  in 
six  weeks  a resistance  to  heat  transfer  meas- 
ured as  0.00045  Btu-lft5-h-°F  with  a standard 
deviation  of  0.00011. 

2.  Light  versus  dark  environment  is  indicated 
to  have  a significant  effect  on  the  rate  of  bio- 
fouling only  in  areas  of  low  water  velocity. 

3.  Tube  sides  tend  to  develop  a smaller  fouling 
coefficient  than  the  tube  top  or  bottom.  The 
tube  bottom  tends  to  have  the  highest  fouling 
coefficient,  especially  if  there  is  no  other  tube 
directly  below  it. 

4.  The  biofouling  accumulation  of  six  weeks  is 
readily  removed  by  light  scrubbing  with  a plas- 
tic sponge  (Scotchbrite). 


Solenoid  latch 


Fig.  2 Schematic  of  an  instrumented  tube. 


89 


Fig.  3 Average  fouling  coefficients  of  tubes  (four 
quadrants)  plotted  versus  time  for  two  test  series, 
based  on  available  data. 

5.  The  cleaning  mechanism  was  effective  in  clean- 
ing the  side  quadrants  of  the  array  but  was 
only  partially  effective  on  the  top  and  bottom 
quadrants.  The  results  are  encouraging  but 
indicate  the  need  for  improvement  of  the 
device. 


6.  The  presence  of  iron  oxide  from  the  pumps 
may  have  affected  the  character  and  rate  of 
growth  of  the  biofouling  layer. 

Subsequent  testing  is  planned  with  the  apparatus 
in  order  to  develop  a cleaning  system  directly  appli- 
cable to  the  APL/OTEC  concept.  Since  several  cycles 
of  fouling  and  cleaning  are  necessary  to  determine 
the  effectiveness  for  design  purposes,  the  first  tests 
give  a preliminary  indication  that  biofouling  can  be 
controlled  cost-effectively. 
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SIMULATION  OF  A COMMUNITY  ANNUAL 
STORAGE  ENERGY  SYSTEM 


The  Community  Annual  Storage  Energy  System 
(CASES)  is  a thermal  utility  that  provides  heating  and 
cooling  services  to  a community  via  water  pipelines 
from  a central  energy-storage  facility.  A simulation  by 
computer  modeling  has  been  developed  for  calculating 
the  requirements  of  an  entire  community  when  heating 
and  cooling  is  supplied  by  means  of  CASES.  It  calcu- 
lates the  size  and  cost  of  all  water-source  heat  pumps 
and  tan  coil  cooling  units  in  each  building  of  the  com- 
munity together  with  the  electrical  energy  they  require 
emcti  hour  of  the  year  to  cope  with  varying  weather 
ondrt’one  It  also  calculates  the  size  and  cost  of  every 
< r»  /serf  to  distribute  heating  and  cooling  water  after 
erst  determining  the  peet  flow  condition  in  each  pipe. 


BACKGROUND 

The  key  element  in  CASES  is  the  storage  of  sur- 
plus thermal  energy’  in  a central  facility.  Much  of  the 
heat  required  by  the  community  is  collected  ''energy- 
free"  at  times  of  heat  surplus  and  is  stored  until 
needed.  I e is  also  collected  energy-free  in  winter 
and  saved  to  he  used  for  summer  coding  Water- 
source  heat  pumps  supply  heating  and  cooling  when- 
ever the  store  of  energy-free  ice  and  thermal  water 
is  inadequate  for  the  community's  needs. 


Since  much  of  the  energy  required  for  heating  and 
cooling  is  collected  as  surplus  summer  heat  or  winter 
ice,  energy  from  other  sources  for  heating  and  cool- 
ing can  he  reduced  by  as  much  as  80Cf.  It  appears 
possible  that  the  total  capital  invested  for  heating 
and  cooling  may  be  less  than  electric  utility  compa- 
nies, building  owners,  and  oil  companies  would  have 
to  invest  for  electrical  air-conditioning  and  conven- 
tional heating  systems.  Building  owners  can  sub- 
scribe to  the  CASES  utility  instead  of  purchasing 
independent  heating  and  cooling  equipment,  thereby 
considerably  reducing  building  costs. 


SUMMARY 

The  CASES  simulation  begins  by  analyzing  the 
heating  and/or  cooling  required  each  hour  for  each 
distinct  building  type  found  in  the  community  (the 
HCLOAD  program  module).  The  annual  peak  heat- 
ing load  is  used  to  size  the  water-source  heat  pump 
in  each  building;  the  annual  peak  cooling  load  simi- 
larly determines  the  size  of  the  fan-coil  cooling 
units.  Next,  the  demand  (Ib/h)  for  distribution- 
pipeline  water  at  each  building  site  corresponding  to 
the  heating  and  cooling  load  (Btu/h)  schedules  is 
calculated,  making  use  of  the  performance  charac- 
teristics of  the  selected  heat  pumps  and  fan-coil  units. 

The  electrical  energy  required  each  hour  to  operate 
this  heating  and  cooling  equipment  in  each  building 
is  also  calculated  (the  CDIST  program  module). 
Finally  the  size  of  every  pipe  used  in  the  distribu- 
tion system  is  determined  so  that  the  peak  annual 
flow  will  be  less  than  a specified  maximum  flow  rate. 
Thus,  by  means  of  the  HCLOAD  and  CDIST  pro- 
gram modules,  the  entire  community  is  reduced  to 
an  hourly  schedule  of  water  and  electric  power  re- 
quirements, and  the  cost  of  all  community  compo- 
nents (heat  pumps,  cooling  coils,  and  pipelines)  is 
known. 


DISCUSSION 

Cooling.  Water  ai  40°  F is  distributed  through 
pipe  that  extend  throughout  the  community  and  is 
then  passed  through  cooling  coils  in  the  various  build- 
ings. Leaving  these  cooling  coils  at  a temperature  of 
about  60°F,  the  water  is  returned  to  a central  ice 
deposit  where  it  is  again  cooled  to  40°  F prior  to  being 
repumped  through  the  flow  loop.  The  heat  imparted 
to  the  water  when  it  is  used  for  cooling  is  saved 


within  the  system  until  needed  or  until  it  can  be 
rejected  without  an  energy-consuming  heat  pump. 
This  is  in  contrast  to  conventional  air-conditioning 
systems  that  remove  heat  from  buildings  and  eject  it 
by  means  of  a heat  pump  (air  conditioner)  into  the 
atmosphere,  where  it  is  lost. 

Heating.  The  community  is  heated  by  water-source 
heat  pumps.  Conventional  air-source  heat  pumps  be- 
come inefficient  when  the  outside  air  temperature 
falls  below  32°F,  but  the  heat  pumps  used  in  CASES 
remain  highly  efficient  because  60°1:  water  distributed 
from  the  central  CASES  plant  is  always  available 
as  a source  of  heat.  The  heat  pumps  extract  ‘hermal 
energy  from  the  water,  lowering  its  temperature  to 
40°  F.  The  40°  F water  is  returned  to  the  central 
facility  in  the  same  40°F  mainline  used  for  cooling 
water,  thus  avoiding  additional  piping.  Likewise,  the 
pipeline  used  to  supply  60°F  heat-source  water  to  the 
heat  pumps  is  the  same  60°F  pipeline  used  to  return 
water  to  the  central  ice  store  during  the  cooling 
cycle.  Thus,  CASES  has  only  two  uninsulated  pipe- 
lines, one  for  40°F  water  and  one  for  60°F  water, 
yet  both  heating  and  cooling  are  always  available  to 
all  users. 

Part  of  the  heat  that  is  used  to  warm  the  40°  F 
water  to  60°  F for  recirculation  to  the  community 
comes  from  the  recovery  of  excess  heat  as  discussed 
above.  The  remainder  comes  from  the  central  plant 
and  will  be  simulated  in  future  CASES  work.  The 
central  plant  feels  only  the  net  demand  for  pipeline 
water.  A significant  fraction  of  the  heating  and  cool- 
ing demand  is  satisfied  by  the  synergistic  interaction 
among  the  diverse  buildings  of  the  community. 

Community.  It  is  necessary  to  specify  the  com- 
munity in  order  to  evaluate  quantitatively  the  extent 
of  synergistic  interaction  among  the  various  buildings 
of  the  community  and  to  predict  the  cost  of  heat 
pumps,  cooling  coils,  and  pipelines.  The  selected 
model  community  is  similar  to  Wilde  Lake  Village 
in  Columbia,  Maryland.  It  has  8000  residents  (2300 
households)  and  866  buildings  in  three  neighbor- 
hoods (NE,  NW,  and  N),  which  share  a common 
village  center,  all  on  2500  acres  of  land.  There  is  a 
job  in  the  community  for  each  worker.  The  physical 
layout  of  the  community  is  shown  in  Fig.  1.  Table  1 
gives  the  demography  and  distribution  of  residential 
units.  Each  neighborhood  has  its  own  elementary 
school  and  civic-commercial  center.  The  middle 
school,  high  school,  office  buildings,  factory,  village 
shopping  center,  and  CASES  facility  are  centrally 
located,  as  shown  in  Fig.  1.  A detailed  study  of  the 
demographics  of  the  community  was  conducted  so 
that  the  relative  number  of  employment  and  residen- 
tial facilities  could  be  determined. 
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Fig.  1 Pipeline  layout  route  for  hypothetical  CASES  community. 


RESULTS 

There  were  interesting  results  from  the  simulation 
of  only  the  community  and  the  distribution  system. 
An  866-building  community  almost  identical  to  that 
described  in  Fig.  I required  3.2  X 10"  Btu  of 
heating  and  8.8  X 1010  Btu  of  cooling  in  Wash- 
ington area  weather  in  1967.  Almost  half  of  the 
cooling  demand  could  be  met  by  "waste-cold  re- 
covery" (chilled  water  produced  by  heat  pumps  used 
to  heat  other  buildings),  and  about  one-eighth  of  all 
heating  was  met  by  waste-heat  recovery.  1 he  heat- 
ing and  cooling  equipment  required  6.3  X 10'"  Btu 
of  electrical  energy.  Thus,  the  combined  seasonal 
coefficient  of  performance  was  6.5,  not  counting  cen- 
tral plant  effects.  The  total  CASES  capital  investment 
required  for  this  community  was  17.4  million  (89 c/c 
in  heat  pumps,  4%  in  cooling  coils,  and  7%  in  the 
distribution  system).  This  is  less  than  half  the  capital 
required  to  expand  a winter- peak  electric  power  plant 


to  provide  electrical  heat  to  the  same  community  (a 
plant  cost  of  $700/kWe  was  assumed). 

No  effort  has  been  made  to  optimize  these  CASES 
results,  and  it  must  be  emphasized  that  the  central 
CASES  facility  is  not  included  in  the  simulation. 
However,  both  the  cost  and  the  losses  of  electrical 
energy  during  transmission  have  also  been  neglected 
in  the  analysis  of  the  electric  heating  alternative. 

CASES  optimization  will  proceed  when  the  central 
plant  simulation  (CAPS  program  module)  has  been 
completed. 
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FUNDAMENTAL  RESEARCH 
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INTRODUCTION 


APL's  continuing  vitality  as  an  institution  stems  in  part  from  its  dual  role 
as  a technology  center  for  solving  national  problems  and  as  a division  of  the 
University.  This  duality  has  provided  a perspective  on  defense  and  civil  con- 
cerns of  national  importance  that  looks  to  the  fundamental  issues  and  trends 
underlying  the  immediate  pressing  problems. 

Institutionally,  APL’s  university  origins  are  reflected  in  its  continuing  sup- 
port for  research  in  the  basic  sciences  of  physics,  chemistry,  biology,  and 
mathematics.  Much  of  the  work  is  relevant  to  the  technologies  used  in  Labora- 
tory programs  concerned  with  national  defense,  energy,  and  biomedicine. 
However,  the  motivation  reflects  a desire  to  understand  the  basic  physical  prin- 
ciples involved.  This  attitude  provides  continuing  links  with  the  academic 
and  medical  divisions  of  Johns  Hopkins  as  well  as  with  the  international  sci- 
ence community.  At  the  same  time,  APL's  competencies  in  these  areas  arc 
a reservoir  for  technological  innovation  and  a resource  for  solving  specific 
problems  where  the  application  of  current  technology  or  engineering  practice 
is  inadequate.  Research  interests  include  topics  in  materials  sciences,  biology, 
chemistry,  physics,  and  mathematics.  Several  notable  contributions  to  these 
fields  are  summarized  in  the  articles  that  follow. 

In  materials  science,  the  work  on  corrosion  processes  in  aluminum  con- 
stitutes an  attempt  to  understand  the  basic  electrochemical  and  surface  reaction 
mechanisms  responsible  for  one  important  type  of  corrosion.  Such  an  under- 
standing should  provide  a guide  for  a more  general  understanding  of  this 
technologically  important  topic.  In  fluid  flow,  a closed-form  solution  is  pre- 
sented for  a problem  in  the  detachment  of  flow  at  a boundary.  This  work  has 
direct  application  in  aerodynamics  and  in  arterial  flow.  In  the  area  of  applied 


mathematics,  large  classes  of  problems  can  be  solved  only  approximately.  The 
next  article  is  an  example  of  new  developments  in  the  construction  of  general 
methods  to  obtain  approximate  solutions.  Methods  now  under  study  at  APL 
could  find  wide  application  in  many  areas  of  science  and  engineering. 

Chemical  physics  is  represented  by  two  articles.  The  work  on  chemically 
induced  polarization  effects  describes  the  development  of  a general  analytical 
representation  of  reactions  that  produce  electron  and  nuclear  spin  polariza- 
tions. These  processes  are  important  conceptually  for  understanding  the  mech- 
anisms of  an  unusual  class  of  reactions.  For  example,  nuclear  spin  polarized 
targets  are  now  used  to  enhance  the  yield  of  certain  reactions  in  accelerator 
physics.  They  could  also  be  important  technologically.  The  next  article  sheds 
new  light  on  energy  levels  and  bonding  in  XeCI.  These  molecules  arouse 
considerable  theoretical  interest  and  also  are  technologically  important,  since 
they  constitute  the  active  medium  in  the  rare  gas  halide  excimer  laser  now  un- 
der development  as  an  ultraviolet  source. 

Finally,  the  article  on  intensity  correlation  spectroscopy  provides  a view  of 
work  now  in  progress  in  the  area  of  biological  physics.  It  describes  the  de- 
velopment of  analytical  tools  needed  to  apply  optical  correlation  spectroscopy 
to  the  measurement  of  the  properties  of  fluid  suspensions. 

These  articles  represent  larger  programs  that  include  problems  in  hydro- 
dynamics, energy  conversion,  reaction  processes  at  solid  surfaces  (including 
corrosion),  and  biomedicine.  Other  major  problem  areas  concern  the  elec- 
tronic and  optical  properties  of  technologically  important  organic  and  metal- 
lic organic  compounds  and  the  development  of  high-energy  laser  devices. 


LOCALIZED  CORROSION  OF  ALUMINUM: 
BLISTER  FORMATION  AS  A PRECURSOR 
OF  PITTING 


Various  localized  corrosion  processes  are  thought  to 
begin  by  similar  mechanisms.  A previously  unk own 
step  in  the  localized  corrosion  of  aluminum  known  as 
pitting  has  been  discovered  in  our  laboratory,  the 
passivating  oxide  on  aluminum  has  been  observed  to 
blister  and  break  prior  to  the  onset  of  pitting  corro- 
sion. 


BACKGROUND 

The  cost  of  corrosion  in  the  United  States  has  been 
estimated  to  be  $20  to  30  billion.  As  much  as  90% 
of  this  loss  has  been  associated  with  localized  corro- 
sion processes  such  as  stress-corrosion  cracking,  pit- 
ting, fatigue,  and  crevice  corrosion.  These  processes 
are  thought  to  have  similar  initiation  mechanisms. 
That  is,  the  first  events  are  probably  the  same,  even 
though  the  eventual  gross  manifestation  of  each  cor- 
rosion type  has  distinct  characteristics  (for  many 
articles  on  localized  corrosion  and  initiation  mecha- 
nisms, see  Ref.  t).  Our  research  has  dealt  with  pit- 
ting corrosion,  which  takes  the  form  of  pits  in  the 
surface  of  the  metal. 

Associated  with  pitting  is  an  electrical  potential 
called  the  pitting  potential.  If  the  potential  of  a 
metal  in  a corrosive  environment,  such  as  aqueous 
chloride  solution,  rises  above  the  pitting  potential, 
pitting  corrosion  begins.  Blisters  in  the  oxide  film 
on  the  metal  have  been  observed  in  our  laboratory 
before  the  onset  of  gross  pitting.  Our  observations 
indicate  that  such  blisters  can  be  initiated  by  short- 
duration  (less  than  1 ms),  low-voltage  pulses  The 
blisters  grow  and  eventually  rupture,  at  which  time 
pitting  corrosion  begins.  Thus,  blisters  are  experi- 
mentally demonstrable  precursors  of  pits  in  this  in- 
stance (Ref.  2). 


The  experimental  configuration  was  similar  to  that 
shown  schematically  in  Fig.  1.  A light  microscope 
was  focused  on  the  electrode  to  observe  any  visible 
phenomena.  The  potential  on  the  aluminum  electrode 
with  respect  to  the  reference  electrode  was  main- 
tained above  aluminum's  pitting  potential.  A very 
short  potential  pulse  was  applied  to  the  aluminum 
electrode.  This  pulse  initiated  the  growth  of  blisters. 


Auxiliary 

electrode 

(stainless 

steel 

beaker) 


Fig.  1 The  experimental  configuration. 


DISCUSSION 

The  study  of  localized  corrosion  in  aluminum  re- 
quires the  investigator  to  characterize  the  oxide  film 
that  covers  and  passivates  the  metal.  This  oxide  can 
be  grown  in  air  as  a nonporous  film  2 to  4 nm  thick. 
Films  of  greater  thickness  can  be  grown  in  solution 
by  a process  known  as  anodization.  The  films  can  be 
of  various  porosities,  depending  on  the  nirrent  and 
the  chemicals.  For  our  investigation,  nonporous  films 
20  to  200  nm  thick  were  grown. 


A typical  unbroken  blister  80  gin  in  diameter  is 
shown  in  Fig.  2.  The  thickness  of  the  specimen's 
oxide  layer  was  80  nm.  Blisters  usually  have  a cir- 
cular shape,  although  their  morphology  can  be  in- 
fluenced by  grain  boundaries  and  other  topographical 
features.  They  appear  under  the  low-power  light 
microscope  a few  seconds  after  the  voltage  pulse  is 
applied.  Blister  growth  typically  is  somewhat  jerky; 
that  is,  the  blister  boundary  does  not  move  steadily 
outward  but  rather  increases  by  irregular  increments 
with  pauses  in  the  boundary  motion.  During  growth 
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Fig.  2 A scanning  electron  micrograph  of  an  unbroken 
blister.  The  beam  voltage  was  20  keV  and  the  stage  tilt 
was  30°. 


there  is  an  observable  frontal  region  extending  out- 
side the  blister  periphery  in  which  there  appears  to 
be  activity  or  stress  of  some  type.  In  the  scanning 
electron  microscope  (SF.M),  the  circumference  of  the 
blister  frequently  shows  signs  of  stress.  The  shape, 
manner  of  growth,  and  stressed  regions  are  all  indic- 
ative of  pressurized  gas  under  the  blistering  oxide. 
The  gas  apparently  is  formed  by  a chemical  reaction 
that  attacks  laterally  along  the  oxide/metal  interface 
as  opposed  to  penetrating  into  the  metal,  although 
the  dark  spots  seen  in  the  blister  are  etch  marks 
caused  by  a minor  penetrating  attack. 

The  micrograph  shown  in  Fig.  2 was  taken  at  20 
keV  beam  voltage.  Micrographs  taken  at  2.3  keV  do 
not  reveal  the  etch  marks;  therefore,  the  marks  are  in 
the  metal  under  the  oxide.  There  are  much  larger 
amounts  of  chloride  in  the  light  contrast  annulus 
about  the  blister  periphery  than  inside  the  blister 
where  chloride  is  occasionally  detectable  or  outside 
the  blister  where  it  is  usually  not  detectable  in  our 
apparatus.  Close  examination  of  the  oxide  over  the 
blister  in  Fig.  2 revealed  no  detectable  cracks  or 
holes.  The  SEM  has  a resolution  of  7 nm,  but  the 
nonconducting  oxide  reduces  this  resolution  some- 
what. 


Several  conditions  must  be  met  if  one  wishes  to 
observe  growing  blisters  initiated  by  a single  pulse 
of  1-ms  duration  or  less.  A smooth  surface  is  favor- 
able. On  rough  surfaces,  the  oxide  frequently  breaks 
soon  after  pulsing.  The  anodized  layer  thickness  also 
plays  a role,  as  might  be  expected.  The  thinner 
layers  break  more  easily;  thus,  blisters  usually  break 
at  smaller  diameters.  The  aluminum  electrode  poten- 
tial must  be  above  the  pitting  potential.  Also,  the 
pulse  amplitude  must  be  large  enough.  The  ampli- 
tude depends  on  oxide  thickness  and  is  somewhat 
variable  from  specimen  to  specimen.  In  general,  2- 
to  3-V  amplitudes  are  required  for  20-nm  oxides; 
4-  to  3-V amplitudes  for  80-nm  oxides;  and  greater 
than  10-V  amplitudes  for  200-nm  oxides.  If  ampli- 
tudes much  greater  than  these  are  used,  the  oxide 
usually  breaks  immediately.  Pulse  duration  was  not 
varied  much  in  our  tests;  however,  if  the  potentiostat 
is  set  at  2 V (referred  to  a saturated  calomel  elec- 
trode) for  a period  of  a few  seconds,  many  blisters 
form  and  grow  under  80-nm  films.  Some  of  the  blis- 
ters break  and  become  pits. 

Both  intragranular  and  grain  boundary  sites  have 
been  observed  as  points  of  blister  initiation  on  the 
aluminum  surface. 

Although  blisters  can  grow  to  sizes  on  the  order 
of  100  g.m  on  smooth  surfaces,  the)’  eventually  rup- 
ture. Upon  breaking,  there  is  a considerable  increase 
in  the  current  (see  Fig.  3).  The  time  at  which  the 
specimen  was  pulsed  is  also  indicated  in  the  figure. 
In  the  case  of  a breaking  blister,  the  current  increases 


Fig.  3 A plot  of  current  versus  time  showing  the  time 
of  pulsation,  the  growth  period,  and  the  time  of  rup- 
ture. 


because  there  is  a greater  area  of  unpassivated  alumi- 
num exposed  to  the  electrolyte. 

A question  of  considerable  significance  is  that  of 
the  transport  of  materials  across  the  oxide  to  the 
aluminum.  No  cracks  or  holes  could  be  found  in 
blistered  oxide  that  had  not  broken  catastrophically. 
A possibility  is  that  such  defects  are  below  the  reso- 
lution of  our  SEM.  The  oxide  over  a growing  blister 
supports  pressurized  gas  that  would  probably  favor 
crack  propagation.  This  does  not  usually  happen  im- 
mediately; when  it  finally  does  occur,  the  cracks  often 
are  outside  the  perimeter  within  which  the  blister 
was  confined  most  of  its  lifetime. 

Suppose  that  the  voltage  pulse  has  created  (or 
opened)  some  holes  of  approximately  circular  cross 
section  in  the  oxide.  Since  the  oxide  is  reasonably 
hydrophilic,  the  following  relationship  holds: 

p-ITJd,  (1) 

where  p is  the  gas  pressure  in  the  blister,  T,  is  the 
surface  tension  of  the  solution  against  the  gas,  and  d 
is  the  pore  diameter.  Using  T„  = 75  dyn/cm,  if 
d < 20  nm,  then  p must  exceed  about  150  atmos- 
pheres to  force  gas  bubbles  through  the  pores  into 
the  bulk  electrolyte.  Thus,  the  gas  remains  trapped 
in  the  blister  until  the  oxide  itself  ruptures,  at  which 
time  many  bubbles  appear. 

There  is  evidence  that  the  same  reactions  that  occur 
in  pitting  corrosion  also  proceed  during  blistering. 
Recall  that  blisters  will  not  grow  unless  the  alumi- 
num's electrode  potential  is  greater  than  the  pitting 
potential.  This  elementary  connection  between  blis- 
tering and  pitting  leads  one  to  suppose  that  all  the 
necessary  reactants  for  the  pitting  reaction  reach  the 
bare  metal  during  blister  growth.  An  important  step 
in  pitting  is  considered  to  be  the  acidification  of  the 
pit  electrolyte  due  to  the  hydrolysis  of  the  metal  ion. 
This  requires  water.  The  pathway  of  the  reactants  to 
the  reaction  site  in  this  hypothesis  must  be  through 
the  holes  and  along  the  inside  of  the  oxide  to  the 


periphery  of  the  blister,  which  then  becomes  the 
prominent  reaction  site  because  elsewhere  the  film 
and  the  electrolyte  are  separated  from  the  metal. 

Although  capillary  forces  probably  play  a role  in 
transporting  materials  along  the  inside  of  the  film  to 
the  reaction  site,  electrical  forces  could  also  be  effec- 
tive in  this  regard  because  the  blister  growth  rate 
increases  with  increasing  bias  potential.  The  pressure 
within  the  blister  may  also  be  a factor  in  breaking 
open  new  pores  near  the  periphery  as  the  blister 
grows.  As  a result  of  the  reaction  occurring  around 
the  perimeter,  the  adhesion  of  the  oxide  to  the  metal 
in  that  region  is  lessened,  allowing  the  blister  to 
grow.  It  does  so  in  an  irregular  incremental  fashion 
because  the  detachment  proceeds  in  a non-uniform 
spatial  manner  until  the  adhesion  weakens  sufficiently. 
The  gas  pressure  then  finishes  the  separation  and 
reestablishes  the  roughly  circular  morphology. 

Many  relevant  questions  remain  to  be  answered 
before  the  full  significance  of  the  blistering  phenom- 
enon can  be  appreciated.  Certainly  a demonstra- 
tion of  generality  to  other  types  of  localized  corro- 
sion processes  would  be  very  important.  Nevertheless, 
the  discovery  of  an  unknown  step  in  a particular 
process  is  significant  because  it  provides  a helpful 
insight  into  a complex  problem  together  with  a new 
direction  for  future  experimentation. 
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SEPARATION  FROM  ALL  CURVED  SURFACES 


Analytic  descriptions  of  steady  two-dimensional 
laminar  separation  from  cylindrical  surfaces,  concave 
or  convex,  reveal  universal  3:2:1  slope  relations  among 
a triplet  of  curves  associated  with  separation.  The 
curves  are  the  streamf unction  * - 0,  the  along-wall 
velocity  component  u - 0,  and  the  vorticity  f - 0. 


BACKGROUND 

The  separation  of  fluid  flow  from  the  bounding 
walls  is  a recurrent  problem  in  many  technical  areas. 
Conventional  approximate  boundary  layer  theory  of 
the  flow  process,  even  for  laminar  two-dimensional 
flow,  is  incomplete,  and  any  definitive  statement  is 
welcome.  The  demonstrated  regularity  of  two-dimen- 
sional viscous  flow  separation  (in  contrast  to  the  sin- 
gularities of  boundary  layer  theory)  generalizes  to  all 
curved  surfaces  the  result  first  enunciated  for  flat 
walls  by  the  German  aerodynamicist,  Oswatitsch,  in 
1957  (Ref.  1). 


DISCUSSION 

It  has  been  nearly  75  years  since  Ludwig  Prandtl 
(Ref.  2)  showed  that  the  separation  of  a viscous 
boundary  layer  from  a body  in  a flowing  fluid  made 
a significant  difference  in  the  expected  inviscid  pres- 
sure field  about  the  body.  Designers  of  aerodynamic 
control  surfaces,  turbine  manufacturers,  and  artificial 
heart  valve  specialists  are  but  a few  of  the  people 
who  are  concerned  about  predicting  and  understand- 
ing the  separation  phenomenon.  It  is  illustrated 
schematically  for  two-dimensional  laminar  flow  in 
Fig.  1.  Although  the  bulk  of  the  fluid  is  flowing 
smoothly  to  the  right,  there  is  a region  of  reversed 
flow  (leftward)  along  the  wall.  It  bends  outward 
and  ultimately  also  flows  to  the  right.  The  curve  that 
contains  this  turning  flow  is  the  separation  stream- 
line, + = 0,  and  the  point,  S,  on  the  wall  is  called 
the  separation  point.  Sometimes  the  separation  stream- 
line turns  downward  and  intersects  the  solid  wall  at 
a reattachment  point,  but  more  usually  the  free  shear 
layer  becomes  turbulent  if  the  flow  Reynolds  number, 
Re,  is  high,  (Re  is  l’L/v,  where  U is  the  oncoming 
flow  velocity,  L a characteristic  length,  and  t>  the 
fluid's  kinematic  viscosity. ) 

Prandtl  showed  that  the  criterion  for  two-dimen- 
sional separation  on  a fixed  wall  was  zero  wall  shear, 


or  (?u/Sy),  — 0,  where  u is  the  velocity  component 
parallel  to  the  surface  and  y is  the  coordinate  normal 
to  the  surface.  Figure  1 shows  sever.  1 velocity  pro- 
files, /<■()  ),  near  the  stagnation  point.  Also  shown  is 
the  curve  u = 0,  which  lies  entirely  within  the  sepa- 
rated flow  region. 

Within  the  boundary  layer,  viscous  effects  are  very 
important.  At  the  solid  smooth  wall,  the  viscous  con- 
dition of  nonslip  must  obtain  for  ordinary  incompres- 
sible flow.  Restricting  the  flow  analysis  to  a local 
region  near  the  separation  point  on  a flat  wall, 
Oswatitsch  (Ref.  1)  showed  the  angular  relationship 
between  the  three  important  curses:  the  two-dimen- 
sional separation  streamline  ♦ = 0,  the  curve  u = 0, 
and  the  curve  £ = 0,  where  the  vorticity,  £,  is  the 
magnitude  of  the  curl  of  the  velocity  vector  (£  = V 
X q,  q = »ij.  + <’{„).  At  S,  where  the  curves  all 
intersect  the  wall,  £ = ?«/0y,  so  this  is  in  accord 
with  the  Prandtl  criterion.  Near  S on  the  flat  wall, 
the  curves  £ = 0,  u = 0,  and  * ■=  0 have  slopes 
1:2:3. 

Although  the  practical  concern  for  separation  often 
is  related  to  high  midstream  Reynolds  number,  sepa- 
ration can  also  occur  at  very  low  Reynolds  number 
(so-called  Stokes  flows).  Two  such  cases  that  might 
occur  near  rotating  shafts  are  illustrated  in  Fig.  2. 
In  Fig.  2a  the  outer  cylinder  of  a pair  of  eccentric 
cylinders  is  rotating  with  constant  angular  velocity, 
T;  the  inner  cylinder  is  stationary.  Note  the  separa- 
tion streamline,  * = 0,  which  intersects  the  convex 
fixed  cylinder  at  the  separation  point,  5,  and  the  re- 
attachment point,  S'.  In  Fig.  2b  the  inner  cylinder  is 
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rotating  and  the  outer  cylinder  is  fixed.  The  stagna- 
tion points,  5 and  S',  now  occur  on  the  outer  concave 
cylinder.  The  exact  closed-form  analytic  expressions 
in  bicylindrical  coordinates  for  the  respective  Stokes 
flows  are  known,  and  the  details  near  S and  S’  have 
been  investigated  (Ref,  3). 

The  two  cases  illustrated  are  complementary,  for 
it  can  be  shown  that  on  the  convex  surface  (Fig.  2a) 


or  the  concave  surface  (Fig.  2b),  the  angular  rela- 
tionships between  each  of  the  three  intersecting  diag- 
nostic curves  and  the  wall  are  exactly  the  same.  The 
angle  of  departure  of  the  separating  streamline  de- 
pends on  the  relative  size  of  the  circular  cylinders 
and  the  eccentricity,  but  the  relationship  of  that  angle 
to  those  of  the  other  two  curves  is  the  universal  rule. 

Locally,  near  the  stagnation  point  (either  separation 
or  reattachment),  let  tj  = In  r/r{  (i  = 1,2)  and  let 
6 be  the  angle  measured  from  the  wall  surface  at  S. 
The  conformal  flow  domain  corresponding  to  Fig.  2a 
is  the  upper  part  of  Fig.  3;  that  corresponding  to  Fig. 
2b  is  the  lower  part,  as  indicated.  Either  way,  the 
diagnostic  curves  have  the  1:2:3  slope  relationship 
near  S (or  S').  Therefore,  the  concavity  or  convexity 
of  the  surface  has  no  effect  on  the  viscous  relation- 
ships near  the  stagnation  points.  Neither  does  the  real 
radius  of  curvature,  since  it  is  used  for  the  normalizing 
length,  L,  in  the  analysis.  Cavitation  can  occur  within 
the  lubricating  fluid  where  separated  regions  are  indi- 
cated in  incompressible  flow. 


Fig.  3 The  three  important  analytic  curves  near  sepa- 
ration. (i|,  8)  are  local  orthogonal  coordinates  (a  con- 
formal transformation). 


Also,  the  two-dimensional  nomalized  pressure 
field,  p',  is  orthonormal  to  the  voiticity  field,  £;  i.e., 

d/  = df.  d / = _ d£ 

dr)  d6'  dd  dr) 

where  p'  — pL/(pU).  Since  the  dimensional  wall 
shearing  stress,  r,  is  proportional  to  vorticity,  other 
relationships  at  separation  (Ref.  1)  generalized  to 
curved  surfaces  are 

f n v = (df/d?).  _ (dr/diy).  = ( [dp'/dd ), 

(df/dd).  (dr/dS),  (d p'/dv),  • V ; 

The  exact  Stokes  analysis  shows  how  the  location 
and  angle  of  the  departing  separation  streamline  from 
a given  rounded  surface  (concave  or  convex)  depend 
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on  the  relative  size  and  eccentricity  of  the  driving 
cylinder.  In  pressure  terms,  the  recirculation  region  on 
a curved  wall  depends  on  the  induced  pressure  gra- 
dients, not  on  the  wall  curvature  per  se.  The  angle 
of  departure  is  determined  by  the  local  ratio  of  normal 
and  tangential  pressure  gradients. 

As  Re  increases  from  nearly  zero,  the  flow  patterns 
in  the  eccentric  cylinder  flows  are  modified.  Within  a 
given  geometry,  one  can  expect  the  two-dimensional 
separation  streamlines  to  change  shape  and  the  attach- 
ment (and  reattachment)  points  to  migrate  from  those 
of  Stokes  flow;  yet  locally  near  the  stagnation  points 
where  the  convection  term  goes  to  zero,  the  balance  is 
still  between  dissipation  and  pressure  forces,  and  the 
local  flow  is  a Stokes  flow'.  Therefore,  on  all  curved 
boundaries  the  local  detaching/reattaching  two-dimen- 
sional flows  resemble  those  discussed  here;  only  the 
intersection  angles  can  change,  not  the  relationships 
among  the  characteristic  curves. 

The  analytic  singularity  found  near  the  separation 
point  by  ordinary  two-dimensional  boundary  layer 
theory  on  flat  or  curved  bodies  under  a prescribed  ex- 
ternal streamwise  pressure  gradient  signals  theoretical 


difficulty.  The  singularity  corresponds  to  a normal 
departure  of  the  separation  streamline,  inconsistent 
with  the  acute  angle  departure  usually  found  when 
the  total  momentum  equation  is  used.  The  analysis  of 
the  local  stagnation  region  must  incorporate  other 
local  flow  information  as  well  as  an  upstream  flow 
having  a similar  high  Reynolds  number  to  ensure  the 
universal  regularity  of  separating  laminar  flow  shown 
by  the  viscous  angular  relations. 
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EXPLICIT  Lz  INEQUALITIES  FOR 
PARABOLIC  EQUATIONS  WITH 
NEUMANN  BOUNDARY  CONDITIONS 


Explicit  integral  inequalities  are  presented  for  sec- 
ond-order diffusion  equations  with  Neumann  boundary 
conditions.  Such  inequalities  are  used  to  approximate 
solutions  to  partial  differential  equations  by  the  meth- 
od ot  a priori  inequalities. 


BACKGROUND 

Suppose  one  wishes  to  determine  the  distribution 
of  temperature,  u(_x,y,t),  over  a region,  R,  subject 
to  the  initial  condition, 

“OosO)  = Kx,y)  CO 


for  time  t = 0,  and  the  Neumann-type  boundary 
condition  for  which  the  rate  of  temperature  change 
on  the  boundary,  dR,  is  given  by 

- cow). 

dn 

C *,y)  effl,  o <t<T,  (2) 

where  d'dn  is  the  outer  normal  derivative  on  R. 
If  the  rate  at  which  heat  is  put  into  R is  H(.v,y,r), 
we  find  the  temperature,  uQx.y.t'),  by  solving  the 
second-order  parabolic  partial  differential  equation, 
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d2u  d2u  du 
dx2  ^ dy2  ~ ~dt 


~ H(x,y,t) 

» 


O ,y)  6 R,  0 < t <T, 


0) 


along  with  the  initial  conditions  (Eq.  1)  and 
boundary  conditions  (Eq.  2).  Equation  3 is  the 
heat  or  diffusion  equation. 

Except  when  R is  very  simple,  the  system  (Eqs.  1, 
2,  and  3)  cannot  be  solved  exactly  but  must  be  ap- 
proximated. The  most  commonly  used  methods  of 
approximation  are  the  finite-difference  method  and  the 
finite-element  method.  In  the  finite-difference  method, 
the  differential  operator  is  replaced  by  a difference 
operator  on  a network  of  points  over  R.  This  results 
in  a system  of  linear  equations  (necessarily  a large 
system  to  obtain  accurate  approximations)  that  is 
solved  algebraically  on  a computer. 


In  the  finite-element  method,  the  solution  is  sought 
in  the  form  of  the  minimum  of  a certain  functional. 
The  functional  is  minimized  over  a finite  set  of  ele- 
ments that  usually  are  piecewise  polynomials  over  a 
network  of  small  triangles  covering  R.  This  also  re- 
sults in  a system  of  linear  equations  that  is  solved  by 
computer.  The  number  of  equations  needed  for  accu- 
racy compared  to  the  finite-difference  method  is  gen- 
erally smaller,  but  the  savings  are  offset  somewhat  by 
the  need  for  more  preliminary  computations. 

We  wish  to  describe  an  alternate  method  of  ap- 
proximating the  solution  to  the  system  (Eqs.  1,  2,  and 
3)  that  is  not  as  well  known  as  the  above  methods 
but  which  has  more  convenient  features.  This  is  the 
method  of  a priori  inequalities. 


S2v  . d2v  dv 

Lv  ~ d?  + df  Jt  ' 

The  constants  «lf  at,  and  a 3 depend  only  on  the 
geometry  of  R and  are  independent  of  the  function 
v.  Explicit  values  for  on,  0%  and  a 3 have  been  ob- 
tained by  us  (Ref.  1). 

Once  we  have  the  a priori  inequality  (Eq.  4),  we 
may  take  any  approximation,  \p,  to  the  solution, 
u,  and  set  v = u — \p  in  Eq.  4 to  get  the  L2  error 
estimate: 


SSSc  («  — 'py  dxdydt  < «i  JJB  (F  — \py  dxdy 

+ J/a  (G-y„J  dsdt 

+ a3  Jffc  (H  - Lpy  dxdydt . (5) 

Notice  that  the  right-hand  side  depends  only  on  the 
known  quantities,  xp,  F,  G,  and  H.  In  particular,  we 
can  let 

m 

i diVi  , 

#-l 

a linear  combination  of  some  simple  trial  functions 
such  as  polynomials,  where  the  coefficients,  ah  are 
chosen  to  minimize  the  right-hand  side  of  the  equa- 
tion. 


DISCUSSION 

As  the  name  implies,  this  method  requires  that  we 
use  a certain  type  of  integral  inequality  related  to  the 
system  and  called  an  a priori  inequality  because  it 
is  to  hold  for  any  sufficiently  differentiable  function, 
v.  defined  on  the  region,  R.  The  appropriate  a priori 
inequality  to  use  for  the  system  is 

J7JC  v 2 dxdydt  < ai  J7fl  v2  dxdy 

+ a,//«(£)2^' 

+ “s  fffc  (£«0J  dxdydt  (4) 

where  C is  the  time  cylinder  R X (0,T);  B is  the 
base  of  the  time  cylinder  R at  / = 0;  and  S is  the 
side  of  the  time  cylinder  SR  X (0,F);  and  where 
v is  any  sufficiently  differentiable  function,  with 
the  notation 


CONCLUSIONS 

Using  the  method  of  a priori  inequalities  to  ap- 
proximate a solution  to  the  system  (Eqs.  1,  2,  and  3) 
also  results  in  a system  of  linear  equations  to  be  solved 
on  a computer.  The  number  of  equations  required  for 
a stated  accuracy  will  generally  be  even  fewer  than  for 
the  finite-element  method.  However,  the  trial  func- 
tions in  the  method  of  a priori  inequalities  do  not 
need  to  satisfy  any  particular  boundary  conditions, 
and  the  amount  of  preliminary  calculations  needed  is 
generally  less  than  for  the  finite-element  method.  The 
method  presented  here  should  be  more  efficient  nu- 
merically in  that  fewer  preliminary  calculations  and 
less  computer  time  will  be  needed  to  solve  a given 
problem  to  an  accuracy  comparable  to  those  of  the 
finite-difference  or  finite-element  methods. 

We  are  currently  in  the  process  of  implementing 
this  method  as  the  basis  of  a computer  software  sys- 
tem. 
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REFINEMENT  OF  THE  RADICAL-PAIR 
MECHANISM  OF  CHEMICALLY  INDUCED 
ELECTRON  POLARIZATION:  VECTOR 
MODEL  AND  AN  ANALYTIC  SOLUTION 


(t  has  been  shown  that  the  evolution  of  the  electron 
spin  state  of  a radical  pair  may  be  described  simply 
as  the  rotation  of  a three-dimensional  vector  repre- 
sentation of  the  spin  state,  where  the  rotation  rate 
and  axis  are  determined  by  the  magnetic  interactions 
within  the  radicals  and  by  the  electron-spin-dependent 
valence  (i.e.,  chemical  bonding)  interaction  between 
the  radicals.  The  model  yields  a clear  qualitative  pic- 
ture of  how  these  interactions  in  chemically  reacting 
radicals  yield  nuclear  spin  polarization  of  the  reaction 
products  and  electron  spin  polarization  of  the  radical s 
themselves  (Ref.  1).  Also,  an  analytic  solution  has 
been  obtained  to  the  Stochastic-Liouville  equations 
for  electron  polarization.  These  equations  are  ob- 
tained by  combining  the  vector  equation  for  temporal 
evolution  of  the  radical-pair  electron  spin  with  a dif- 
fusion term  describing  the  spatial  evolution  of  the 
pair  (Ref.  1). 


BACKGROUND 

The  radical-pair  model  of  chemically  induced  mag- 
netic polarization  (Ref.  2)  shows  that  the  polariza- 
tions result  from  the  combined  effects  of  the  interac- 
tions described  in  the  following  Hamiltonian : 

3C (r)  = -/(r)  (K  + 2S,  • S2) 

+ «C?iSi  • H + gt S2  • H) 

+ Ath  • Si  + /f2I2  • S2 . (1) 

The  first  term  in  Eq.  1 is  the  valence  interaction 
where,  as  shown  in  Fig.  1,  J(r)  is  a short-range  func- 
tion of  inter-radical  separation,  r,  and  Si  and  S2  are 
the  spins  of  the  two  radicals.  The  remaining  terms 
are  the  Zeeman  interaction  of  the  electron  magnetic 


Fig.  1 Typical  exchange  interaction  with  a rigid  sphere 
barrier  at  r — a. 


moments,  mtg.S,  and  upgSt.  ’’.'th  the  external  mag- 
netic field,  H.  and  the  electron- nuclear  hyperfine  in- 
teractions where  I,  and  I.,  are  nuclear  spins  on  radicals 
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1 and  2 and  A,  and  A.,  are  the  corresponding  hyper- 
fine  constants. 

These  nuclear-spin-dependent  magnetic  interactions 
can  mix  the  chemically  reactive  singlet  and  chemically 
unreactive  triplet  electron  spin  states  of  the  radical 
pair.  This  leads  to  nuclear  polarization  since  certain 
nuclear  spin  states  acquire  or  lose  singlet  character 
faster  than  others  and  are  correspondingly  more  or 
(ess  abundant  in  the  radical-pair  reaction  products. 
This  singlet-triplet  mixing  combined  with  the  spin- 
selective  valence  interaction  also  leads  to  electron  po- 
larization; however,  no  suitable  qualitative  picture  of 
this  process  was  available  until  the  development  of 
the  vector  model  described  in  this  report. 

Diffusion  plays  an  essential  role  in  the  radical-pair 
polarization  mechanism.  Specifically,  it  enables  some 
pairs  to  separate  and  subsequently  reencounter  each 
other  after  an  intervening  period  of  diffusion  that  is 
long  enough  to  allow  for  significant  singlet-triplet  mix- 
ing by  the  weak,  slowly  acting  magnetic  interactions 
(Refs.  2 and  3).  Quantitative  treatment  of  diffusion 
is  difficult  because  it  requires  the  time-dependent 
Schrodinger  equation  to  be  solved  for  the  radical-pair 
spin  function  with  J(r)  a random  function  of  time 
because  of  the  diffusive  motion.  Previous  treatments 
of  this  problem  were  either  highly  approximate  or 
required  numerical  solution  and  none  gave  a really 
clear  qualitative  picture  of  the  electron  polarization 
process. 


DISCUSSION 

The  radical-pair  spin  state  is  conveniently  de- 
scribed by  a density  matrix  whose  components  are 
PiiC.0  ~ where  the  C's  are  the  coefficients  of 

the  singlet  and  triplet  states  of  the  radical  pair 
wave  function,  i.e.,  yf>ap  = CsCOlJ-)  + Cr(/)]T). 
Thus,  p is  a 2 X 2 matrix  where  Cs  and  Cp  are 
complex  quantities.  yf/Rp  is  determined  by  the  time- 
dependent  Schrodinger  equation  ’X.yj/Rp  = idyf/Rp/dt, 
where  3C  is  given  by  Eq.  1;  from  this  it  can  be 
shown  that  X p — pX  = idp/dt. 

It  is  convenient  to  transform  p so  that  its  com- 
ponents have  the  following  real  and  physically 
significant  forms  at  all  times  (since  the  radical  pair 
is  always  a singlet,  a triplet,  or  some  combination 
thereof);  po  = pss  + Ptt,  x = par  + Pts,  y = 
—Kpst  — Pts),  z = pss  — Ptt,  po  = 1.  The  term 
z gives  the  relative  amounts  of  singlet  and  triplet 
character.  It  can  be  shown  that  x is  the  electron 
spin  polarization  of  the  pair,  where  polarization 
means  an  excess  of  "up”  spin  on  one  radical  and 
a matching  excess  of  "down"  spin  on  the  other. 


In  the  absence  of  diffusion,  the  time  evolution  of 
the  components  of  this  density  matrix  is  given  bv 
the  following  simple  vector  equation : 


(2) 


Fig.  2 Vector  model  of  chemically  induced  nuclear  polar- 
ization (CIDNP)  and  chemically  induced  electron  polarization 
(CIDEP).  (a)  S - Tq  mixing  quenched  when  ) >>«.  (b)  S - 
To  mixing  when  a (c)  CIDEP  results  from  exposure  to  a 
sizable  exchange  interaction,  J > a,  after  a period  of  S~T0 
mixing  as  shown  in  (b).  (d)  Electron  polarization  destroyed 
by  spin  exchange  if  radical  pair  remains  too  long  in  region  of 
large  J. 
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Here,  X denotes  the  vector  cross  product  and  a is  the 
singlet-triplet  mixing  rate  given  by 

J “ V’PfiH (gx  - gt)  + \/lAlMx  - l/yAjM, 
where  dl,  and  At,  denote  the  spin  states  of  the  nuclei. 

According  to  Eq.  2,  the  rate  of  change  of  p at  any 
time  is  a rotation  of  magnitude,  about  an 

axis  in  the  direction  of  fj.  As  depicted  in  Pig.  2a, 
when  the  radicals  are  close  together,  / is  much  greater 
than  .«  and  the  rotation  of  p is  about  the  z axis;  thus, 
the  2 component,  which  determines  the  relative 
amounts  of  singlet  and  triplet  character  of  the  pair, 
is  unvarying. 

As  the  radicals  separate,  /(  r)  approaches  0,  and  p 
begins  to  rotate  above  the  v axis,  as  shown  in  Pig.  2b. 
Phis  changes  the  2 component  anil  corresponds  to 
singlet  triplet  mixing  at  rate  a.  However,  this  process 
alone  cannot  give  p a component  along  the  electron 
polarization  axis,  ,v.  As  shown  in  Pig  2c,  an  x com- 
ponent can  be  achieved  from  the  usual  starting  point 
of  a pure  singlet  or  triplet  state  ( i.e  , p - (0,0,  ± 1 ) ) 
if,  after  a period  of  singlet-triplet  mixing,  the  radicals 
come  close  together  again,  whereupon  the  p vector 
(which  is  now  in  the  )z  plane)  begins  to  rotate  about 
the  2 axis  at  rate  J(r),  thereby  acquiring  an  x compo- 
nent. 


Finally,  as  shown  in  Fig.  2d,  if  the  radicals  ap- 
proach too  closely,  ](r)  becomes  very  large  and  p 
makes  a number  of  complete  revolutions  about  the 
2 axis,  thereby  yielding  an  average  value  of  zero  for 
the  x component.  This  last  effect,  known  as  spin  ex- 
change, occurs  because  the  electrons  become  delocal- 
ized between  the  strongly  interacting  radicals  and  any 
electron  spin  polarization  on  one  radical  is  distributed 
between  the  two  radicals.  Since  the  polarization  in  our 
case  is  an  excess  of  up  spins  on  one  radical  and  a 
matching  excess  of  down  spins  on  the  other,  the 
polarization  is  averaged  to  zero  by  the  spin  exchange 
process. 


Equation  2 is  converted  to  the  Stochastic-I.iouvillc 
equations  by  adding  a diffusion  term  (Ref.  4): 

3*0,0 


dt 


D 


m 


+ 9(0  X p 


CO 


Since  we  are  interested  in  the  steady-state  (t -*  <x>) 
polarization,  the  time  dependence  can  be  eliminated 
from  this  equation  by  taking  the  Laplace  transform 
and  using  the  relation 

lint  /(r)  = lint  r/0) 

f—"  i-«0 


where  /(. i)  is  a Laplace  transform.  The  resulting  sys- 
tem of  coupled  differential  equations  can  be  trans- 
formed into  a single  integral  equation  for  the  electron 
polarization  component  v.  For  slow  singlet -triplet 
mixing,  which  is  the  usual  case,  the  kernel  of  this 
integral  equation  is  positive-definite,  symmetric,  and 
compact.  The  kernel  therefore  has  a complete  set  of 
eigenfunctions  that  happen  to  lx-  zero-order  Bessel 
functions,  and  the  integral  equation  may  lx-  solved 
by  expanding  in  these  eigenfunctions.  The  solution 
for  the  usual  case  of  a strong,  exponentially  decaying 
exchange  interaction,  J(r)  - /„<■  Xr,  has  the  simple 
form: 


Polarization  = — 

This  result  agrees  well  with  previous  numerical  cal- 
culations  (Ref.  ">). 
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VALENCE-BOND  STUDY  OF  HYPERFINE 
INTERACTIONS  AND  STRUCTURE  OF 
THE  NOBLE  GAS  MONOHALIDES 
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It  has  bean  shown  that  a samiempirical  valence- 
bond  (VB)  model  of  the  noble  gas  monohalides  that 
includes  electron  correlation  effects  can  account  tor 
the  observed  electron-nuclear  magnetic  hypertine 
structure  (hfs)  interactions  in  these  molecules  (Ref.  1). 
The  VB  wave  function  is  a linear  combination  of  neu- 
tral and  ionic  VB  structures,  that  is, 

+nx(*2)  - x*(N  *•)  + *("'•  X ) 0) 

when  N denotes  the  noble  gas,  X is  the  halogen,  x is 
the  charge  distribution  parameter,  and  typically 
VI  - x1  i*  much  less  than  1 . The  dot  denotes  the 
location  of  the  unpaired  electron  in  a given  VB  struc- 
ture, Electron  correlation  is  included  in  this  wave  func- 
tion by  using  different  atomic  orbitals  for  a neutral 
atom  and  the  corresponding  ion,  and  by  considering 
the  polarizing  effects  of  van  der  Waals  interactions  in 
the  neutral  VB  structure  (N  X ) and  coulomb  inter- 
actions in  the  ionic  structure  (N'  - X ).  The  bond  dis- 
tance, ft,  and  x in  fhe  molecules  XeF,  KrF,  and  XeCI 
are  determined  as  those  values  that  give  the  best 
agreement  between  the  observed  and  the  calculated 
hypertine  structure  constants  (hfc). 


BACKGROUND 

The  noble  gas  monolulides  are  a very  interesting 
example  of  a class  of  molecules  known  as  hypervalcnt 
radicals,  so  called  because  they  contain  an  unpaired 
electron  that  cannot  lx-  accommodated  by  VB  struc- 
tures satisfying  the  classical  I.ewis  octet  rule.  I bis  rule 
predicts  that  noble  gas  atoms  will  not  form  chemical 
bonds  because  they  already  have  a filled  outer  shell 
of  eight  electrons.  In  addition  to  their  importance  in 
terms  of  chemical  bonding  theory,  their  strong  ultra- 
violet emission  from  an  excited  state  corresponding 
to  the  (N'-  X ) VB  structure  gives  them  practical 
importance  as  lasers. 

The  unpaired  electron  of  these  modecules  makes 
them  paramagnetic;  consequently,  much  important 
structural  information,  especially  the  hfs  interactions, 
can  be  obtained  from  experiments  in  electron  spin 
resonance  spectroscopy.  Transitions  are  observed  be- 
tween energy  levels  associated  with  the  interaction  of 
the  electron  spin  magnetic  moment  with  an  external 
magnetic  field  and  with  the  magnetic  nuclei.  In  order 
to  translate  the  observed  hfc  into  information  about 
the  electron  charge  distribution  and  chemical  bonding, 
it  is  necessary  to  construct  a model  wave  function  for 
the  system. 


DISCUSSION 

The  hyperfine  data  for  a given  nucleus  consist  of 
an  isotropic  and  an  anisotropic  hfc,  denoted  ,i  and  Rt 
respectively  (Refs.  2 to  4).  The  value  of  ,/  is  deter- 
mined by  the  unpaired  electron  density  in  the  spheri 
cally  symmetric  s orbitals  of  the  atom  to  which  the 
nucleus  belongs,  whereas  R is  determined  by  the  un- 
paired electron  density  in  the  nonsphcrically  sym- 
metric atomic  orbitals,  in  this  case  the  p orbitals.  In 
molecules  such  as  the  noble  gas  monohalides,  which 
contain  heavy  atoms,  the  spin-orbit  interactions  arc 
large  enough  that  the  hfc  contain  a contribution  from 
the  orbital  magnetic  moment  of  the  electron  as  well 
as  from  its  spin  moment.  It  is  fairly  straightforward 
to  determine  the  pure  spin  hfc  by  subtracting  an  esti- 
mate of  the  orbital  contribution  from  the  experimental 
values  (Ref.  1). 

The  calculation  of  R,  from  the  wave  function  of 
Eq.  1 is  relatively  straightforward  because  the  un- 
paired electron  is  primarily  fixated  in  valence  f>  or- 
bitals of  the  noble  gas  cation  N'-  and  the  halogen 
atom  X-  Accordingly,  R v is  approximately  propor- 
tional to  ( 1 — x’).  which  is  the  weight  factor  of  the 
(N*-  X ) structure  in  Eq.  1,  whereas  B j is  ap- 
proximately proportional  to  y2,  which  is  the  weight 
factor  of  the  (N  X-)  structure.  There  is,  however, 
an  additional  contribution  to  Kith  R , and  R >■  from 
a cross  term  involving  the  two  VB  structures.  This 
term,  whose  weight  factor  is  y V 1 — x‘  (where  S 
is  an  overlap  integral  with  an  order  of  magnitude  of 
0.1  to  0.2),  is  especially  important  for  B.  v because 
the  direct  term  proportional  to  1 — y-  **  usually 
small.  Thus,  the  parameter  y can  be  determined  with 
good  accuracy  from  the  anisotropic  hfc  alone. 

The  origin  of  the  isotropic  hfc  is  more  subtle  lie- 
cause,  at  the  level  of  approximation  denoted  by  Eq.  1, 
neither  of  the  VB  structures  provides  an  unpaired 
electron  density  in  the  atomic  .t  orbitals.  The  .i-orbital 
density,  and  consequently  the  isotropic  hfc,  results 
from  the  following  refinements  to  the  model  depicted 
by  Eq.  1. 

Overlap.  In  the  VB  structure  (N  X - ).  a small 
amount  of  the  unpaired  electron  density  in  the  valence 
pz  orbital  of  X (i.c.,  a p orbital  directed  along  the 
molecular  bond)  is  transferred  to  the  > orbitals  of  N 
via  the  overlap  of  the  X,,.  and  N,  orbitals.  Similarly 
the  overlap  of  the  N'p;  and  X , orbitals  in  the  VB 
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structure  (N'  ■ X ) transfers  unpaired  electron  dens- 
ity from  the  p orbital  of  the  noble  gas  cation  to  the 
halogen  i orbital. 

Polarization.  In  the  neutral  VB  structure 
(N  X-),  the  polarizing  interactions  involve  the 
effect  of  the  instantaneous  dipole  moment  of  one  atom 
on  the  other.  This  van  der  Waals  interaction  exists 
even  though  the  average  dipole  moment  of  the  atoms 
(obtained  by  averaging  over  the  instantaneous  mo- 
ments associated  with  various  electron  positions)  is 
zero.  In  the  ionic  VB  structure,  the  polarization  is  due 
to  the  coulomb  interaction  between  the  electrons  of 
one  ion  and  the  charge  on  the  other.  The  polarization 
of  the  electronic  charge  clouds  of  the  open-shell 
species  N*-  and  X-  is  spin  selective  and  leads  to  the 
following  isotropic  hfs  interactions.  First,  the  polari- 
zation of  the  pz  orbitals  containing  the  unpaired  elec  - 
tron involves  a deformation  of  these  orbitals,  which 
gives  them  a small  amount  of  spherically  symmetric 
or  r-orbital  character.  Second,  the  polarization-induced 
deformation  of  the  closed-shell  s orbitals,  occupied  by 
a pair  of  electrons  of  opposite  spin,  is  spin  dependent, 
thereby  creating  unpaired  electron  density  in  these 
orbitals.  This  occurs  because  the  Pauli  exclusion  prin- 
ciple, which  states  that  electrons  of  the  same  spin  may 
not  occupy  identical  orbitals,  requires  that  the  defor- 
mation of  the  s orbital  with  the  same  spin  as  the  pz 
orbital  not  give  it  any  of  the  shape  characteristics  of 
the  pz  orbital,  whereas  the  deformation  of  the  other 
s orbital  is  unrestricted. 

Overlap-Polarization.  This  is  a contribution  to 
the  isotropic  hfs  from  a combination  of  the  aforemen- 
tioned overlap  effect  and  the  polarization  of  the 
closed-shell  r orbitals. 

Atomic.  This  is  the  isotropic  life  of  the  isolated 
N‘-  ion  and  the  X-  atom.  The  hfs  term  results  from 
spin-dependent  interactions  between  the  closed-shell 
s orbitals  and  the  unpaired  electron  in  the  valence  pz 
orbital.  The  atomic  hfc  are  known  from  experiment 
for  the  halogen  atoms  and  have  been  estimated  for 
the  noble  gas  cations  by  extrapolation  from  the  halo- 
gen values. 

The  calculated  results  are  illustrated  for  XcF  in 
Fig.  1,  which  gives  the  isotropic  Xe  and  F hfc  as 
functions  of  R and  y.  As  expected,  the  calculated 
values  of  hfc  decrease  rapidly  with  increasing  inter- 
atomic distance,  which  reduces  the  overlap  and  polari- 
zation interactions.  Comparison  of  the  calculated  and 
experimental  isotropic  hfc  for  XeF,  also  shown  in 
Fig.  1 , shows  that  theory  and  experiment  agree  well 
for  R “ 4.70a„  and  y — O.H(X).  R and  y were  simi- 
larly determined  for  KrF  and  XeCI;  all  results  and 
a comparison  of  theory  and  experiment  are  given  in 
Table  1.  The  resulting  values  of  interatomic  distance 


0 0.5  10 

X 


Fig.  I Isotropic  hyperfine  constants  for  XeF  as  func- 
tions of  x and  R. 


are  consistent  with  the  relative  sizes  of  the  atoms,  i.e., 
Xe  > Kr  and  Cl  > F.  The  values  of  y show  that 
the  weight  of  the  (Nf-  X ) state  is  greater  in  XeF 
than  in  either  KrF  or  XeCI.  It  is  reasonable  that  the 
Xe'F  state  plays  a greater  role  in  XeF  than  does 
Kr*F-  in  KrF  because  the  ionization  potential  of  Xe 
is  less  than  that  of  Kr.  Similarly,  the  Xe’F  state  is 
more  important  in  XeF  than  is  the  Xe'CJ  state  in 
XeCI  because  the  F anion  is  considerably  smaller  than 
the  Cl  anion;  hence  the  coulomb  stabilization  of 
Xe*F  is  greater  than  that  of  Xe‘CI  . 
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COMPARISON  OF  COMPUTED  HYPERFINE  CONSTANTS  OF  THE 
NOBLE  GAS  MONOHALIDES  WITH  EXPERIMENTAL  VALUES 


R 

r 

4<Tk> 

B 

B (Mur) 

Molecule 

Ui> 

X 

Nucleus 

(MHz) 

(MHz) 

(MHz) 

(MHz) 

‘-••Xc 

1310 

1330 

1227 

1290 

XeF 

4.70 

0.800 

F 

509 

52} 

2125 

2126 

KrF 

4.00 

0 98 

F 

552 

549 

2925 

309) 

u.Xc 

270 

265 

425 

424 

XeC! 

6.  SO 

0 955 

•’•'Cl 

37 

37 

260 

305 

Note'  Th  - theory,  Exp  — experimental 


3.  W E.  Falconer.  J R.  Morton,  and  A.  C».  Streng.  "Electron 
Spin  Resonance  Spectrum  ot  KrF."  /.  Cktm  Pt>yj.  41.  1964. 
p.  90 2. 

4.  F.  J.  Adrian  and  V.  A.  Bowen.  *'ESR  Spectrum  of  Xetl  in 
Argon  at  4.2K."  /.  C.btm.  Pkyt.  65.  1976.  p.  4316. 
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LIGHT-SCATTERING  MEASUREMENT  OF 
PARTICLE  SIZE  DISTRIBUTIONS 


Intensity  correlation  spectroscopy  (ICS)  is  the  most 
rapid  and  accurate  method  ot  obtaining  macromolec- 
ular  diffusion  coefficients . In  addition,  it  provides  valu- 
able measures  of  sample  polydispersity.  A method  has 
been  developed  to  extract  size  distribution  parameters 
from  ICS  data  and  is  applied  to  data  on  very  low 
density  lipoproteins  (VLDL)  from  human  serum. 


BACKGROUND 

One  goal  of  our  present  research  is  to  determine  the 
size  distribution  of  various  classes  of  human  scrum 
lipoproteins  in  certain  patients.  We  a'e  examining 
these  lipoproteins  using  ICS.  Molecular  diffusion  coef- 
ficients and  therefore  particle  sizes  are  determines! 
from  the  autocorrelation  function  of  the  time-depend- 
ent intensity  fluctuations  in  the  light  scattered  from 


the  particles  as  they  diffuse  in  solution.  This  autocor- 
relation function  can  he  analyzed  by  an  expansion 
technique  in  which  the  coefficients,  called  cumulants, 
provide  a general  description  of  sample  polydispersity 
(Ref.  I).  If  the  form  of  the  size  distribution  is  speci- 
fied, then  the  values  of  the  parameters  that  charac- 
terize this  form  can  be  obtained  directly  from  the 
cumulants  (Refs.  2 to  4).  Although  this  approach  is 
less  general,  it  is  of  great  potential  use  in  our  re- 
search. 

Our  preliminary  studies  of  VLDI.,  which  are 
spherical  particles  with  diameters  ranging  from  30  nm 
to  more  than  90  nm,  indicated  that  the  Gaussian  dis- 
tribution we  deduced  from  the  cumulants  was  phys- 
ically unrealistic  and  that  a rectangular  distribution 
was  impossible  (Ref.  2).  This  finding  established  a 
clear  need  to  extend  our  method  to  other,  more  useful, 
size  distributions. 


DISCUSSION 

The  general  form  of  the  intensity  correlation 
function  from  a system  of  scattcrers  is  C(r)  = 

A 5 + fl.  where  A,  B,  and  r are  the  ampli- 

tude, background,  and  delay  time,  respectively , 
For  a single-sired  molecular  species,  g"’(r)  = 
exp (— /ty!r),  where  the  translational  diffusion 
coefficient,  />,  is  inversely  proportional  to  the  sire 
of  the  scattcrers.  The  factor  ^ is  an  instrumental 
parameter  related  to  the  laser  wavelength  and 
scattering  angle 

In  polydisperse  systems,  is  a sum  of  many 

such  exponentials.  Gumulznt  analysis  provides  the 
most  general  description  of  the  autocorrelation 
function  in  such  systems.  In  theory,  the  cumulants, 
k„,  arc  found  from  the  Maclaurin  scries  for  In 

Xm(r): 

lng‘l’0)  = In  A 4-  Y . 

,-i  m 

In  practice,  the  cumulants  arc  computed  approxi- 
mately from  a truncated  scries,  where  m is  limited 
to  2 or  ) due  to  scatter  in  the  experimental  data 
The  first  cumulant  is  related  to  the  ^-average 
diffusion  coefficient  of  the  scattcrers  hv  ki  = 
(D, ) if1,  and  k j is  related  to  the  normalized  variance 
of  the  distribution  of  diffusion  coefficients  by 
k;  kr  = «D.*>  - </>,>’)  (O.l*  (Ref  1). 

In  addition  to  providing  this  general  characteriza- 
tion of  polydispersity,  we  have  shown  that  cumulants 
can  also  be  used  to  find  the  parameters  of  postulated 
distributions  of  model  size.  Three  statistical  functions 
that  were  useful  in  characterizing  VLDL  in  this  way 
are  known  as  the  Pearson  V,  the  Pearson  III,  and  the 
Schulz  distributions  (Ref.  4).  Path  of  these  functions 
has  the  form  of  a skewed,  bell  shaped  curve  (Fig.  1) 
and  is  a function  of  parameters  that  arc  related  rather 
simply  to  such  quantities  of  interest  as  average  size 
and  the  standard  deviation  of  sizes.  Using  the  explicit 
expressions  for  the  first  two  cumulants,  K,  and  k%,  we 
also  showed  that  the  ratio  k'.VkV  has  distinct  upper 
bounds  for  each  of  these  distributions  as  well  as  for 
the  Gaussian  and  rectangular  distributions  treated  in 
Ref.  3.  This  property  can  often  be  used  to  rule  out  a 
distribution,  as  in  the  case  of  the  rectangular  distribu- 
tion for  VLDL.  The  usefulness  of  the  method  for  de- 
scribing VI.DI.  sizes  is  illustrated  in  Fig.  1.  The  three 


Fig.  I The  experimentally  determined  size  distributions 
tor  VLDL.  The  expected  range  ot  VLDL  sizes  based  on 
EM  studies  is  shown  on  the  baseline. 


size  distributions  compare  very  favorably  with  the  re- 
ported electron  microscopic  (EM)  size  range  of 
VLDL.  Experiments  are  presently  under  way  to  com- 
pare the  size  distributions  obtained  with  EM  and  ICS 
measurements  on  the  same  samples. 
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JHU  EVENING  COLLEGE  CENTER  AT  APL 


THE  APL  GRADUATE  EDUCATION  CENTER 


During  lis  thirty- five  years  of  existence,  API.  has 
used  various  methods  to  help  its  staff  members  main- 
tain and  improve  skills  needed  to  carry  out  their 
technical  assignments  in  the  most  effective  manner 
The  met  hods  include  financial  support  for  pertinent 
pn’grams  at  educational  institutions  and  specialized 
short  courses  taught  both  inside  and  outside  the  Lab- 
oratory. 

In  1964.  m order  to  make  some  of  APL's  m house 
courses  acailable  to  all  qualified  persons  regardless  of 
whether  they  were  employed  by  the  Laboratory, 
the  Fvening  C ollege  of  The  Johns  Hopkins  Univer- 
sity and  APL  cooperated  to  establish  a Graduate  ( en- 
ter at  API.  oriented  to  the  educational  needs  of  em- 
ployed scientists  and  engineers 

Under  the  academic  sponsorship  of  the  Evening 
College  and  with  faculty  members  drawn  primarily 
from  APL's  Senior  Staff,  the  ( enter  offers  a variety  of 
courses  in  scientific  and  engineering  areas  ( urrent 
programs  lead  to  master  of  science  degrees  in  Applied 
Physics.  Computer  Science,  Electrical  Engineering. 
Numerical  Science,  and  Space  Technology.  Classes 
meet  in  the  R E.  Cobson  Library  (Fig.  1 1 . 


Fig.  I The  R.  E.  Gibson  Library  where  Graduate  Ou- 
ter classes  meet. 


During  Fiscal  Year  1977.  class  registrations  lor 
API.  ( enter  programs  were  as  follows:  Fall  1976. 
748;  Spring  1977,  700;  Summer  1977,  224;  and  Fall 
1977,  868. 

Because  some  students  take  more  than  one  course, 
the  number  of  individual  students  is  slightly  less  than 


the  number  of  class  registrations,  l or  example,  in  the 
Fall  of  1977,  671  individuals  made  up  the  868  class 
registrations  shown.  Of  these.  338  1870  y were  not 
APL  staff  members  The  individuals  were  employed 
as  follows: 


Private  industry 

234 

Government  institutions 

17.3 

Eduiational  institutions 

71 

Not  employed  or  did  not  reply 

40 

It  is  evident  that  the  Center  is  a substantial  resource 
in  the  technical  and  educational  community. 

Eighty-four  persons  completed  their  programs  and 
received  master's  degrees  in  1977.  Seven  hundred  and 
forty-seven  APL  Center  students  have  received  mas- 
ter's degrees.  Of  that  number,  139  were  APL  staff 
members  and  608  were  not.  These  individuals  con- 
stitute 720  of  all  those  who  received  master  of  sci- 
ence degrees  in  technical  areas  from  the  various  JHU 
Evening  College  centers  during  that  time.  The  num- 
ber of  recipients  in  each  area  of  specialization  is: 


Numerical  Science 

280 

Electrical  Engineering 

271 

Computer  Science 

137 

Applied  Physics 

34 

Space  Technology 

25 

When  the  Center  1 ->egan  operation,  all 

courses  were 

taken  from  the  existing  JHU  curriculum.  In  Fiscal 
Year  1977,  however,  only  160  of  the  offerings  were 
derived  from  that  curriculum.  In  fact,  most  courses 
(620)  were  originated  and  taught  by  APL  staff  mem- 
bers. Another  160  were  created  expressly  for  the 
Center  by  JHU  faculty  members  of  the  Baltimore 
campus,  and  60  were  originated  by  Center  instructors 
who  are  primarily  employed  by  other  organizations. 

In  our  continuing  ctfort  to  maintain  the  leadership 
that  the  APL  Center  has  attained,  we  arc  constantly 
analyzing  curricula  and  making  revisions.  In  Fiscal 
Year  1977,  APL  Center  instructors  originated  and 
taught  six  courses  for  the  first  time:  Detection  and 
Estimation  Theory.  Hardware  Architecture  of  Digital 
Systems,  Digital  Simulation,  Computer  Systems  Engi- 
neering Management  Applications.  Microprocessor 
Applications,  and  Fault  Tolerant  Digital  Systems. 
Three  new  courses  are  planned  for  Fiscal  Year  1 9~8. 


Figure  - (incloses  thjt  records  were  set  in  Fiscal 
Year  tl)77  for  class  roll  totals,  number  of  individuals 
enrolled,  number  of  outside  participants,  number  of 
participants  from  private  industry,  number  of  partici- 
pants from  Government  facilities,  and  number  of 
candidates  for  master's  degrees  offered  at  API.. 


Although  the  physical  facilities  arc  modern  and 
the  audiovisual  aids  arc  excellent,  the  primary  strength 
of  the  (enter  is  its  faculty.  (Figure  } shows  a typical 
class.)  Some  full-time  faculty  members  from  the 
Baltimore  campus  instruct  in  this  program,  but  70G 
of  the  tenter's  faculty  are  members  of  the  API.  Senior 
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Fig.  2 Data  on  the  APL  Graduate  Education  tenter. 


Stiff  and  15%  arc  members  of  similar  staffs  of  other 
scientific  or  engineering  organizations.  The  students 
appreciate  and  appear  to  benefit  from  studying  under 
academically  qualified  instructors  who  arc  actively  in- 
volved with  the  subjects  they  teach  and  who  under- 
stand the  needs  of  part-time  students  who  work  in 
research  and  development  organizations. 


Fig.  ) A typical  class  in  session. 


Author:  P.  B.  Edwards 


PATENTS 


PATENTS  ACTIVITIES 


lie 


The  API.  Patents  Office  is  responsible  for  ensuring 
compliance  with  contract  and  grant  requirements  rcla 
live  to  patent  and  data  rights,  as  imposed  by  the  vari- 
ous governmental  agencies  that  sponsor  work  at  the 
Laboratory.  In  addition  to  preparing  formal  disclo- 
sures of  inventions  for  the  appropriate  sponsors,  the 
Patents  Office  prepares  and  prosecutes  patent  applica- 
tions on  behalf  of  both  the  University  and  the  Navy 
L>epartment 

The  following  lists  indicate  the  invention  disclo- 
sures submitted  to  sponsors,  the  patent  applications 
prepared  and  filed  in  the  United  States  Patent  Office, 
and  the  previously  filed  applications  that  were  success- 
fully prosecuted  to  issuance  as  patents,  during  Fiscal 
Year  1977, 


INVENTION  DISCLOSURES 

Minimal  Distortion  ViJeo  bandpass  Filter 
F N.  Sanione 
Quick  Release  Hook 
/ D Sternberg 

Trigger  Circuit  with  Noise  Controlled  Adaptive  Threshold 
for  Detection  of  Tone  Hurst  Signals 
/.  H Kick 

Circuit  for  Sorting  Direct  Path  Signals  from  Reflected  Signals 
/ H Kuck 
Pattern  Match  Circuit 
J H Kick 

Tone  Burst  Signal  Detection  System 
/ H KuJt 

Submersible  Non-Explosive  Acoustic  Generator  (SNAG) 

G.  S Bitnngs 

Cell  for  Measurement  of  Passive  Membrane  Transport  Co- 
efficients 

R A Meyer  and  M.  H.  Friedman 

Oblique  Folding  Wing 

R M.  Rttello  and  A S.  Polk 
Velocimeter  with  Acoustical  Sensor 
G 5.  Bitting i and  VT.  J.  Eresian 

Optical  Detector 

G.  S.  Binmgs , W J.  Eresian,  and  R.  E.  I X’alker 
Toroid  Turns  Tester 
F.  T.  Mar  cell  mo 

RF  and  IF  Circuitry  of  an  Uplink  Receiver 

R.  C.  Gipe 

Bonding  Agent  for  HMX 

/.  F Kincaid  and  R Reed 

Low  Level  (Dry)  Circuit  Continuity  Tester 
F.  T Marcellmo 

Switching  Voltage  Regulator 
D.  Qualkinbush 

Tape  Recorder  Sequencer 
/.  5.  Lombardo 


Programmable  To  ic  Burst  (PTB)  Generator 

S.  R.  Osborne 

Planned  View  Display  (PVD)  Raster  Scan  Generator 
5.  R Osborne 
Fire  Detector- Extinguisher 
R Af  Fn  strom 

Programmable  Mechanical  Associative  Processor 
L R.  Gieszl 

Clock  Invariant  Synchronization  for  Random  Binary  Se- 
quences 
C.  Philip pidet 

Harmonic,  Digital  Sine  Generator 
5.  R.  Osborne 

Programmable  Linear  FM  (LFM)  Generator 
5 R.  Osborne 
Thermistor  Chain 
G . 5.  Bttnngs 

Time  Transfer  Receiver  RF  Circuitry 

G . G W"  hit  worth 

Digital  AGC  Amplifier 
/ H Kick 

Method  of  Fabricating  Toroid  Inductors  and  Transformers 
LTsing  Torotdalgrams 
/.  F Adams,  Jr. 

Code  Converter 

P.  Klose  and  T.  G.  Calhoon 

Method  of  Determining  Sea  Forces  Acting  on  a Hovering 
Submarine 

5.  A.  Fogel,  W”.  F.  Knraua.  1.  D.  Rapport,  P.  C.  Smith, 
and  J.  J.  I X'osmak 

Firing  Control  Circuit  for  Resistance  Welder 
F.  T.  Marcellino 

Fire  Atmosphere  Badge  Structure 

R.  M.  Fnstrom,  D.  Shapiro,  and  V.  Esch 

High-Speed,  Glitch-Free,  Digital -to- Analog  Converter 
P.  J.  Grunberger 

Integral  Rocket -Scram  jet 
/.  L.  Ketrsey 

Time  Optimal  Geomagnetic  Spacecraft  Maneuvering  Tech- 
nique 

B.  E.  T os  smart 

Recovery  of  a Remotely  Piloted  Vehicle  (RPV)  Along  a 
Kite  Suspended  VC' ire 
M L.  Hill 

Vertical  Relative  Landing  of  a Remotely  Piloted  Vehicle 
(RPV)  lTsing  Wire  Terminal  Guidance 
M L.  Hill 

Net  Retrieval  of  a Remotely  Piloted  Vehicle  (RPV) 

AL  L.  Hill 

War- At-Sea  Graphical  Analysis  Model  (WASGRAM)  Game 
Including  an  Air  Battle  Package 

T.  P.  Mod  el  ski,  R.  A.  Freeman,  and  Jen-Yih  W "ang 

Vacuum  Probe  for  Water  Pickup  and  Handling 
J.  L.  Abita  and  H.  K.  Charles,  Jr. 

Miniature  Probe  Block 

H.  K Charles,  Jr. 

Apparatus  and  Method  for  the  Simultaneous  Measurement 
of  Passive  Membrane  Transport  Properties 
AL  H.  Friedman  and  R.  A.  Meyer 

-Sensor  for  Measuring  Pitch  and  Yaw  Angular  Orientation 
Between  a Spacecraft  and  an  Attached  Magnetometer 
F.  V.  Scheniel  and  A.  Finkel 


Sensor  for  Measuring  the  Relative  Twist  Between  a Space- 
craft and  an  Attached  Magnetometer  by  Means  of  a 
Dihedral  Mirror 

F.  IF'.  Schenk  el  and  A.  Finkel 

Continuous  Loop  Flywheel 
D.  IF'.  Rabenhorst 

Apparatus  for  Identifying  Missile  Codes  Without  Internal 
Clock  Synchronization 

C.  Phtlippides  and  IF'.  H.  Zinger 
Missile  Identification  and  Data  Entry  Apparatus 

C.  Phtlippide  i 
Rotor  Balancing  Apparatus 

D.  IF'.  Rabenhortt 

Kite  String  Wound  Flywheel  Rotor 
D.  IF'.  Rjbenhont 

Laser  Interferometric  System 

/.  C.  Murphy  jnd  R.  C.  Cole 

Logarithmic  Time  Variable  Gain  Circuit 
/.  H Knck 

No  Pump  Flow  System 
/.  H.  Pjnesci 

Submarine  Launched  Expendable  Flashing  Buoy 

G.  S.  Billing i 

Laser  Beam  Folding  Device 
G.  S.  Bitting s 

Communication  Line  Switch 
A.  IF".  Currano 

Natural  Ice  Collection  Device  and  System 
IF”.  R.  Powell 

Spoked  Multi-Ring  Filament  Rotor 
D.  IF'.  Rjbenhont 
Antenna  Alignment  indicator 
J.  R.  Butler 

Weight  Release  Mechanism 

R.  A.  Mottey  and  D.  C.  Wenslrand 

High  Density  Digital  Magnetic  Recording  Using  a (5,  6) 
Alternating  Disparity  Block  Code 
AL  Daridtun  and  /.  L.  Machamer 
Differential  Screw  Locking  Mechanism  for  Finely  Adjustable 
Optical  Mount 
R.  H.  Lapp 

Cartridge  Tape  Formatter  and  Interface 
T.  G.  Boland 

Long-life  Source  for  Thermal  Vacuum  Deposition 
/.  L.  Abita  and  J.  G.  Bebee 


PATENT  APPLICATIONS 

Human  Tissue  Stimulation  Electrode  Structure 
R E.  Fischell 

Epidural  Lead  Electrode  and  Insertion  Needle 
R.  E.  Fischell  and  IF'.  R.  Powell 

Low  Axial  Force  Servo  Valve  Spool 
IF'.  Seamone 

Digital  Beamsteering  for  a Parametric  Scanning  Sonar  System 
/.  IF'.  Follin  and  R.  E.  Miller 

Pressure  Responsive  Switch 

R.  T.  Cusick  and  IF'.  J.  Pleagle 

Ionic  Air  Speed  Indicator 
Af.  L.  Hill 

Rechargeable  Body  Tissue  Stimulator  with  Back-Up  Battery 
and  Pulse  Generator 
R.  E.  Fischell 


Solar  Panel  Deployment  Mechanism 
L.  E.  Stillman  and  T.  B.  Coughlin 
Photoacoustic  Radiometer 

/.  C.  Murphy  and  L.  C.  Aamodt 
Single-Axis  Disturbance  Compensation  System 

R.  E.  Fischell,  A.  C.  Sadilek,  F.  F.  Mobley,  and  G.  H. 
Fountain 

Method  for  Correcting  Navigation  Errors  Due  to  Water 
Currents 

R.  H.  Bauer  and  Af.  C.  Osborne 
A Clock  Invariant  Synchronization  Technique  for  Random 
Binary  Sequence 
C.  Philippidet 


PATENTS  ISSUED 

No.  4,000,665 

Woven  Filament  Rotor  Structures 
D.  IF'.  Rabenhorst 

No.  4,001,728 

Digital  Method  of  Pulse  Width  Modulation 
IF'.  Schneider 
No.  4,005,415 

Automated  Radar  Data  Processing  System 
A.  Kossiakoff  and  ].  R.  Austin 

No.  4,011,559 

Universal  Binary  Code  Converter 

D.  L.  Sharp  and  E.  A.  Frekko 

No.  4,011,564 

Phase  Modulated  Monopulse  System 
/.  F.  Gulick 
No.  4,016,560 

Fractional  Binary  to  Decimal  Converter 
R.  IF'.  Fowler 

No.  4,019.032 

X-Y  to  Range-Bearing  Converter 

E.  J.  McDevitt 
No.  4,020,714 

Filament  Connected  Rim  Rotor 
D.  IF'.  Rabenhorst 

No.  4,023,437 

Filament  Rotor  Having  Elastic  Sheaths  Covering  the  Fila- 
mentary Elements  of  the  Structure 

D.  IF'.  Rabenhorst 

No.  4.025,917 

Simplified  Time  Code  Reader  with  Digital  PDM  Decoder 
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